SN-230 



TRANSPIRATION AND FILM COOLING 


BOUNDARY LAYER COMPUTER PROGRAM 


J 




VOLUME n 


Computer Program and User's Manual 


FINAL REPORT- CONTRACT NAS7-79I 


• IIINF 1971 

(NASA-CR- 126728) TRANSPIRATION AND FILM 
COOLING BOUNDARY LAYER COMPUTER PROGRAM. 

VOLUME 2: COMPUTER PROGRAM AND USER'S 

MANUAL Final Report R.J. Gloss (Dynamic 
Science) Jun. 1971 1U7 p CSCL 2CM G3/12 

PREPARED BY- ROGER J. GLOSS 


N72-253C 1 


Unclas 

15425 



PREPARED FOR : NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 


DYNAMIC SCIENCE 

A Division of Msrshatt Industries 



2400 Michel son Drive, Irvine, California 92664 (714) 833-2670 


r ** ’i* **- • f\>i i. J ! ' {' I ^,'*.1'!. 


i 


*\ 

t 


SN-230 


— 

TRANSPIRATION AND FILM COOLING 
BOUNDARY LAYER COMPUTER PROGRAM 




J 


VOLUME n 

Computer Program and User's Manual 


FINAL REPORT- CONTRACT NAS7-79I 
JUNE 1971 

PREPARED BY : ROGER J. GLOSS 

PREPARED FOR : NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 


DYNAMIC SCIENCE 

A Division of Msraho* Indvs trios 




2400 Miche*»on Drive, Irvine, California 92664 (714) 833-2670 




FOREWORD 


The present report is part of a two volume set which describes a turbulent 
boundary layer analysis and computer program that includes the effects of wall 
cooling and equilibrium chemistry. Volume I contains the analytical basis for 
the computer program and a discussion of the results obtained to date: Volume II 
of the set describes the computer program and serves as a user's manual. 

This investigation is entitled TRANSPIRATION AND FILM COOLING 
BOUNDARY LAYER COMPUTER PROGRAM. The two volumes are additionally 
subtitled as follows: 

Volume I - Numerical Solution of the Turbulent Boundary Layer Equations 
With Equilibrium Chemistry - by Jay N. Levine 

Volume II - Computer Program and User’s Manual - by Roger J. Gloss 

Volumes I and II, and the computer program, have been distributed 
according to the attached distribution lists. Additional copies of the two 
volumes and the computer program (UNIVAC 1108 and IBM 360 versions) may 
he obtained from T. Reedy, CPIA, APL/JHU, 8621 Georgia Avenue, Silver 
Spring, Maryland 20910. 

This investigation was conducted for the Jet Propulsion Laboratory, 
National Aeronautics and Space Administ ation under contract No. NAS7-791 with 
Walter B. Powell as technical monitor. Jay N. Levine of Dynamic Science was 
Program Manager. 
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NOMENCLATURE 


A 



L 

F 

G 

G\ G" 
H 
h 
j 

K 

k 


L 

Le 


Le T 

l 

M 


m w 

m,n 

n 

P 

^amb 

Pr 

Pr»p 

q w 


w 


Van Driest damping constant in eddy viscosity formula 
coefficients in the difference equations 

value of y at the initial station 
e 

specific heat 

mass fraction of the ith species 
local shear stress coefficient 

right hand side of difference equations 

eddy diffusivity 

pG' r y/C 

G* /£ 8 Re , also thrust 

function used to transform y 

1st and 2nd derivatives of the G function 

total enthalpy, also shape factor = 6 * / 0 

static enthalpy 

equals 1 for axisymmetric flow, 0 for planar flow 
acceleration parameter 

thermal conductivity, also common ratio in geometric 
stepsize progressions 

reference length 

molecular Lewis number 

turbulent Lewis number 

mixing length 

molecular weight 

surface mass transfer rate = p w v w 

mesh indicies 
parameter in G function 
pressure 

ambient pressure at nozzle exit 
Prandtl number 
turbulent Prandtl number 

local heat transfer rate 
nozzle throat radius 
local radius of the wall 


/I <w vf«h n 




R 

Re 

S 

s 

St. 

T 

u,v 

+ 

u 

V 

+ 

v w 

X 

y 

■+ 

y 

Y 


Universal gas constant 
Reynolds number 

- entropy 

- wetted length along wall 

- Stanton number 
temperature 

velocity in s and y directions, respectively 

- dimensionless velocity ratio, u/ u* 

w 

1/2 

friction velocity, (t w /p w ) 

velocity in y direction 
dimensionless velocity ratio, v w /u^ 

axial distance 
normal coordinate 

dimensionless normal coordinate, y u^/j7 
nondimensional y coordinate =y/L£ 


trans 


constant in G function 

- mass fraction of the ith element 

mass fraction of the ith element in the transpiration coolant 

- constant in G function 
ratio of specific heats 

- boundary layer thickness 
displacement thickness 

- incompressible displacement thickness 

kinetmatic eddy viscosity 
function used to scale y 
d£ / ds 

- momentum thickness 

- eddy conductivity 


iv 
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molecular viscosity 
kinematic viscosity 
density 

used to represent various transport property terms 
shear stress 


Subscripts 

e 

i 

r 

w 

y 

yy 


at the edge of the boundary layer 
pertaining to the ith species 
evaluated at reference conditions 
at the wall 

partial derivative with respect to y 
2nd partial derivative with respect to y 


Superscripts 

i - equals 1 for axisymmetric flow, 0 for two dimensional flow 

1 - denotes a fluctuating term 

denotes a dimensional variable and/or a time averaged 
quantity 


v 


1. INTRODUCTION 


This report describes a new finite difference turbulent boundary layer 
computer program which allows for mass transfer wall cooling and equilibrium 
chemistry effects. Volume I of this set describes the analytical foundation of 
the program and contains a discussion of many of the solutions achieved to date. 

The program is capable of calculating laminar or turbulent boundary 
layer solutions for an arbitrary ideal gas or an equilibrium hydrogen oxygen 
system. Either two dimensional or axisymmetric geometric configurations may 
be considered. 

The equations are solved, in nondimensionalized physical coordinates, 
using the implicit Crank-Nicolson technique. The finite difference forms of the 
conservation of mass, momentum, total enthalpy and elements equations are 
linearized and uncoupled, thereby generating easily solvable tridiagonal sets 
of algebraic equations. 

The computer program is written in A.S.A. standard FORTRAN IV, and versions 
are now operational on UNIVAC 1108 and IBM 360 computers. 

This volume is designed to serve as a detailed description of the computer 
program, as well as a program user's manual. Detailed descriptions of all 
boundary layer subroutines are included, as well as a section defining all program 
symbols of principal importance. Instructions are then given for preparing card 
input to the program and for interpreting the printed output. Finally, two sample 
cases are included to illustrate the use of the program. 


V 


2. PROGRAM DESCRIPTION 


The Mass Addition Boundary Layer Program is actually composed of two 
distinct modules or collections of subroutines: (1) a boundary layer analysis 
which contains all subroutines necessary for a perfect gas case, and (2) the One 
Dimensional Equilibrium Program (0DE) , modified to operate in a subroutine mode 
and to solve, in particular, the hydrogen-oxygen chemical system. These two 
modules interface, in a hydrogen-oxygen equilibrium case, vi<? a single subroutine 
H00DE. 



In this section, all program subroutines are described briefly. Detailed 
documentation then follows for all boundary layer subroutines and for subroutine 
HOODE. Although the equilibrium subroutines are not discussed in detail, sub- 
routines CPHS, EQLBRM, MATRIX, MGAUSD, REACT, and SEARCH are described 
in the Two Dimensional Kinetic Reference Program (TDK) document (Reference 1) . 




Boundary Laver Subroutines 


AD’ *>T - 
BN JJND - 
C0NTNU - 
DEBUG - 
DUMPIT - 

EDDY - 
ELEMTS - 

ENERGY - 
EXECUT - 
GFUNC - 

H00DE - 

IG0DE - 

HERAT - 

LCURV - 
M0MNTM - 
NL0UT - 
PARAMS - 
PH0ENX - 

PRINT - 

PR0FIL - 


Adds mesh point at edge of boundary layer. 

Calculates wall and edge conditions at forward station. 

Integrates continuity equation to obtain pv. 

Error subroutine which indicates origin of program-detecteu error. 

Dumps matrix coefficients of a given difference equation at 
selected stations. 

Calculates turbulent transport properties ( , Pr^, Le^. 

Solves element conservation equation for mass fraction of 
hydrogen <* . 

Solves energy equation for enthalpies H and h. 

Execution control subroutine. 

Generates boundary layer stretching function G and calculates 
boundary layer mesh. 

Interface subroutine which calls One Dimensional Equilibrium 
subroutines to calculate equilibrium properties for hydrogen-oxygen 
system. 

Calculates thermodynamic and laminar transport properties for 
perfect gas case. 

Updates average properties and calculates auxiliary quantities 
for a given iteration on boundary layer equations. 

Linear interpolation subroutine. 

Solves momentum equation for velocity u. 

Prints out program input. 

Calculates gross boundary layer parameters at a given station. 

Interpolates for missing values in a table in which every nth value 
is known. 

Stores items in summary table at each station and prints profiles 
at a given station if requested. 

Generates initial dependent variable profiles across boundary layer 
at start of case. 


RDTAPE - 
STEP - 

SUMTAB - 
TABLES - 
TFCBL - 

TRIM- 

XNTERP - 


Reads restart data from tape for a restart case. 

Determines stepsize and calculates contour properties at 
forward station. 

Writes summary table at end of case. 

Normalizes 'nput tables and initializes wall and edge conditions. 

Main prog? which reads input data and controls calculations 
at initial siacion. 

Solves system of linear equations whose matrix of coefficients 
is tridiagonal. 

Averaged parabolic interpolation subroutine. 

Evaluates boundary layer normalization function £ every iteration. 


ZFUNC - 


One -Dimensional Equllibrlun. Subroutines 


ANSWER - 

BL0CK DATA 
CPHS - 

CPSPEC - 

EQLBRM - 
HPGALC - 

MATRIX. - 

MGAUSD - 

0DE - 
REACT - 
SEARCH - 
SPCALC - 

TPCALC - 
VISCX - 


* 


Calculates total number of moles, sound speed, and density 
based on equilibrium solution. 

Stores atomic symbols, weights and valences. 

Calculates species enthalpies and total specific heat from 
temperature curve fits. 

Calculates specific heat for each species from temperature 
curve fits. 

Calculates equilibrium composition and properties. 

Performs a series of up to thirteen enthalpy-pressure equilibrium 
calculations. 

Constructs the iteration equations required for the equilibrium 
solution. 

Solves linear set of up to twenty simultaneous equations using 
modified Gaussian elimination. 

Initialization subroutine for equilibrium portion of program. 

Reads and processes REACTANTS data. 

Reads, processes, and stores THERM0 data. 

Performs a series of up to thirteen entropy-pressure equilibrium 
calculations. 

Performs a single temperature-pressure equilibrium calculation. 

Calculates viscosity and Prandtl number for hydrogen-oxygen 
system using mixture formulas. 



SUBROUTINE ADOPT 


This subroutine adds another mesh point at the edge of the boundary 
layer and prepares for recalculation of the difference equation coefficients for 
the last two equations in the matrix solution. ADDPT may be called from the 
following subroutines under the conditions specified: 

(1) from M0MNTM, if, after solving the momentum equation, 


( 2 ) 


u 


du 

dy 


EPSLN1 


from ENERGY, if, after solving 


(input) ; 

the energy equation, 


(3) 


H 


dH 

dy 


> EPSLN2 


from ELEMTS, if, after solving 


(input) ; 

the element conservation equation. 



da 

dy 


e 


> EPSLN3 (input). 


Subroutine ADDPT is called as follows: 


CALL ADDPT (IFLAG) 


where 


IFLAG is a flag (integer) indicating the subroutine 

from which ADDPT is called: 

= 1, if ADDPT is called from subroutine 
M0MNTM (i.e. the momentum equation 
does not satisfy the edge criterion); 

= 2, if ADDPT is called from subroutine 
ENERGY (i.e. the energy equation does 
not satisfy the adge criterion); 

= Z, if ADDPT is called from subroutine 
ELEMTS (i.e. the element equation does 
not satisfy the edge criterion) . 


* 


Subroutine ADD PT first increments NY, the number of points in the 
boundary layer by one. Accordingly, 

NY = NY + 1 

NY1 = NY - 1 

NY 2 = NY - 2 

During initialization, the arrays for y, y, G', and G" , were already 
evaluated for the maximum possible number of points in the boundary layer, 
so that these quantities are already stored for the newly added point. All 
other properties which vary across the boundary layer, however, must be ex- 
tended to the new point. For all properties except those noted here, this 
amounts to setting the value at the new point equal to the value at the old 
edge point. For example, for velocity, 

U (NY, J) = U (NY1 , J) for J= 1,2,3. 


Note that NY1 is the subscript for the point which was previously the edge of 
the boundary layer but is now the first interior point. Note also that back, 
average, and forward values are extended. 


The quantity RH0V (pv) is not constant at the edge but is integrated 
to obtain the value at the new point: 


> ;V NY ' V NY1 + Ay d dy^ 


NY1 


where 


d (py) 

dy 


NY1 


is parabolic derivative at the old edge point. 


The auxiliary quantities E and F are then evaluated at the new point NY. 

Following return from subroutine ADDPT, the equation which necessitated 
the addition of the point will be solved again and the solution will be checked 
to see if it satisfies the edge criterion. Before doing this, however, the difference 
equation coefficients for ‘he points near the edge must be recalculated due to 
changes resulting from addition of the new pomt. The following discussion refers 
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to the diagram below. 




Before point addition After point addition 


NY2 0 


o 

0, 

4 

Coefficients here un- 
affected by point addition 

NY1 0 

Coefficients hare were 
mo< .fied by boundary 
condition 

0 

NY2^ 

n 

Calculate o. (1) and 
ctj (2) here 

NY 0 


o 

NY1 

n 




o 

NY 



After point addition, all coefficients for points interior to NY2 n are 

unchanged and need not be recalculated, at NY2 the coefficients were modified 

n 

to account for the boundary condition at NY q , but NY2^ is now merely an interior 
point whose coefficients depend on other interior points. Thus the coefficients 
must be recalculated at NY2 n by the calling routine, and to do this, the (1) 
and Oj (2) must be calculated in ADDPT at NY2 n - 1 and NY2 r respectively. 
Accordingly, if IFLAG = 1, ADDPT calculates a ^ at these points for the momentum 
equation; if IFLAG = 2, ADDPT calculates Og , and a ^ for the energy equation; 
if IFLAG = 3, ADDPT- calculates u 5 for the element equation. 


Finally, before returning to the calling routine, subroutine ADDPT prints 


the following message: 


POINT WAS ADDED TO BOUNDARY LAYER IN SUBROUTINE 
STATION AND ITERATION . 


momntm] 

ENERGY } AT 
ELEMTS J 
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SUBROUTINE BNDCND 


w 


This subroutine obtains values at the forward station (m+1) of quantities 

necessary for the difference equation boundary conditions. Previous forward 
dP_ du 


values of 


ds 


■jj- , and m w , are first stored as the back values for the new 


station. The following quantities are then obtained: 
Name 
D TT ED Q N 


DUEDS 


UEDGE 


dv^ 

ds 


u 


'm+1 


Quantity 

How Obtained 

du e 


du 

e 


ds 

m+1 

ds 

m+1 ^ 



du e 1 

where : 


cos ft 


m+1 w m+l 

is obtained 
by interpolation using LCURV. 

Average of m and m+1. 


Euler integration: 


u 


u 


'm+1 


+ A s 


m 


ds 


m 


SMDWN 


SMDW 

TWALL 

PEDGEB 

DPEDSN 


DPEDS 


m 


w 


m+1 


m 

w 

T 

w 

P 

e 

dP 

i 

"ds 


e 


-fe 

ds 


w 


m+1 u 

m w " n 'u"r ' Whete 

_ w m+l u r ^ m+1 

m is obtained by linear 
m+1 

interpolation using LCURV. 
Average of m and m+1 . 

Linear interpolation using LCURV. 

Interpolation using XNTERP. 



l!e 

i 

dP 

e 

m+1 

ds 

m+1 ^r ^ B 

"Sc 


cos 


dx 


where 

m+1 

using XNTERP. 


m+1 
is obtained 


% 


Average of m and m+1. 



m+1 



SHEDGE 


SHWALL 

HWALL 



h_ = H - u V2, whore H 
e e e e 

was calculated initially and is 
constant. 

For an equilibrium case, subroutine 
H00DE is called to obtain h w from 

a temperature^pressure calculation 
using T w and P g . For a perfect gas 

case, subroutine IG0DE is called 

with T , and P to obtain h . Then 
we w 



In both cases H = h . 

w w 


SUBROUTINE CflNTNU 


$ 


Subroutine C0NTNU integrates the continuity equation to determine 
pv vs y at each mesh point across the boundary layer. 


’"he continuity equation is solved at m + 1/2, that is, at a distance 
halfway between the back station m and the forward station m + 1. The derivatives 
and coefficients are evaluated at m + 1/2, n - 1/2. (See diagram.) 


Thus, for example, 

Sn+1/2 

Vl/2 = 



y (r w +r » > 

m m+1 


s 


We also have 


° u m+l/2, n-1/2 

dpu _ 1 

as 2 as 

&Pu _ 1 

ay 2 Ay 

apy _ l 

57 Ay 


“ 4" ^ DU m,n + P u m,n-1 + P u m+l,n + ® u m+l, 
^ ° u m+l,n + ° u m+l, n-1 ” ^ u m,n ° u m,n-P 
^ Du m+l,n + ^ u m,n " Pu m+1, n-1 ” pu m,n-P 
(pV m+l/2,n " P v m+l/2, n-1* 


n- 


1 


) 
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With the properties and denva es established in this manner, the continuity 
equation is integrated as follows: 


ov rr.+l / 2 , 1 
DV ir.+l/2,n 


rr. 


w 


m+I 2 


at the wall; 




rr.+i '2 .i-l 


1 i? ou 

G* i //j 
n-1/2 


Pu m+l/2. n-1/2 
G n-l/2 r » m+1/2 



m+l/2 


( -^) 

* y m+l/2. n-1/2 


SUBROUTINE DEBUG 


This subroutine is called when the program detects an errm; it prints a 
message indicating where and when the.* error occurred, causes the summary table 
up to this point to be printed, and th> ,1 terminates program execution by calling 
EXIT. Subroutine DEBUG is called as follows: 

CALL DEBUG (SNAME) 

where SNAME is a single cel! containing the Hollerith repio entation of the name 
of the calling routine, where the error was encountered. 

Subroutine DEBUG prints the following message: 

EXIT CALLED FROM SUBROUTINE [ SNAME ] AT STATION 

AND ITERATION . 

Subroutine DEBUG then calls subroutine SUMTAB to print the summary 
table (See description of subroutine SUMTAB.) and finally calls EXIT. 



0 


SUBROUTINE DU M PIT 

This subroutine optionally prints program debugging information over a 
specified range of stations for a given case. Either the momentum equation or 
the energy equation matrix coefficients, as well as the values u, H, h, and pv 
are printed for the first five and last five mesh points in the boundary layer. The 
output numbers are not labelled; however, they may be identified by comparison 
with the WRITE statements in subroutine DUMPIT which produce them. 

See the section on Additional Program Features for instructions for imple- 
menting this debug capability. 
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SUBROUTINE EDDY 


$■ 


SUBROUTINE EDDY 

This subroutine is called every iteration to calculate the eddy viscosity 
e, turbulent Prandtl number Pr,p, and turbulent Lewis number Le^, at each mesh 
point, based on the latest iteration values of all the other boundary layer properties, 

The eddy viscosity calculation begins at the wall point and proceeds toward 
the edge. A quantity and a quantity c Q are calculated at each mesh point, and 
€ = at each point until e i ^ for some point. From that point on out to the edge 
of the boundary layer, only ( Q is calculated, and * = « 0 . 

The equation for e i is 


;i ■ 


e i = ^1 R e r C G' y 2 


where 


9 u 
9y 


(l-e 0 ^ 


Kj = 0.4, 


a = 


Re CY N 
r * 

26 u 


P4-G’ 


w w 


R e r C 


9u_ 

by 


1/2 


w 


In the expression for a , N is calculated in one of two ways. If m w ^0, 


N 


4 - 


ils 

as 


qi b u 
w Sy | w w 


( 1 - efy + e f 


m 


1/2 


where 


0 = 


11.8 fn 

w 


Re r 


1 1/2 


r . a_u_ 

°W aF 1 


w 


If m w = 0, then N is constant across the boundary layer, and 

- Re C 3 ^ 1/2 ' S 


N = < 1 + 


11.8 


a p. 


P w l/2 as 


l_Mw (G w "tyU j 
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e Q is calculated as follows: 

*o =K 2 u e Re r 6 inc [ 1+5 * 5 ("|“ ) ] 

where Kg = .0168, 6 is the y - location at which u/u g = .995, and 

0 e 

The correspondence between the quantities defined above and the local subroutine 
variables is given below: 

Physical Quantity Program Variable 


«i 

EPSI 

c o 

EPS# 

K 1 

BK1 

K 2 

BK2 

a 

BRKT 

e^ 

TERM 

N 

BN 

6 

DELTA 

«fnc 

DLSINC 


The eddy viscosity calculation begins by finding 6, where r.Ai = .995. 

6 £ is then evaluated using trapezoidal integration. The eddy visclsity loop 
starts near the wall with the auxiliary multipliers FIN = 1, and F#UT “0. Both 
EPSI and EPS# are evaluated and compared. If EPSI <EPS#, FIN and F#UT remain 
unchanged; however, if EPSI >. EPS#, the routine sets FIN = 0. and F#UT =1. 
through the rest of the boundary layer and EPSI is no longer calculated. The eddy 
viscosity at every mesh point can thus be set as follows: 

EPS(I, JN) = FIN*EPSI + F0UT*EPS# 

Finally, in order to assure a smooth eddy viscosity profile across the boundary 
layer, e at each interior point is recalculated as the average of the values cal- 
culated for that point and the four neighboring points, two on each side. 
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At the wall, the turbulent Prandtl number is given by 


Pr„ 



34 

26 



Elsewhere in the boundary layer. 



1 - exp 


1 - exp a 

ttt 


26 


r) 


where a is as defined previously for 

Presently, the turbulent Lewis number Le^, is set to unity throughout the 
boundary layer. 



SUBROUTINE ELEMTS 


Subroutine ELEMTS solves the element conservation equation for a hydrogen- 
oxygen equilibrium case to obtain a , the mass fraction of hydrogen, at each mesh 
point in the boundary layer. The mass fraction of oxygen is then equal to 1 - a at 
each mesh point. In normalized, finite-differenced form the element conservation 
equation is as follows: 


A i a. 


1 °m+l,n-l + A 2 a m+l,n T A 3 tt m+l,n+l 


+ A~ «. 


= B 


where 


B 


Eu _ 

m. n _ P 

~ 4 


A S 


= - A, - 


FG* o c 
FG* (r c 


+ -jj- [ ( G' 
Ay 2 


QH J 

°sy + 'g r<, s ) TK7 ‘ ay 




1 £ y* 


pUOf m,n p vG* . F L.. . G" , 

iT* Ji T ®y „ T (g °5y + '5 r *s> a y m „ 

m # n L * *m,n 

+ G,(T 5*yy I + f a y u m+l,n 

r7 m,nJ *m,n 


Here 

PG'C’y p _ G' 1 
C ' I s ” Re r ' 

= MLe . Ie T 
ff 5 Pr p€ Pr T * 

The velocity u in the expression for B has already been obtained by 

subroutine MJ0MNTM. 

The hydrogen element equation boundary conditions are as follows. At 
the edge of the boundary layer. 


Of = Cf e * 
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' w * w.- -v, ■* r~s aa^-sais-s HL 







At the wall. 




w 


- A 0 ( “w' a tr»ns ) - 


where is the mass fraction of hydrogen in the transpiring gas (in 

troiis 

°" G ‘«w Le w ^ 


a d 


Using a three-point derivative formula for d - Qi 
becomes 


, the wall boundary condition 


w 


a 


2a y * ” 3a l + 4 “2 " a 3 ) A 0 ® 1 " A 0 “trans 

4a 2 -a 3 + 2 ^ y A 0 «trans 


2Ay A q +3 


In subroutine ELEMTS, the correspondence between the program subscripts 
and the physical points in the boundary layer is the same as in subroutine M0MNTM. 

The push-down storage technique explained in the description of subroutine 
M0MNTM is applied to the auxiliary array SIG5 (oj.) in subroutine ELEMTS, 

Applying the differenced element equation at each interior point in the 
boundary layer results in a tridiagonal matrix of NY - 2 equations as in subroutine 
MOMNTM. To eliminate aj and as unknowns in the system of equations, the 
boundary conditions are applied. For the first equation in the matrix. 


A( 1 , 1) ttl +A(1,2) a 2 +A(1,3) « 3 = B(1) 

Substituting the expression for a j obtained above from the wall boundary condition 
yields 

|a( 1,2) + 2 Ay aJ+3 a 2 + [ A(1,3) " 2 Ay Aq+ 3 a 3 = 


l 


B(l) - 


A(l,l) 2 A yA 0 0^3 

2AyA Q + 3 




i i 
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Therefore we can set 


§ 


A( 1 , 2) 

A(l,2) 

+ «U4L 

2AyA Q +3 

A( 1 , 3) 

A(l,3) 

- *(i.D 

2AyA Q +3 

B(1) 

B( 1) - 

Ml.l) 24yA oatrans 

2Ay A 0 + 3 

A( 1 , 1) 

0 . 



For the last equation in the matrix 

A(NY2,1) a NY _ 2 + A(NY2, 2) + A(NY2,3) a m = B(NY2) 

The edge boundary condition is a m - a e so that we can set 

B(NY2) = B(NY2) - A(NY2, 3) WEDGE 

A(NY2,3) = 0. 


Subroutine TRIM is now called to solve for a 2# tt;J , 
tt NY are set according to the boundary conditions. 


• • • •» ft nY“ 1 ant ^ ft ^ and 
Subroutine ELEMTS finally tests 


whether the hydrogen mass fraction profile obtained from the matrix solution 
asymptotically approaches a e at the edge of the boundary layer. If 


1 dot I 

of e ay | 

e 


*:EPSLN3 (input). 


control returns to the calling routine; otherwise another mesh point is added at 

the edge. The mechanism of mesh point addition in subroutine ELEMTS is the same 
as in subroutine M0MNTM. 



SUBROUTINE ENERGY 




Subroutine ENERGY solves the energy equation to obtain the total enthalpy H 
and the enthalpy h ^t each mesh point in the boundary layer. In normalized, finite- 
differenced form the energy equation is as follows: 

A l H m+l,n-l 4A 2 H m+l,n' fA 3 H m+l,n+l = B ' ' C 1 u m+l,n-l “ C 2 u m+l,n “ C 3 u m+I,n+1* 

where the velocities u m+ ^ on the right-hand side have been obtained already by 
subroutine M0MNTM, and 
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Here, Z is the summation over all species, and 
i 

E = PC C-y F , g: j_ 

C ' F 

°2* nT + + P< (1 -P^>' 

°4 = ft - (Le ‘ U + (le T * 11 

The energy equation boundary conditions are 

H = H at the wall; 
w 

H = H at the edge, 
e 

In subroutine ENERGY, the correspondence between the program subscripts 
and the physical points in the boundary layer is the same as in subroutine M0MNTM. 

The push-down storage technique explained in the description of subroutine 
M0MNTM is applied to the auxiliary arrays SIG2 ( CTg ) » SIG3 (a 3 ) and SIG4 in 
subroutine ENERGY • 

Applying the differenced energy equation at each interior point in the 
boundary layer results in a tridiagonal matrix of NY-2 equations as in subroutine 
M0MNTM. To eliminate Hj and Hj^y as unknowns in the system of equations, the 
boundary conditions are applied. For the first equation in the matrix 

A( 1 , 1) Hj + A( 1 , 2) H 2 + A(l, 3) H 3 = B(l) 

(B(l) includes u m+ j terms as well as B‘.) 

The boundary condition at the wall is H = H , so that 

X w 

A( 1 , 2) H 2 +A(l,3) H 3 =B(1) -A(l,l) H w 
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Therefore, in the coefficient matrix we can set 

B(l) = B(l) - A(1,1)*HWALL 
A( 1 , 1) =0. 

For the last equation in the matrix 

A(1MY2, 1) H ny _ 2 + A(NY2, 2) + A(NY2,3) = B(NY2) 

The edge boundary condition is = H g , so that 

A(NY2, 1) H ny _ 2 + A(NY2, 2) = B(NY2) - A(NY2, 3)*HEDGE 

We can then set 

B(NY2) = B(NY2) -A(NY2,3) *HEDGE 
A(NY2, 3) =0. 

Subroutine TRIM is now called to solve for H 2 , H^, . . . h ny_i» and H i and 
Hny are set according to the boundary conditions. The enthalpy h at each mesh 
point is then calculated as follows: 

h =w 1 2 

m+l,n m+l,n 2 u m+l,n 


Subroutine ENERGY finally tests whether the total enthalpy profile obtained 
from the matrix solution asymptotically approaches H at the edge of the boundary 
layer. If 


i aH 


H e 171 

e 


s EPSLN2 (input) , 


control returns to the calling routine; otherwise, another mesh point is added at 
the edge. The mechanism of mesh point addition in subroutine ENERGY is the 
same as in subroutine M0MNTM. 
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SUBROUTINE EXECUT 


Subroutine EXECUT is the execution control routine for the Mass Addition 
Boundary Layer Program. Once the main program has completed all input pro- 
cessing and initialization, subroutine EXECUT is called to solve for the boundary 
layer properties from X = XINIT to X = XMAX. Subroutine EXECUT tests the input 
flags to call the appropriate subroutines for a given problem (perfect gas or equili- 
brium, laminar or turbulent, compressible or incompressible) and checks for 
convergence of the iterations at a given station. Following is a flowchart showing 
the logic in subroutine EXECUT in detail. 



BEGIN NEW 

STATION 

CALCULATION 


BEGIN 

ITERATION 

LOOP 



MADFLG = 0 



zn 



► 

ISTATN=ISTATN+ 1 
ITER=0 







CALL STEP 

to determine new r . 0 . tfi. Ax 
w w 




^new 

= C old + AsC 

^ave 

- ^old + ^»nevP^ 


CALL BNDCND 
to evaluate wall and edge conditions 
at forward station 




NO 


CALL H00DE 

to calculate h w = H w based 
on latest O/F at wail 



y 


CALL ITERAT 
to update average properties 
and auxiliary quantities for 
difference equations 




CALL M0MNTM 
to solve momentum equation 
for velocity u. 


T 
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<Qncomp > o> 
|no 





CALL ENERGY 
to solve energy equation 
for enthalpies H and h. 



CALL ELEMTS 

to solve element equation for 
mass fraction of hydrogen a 




CALL H00DE \ 
to evaluate p, p, Pr # T 
c k and h^i k™lf 2| • • • $ 6y 


CALL XG0DE 
to evaluate p, p, Pr t T 



<• 


CALL ZFUNC 

u P° Jte c„„„ c 


new *ave 




LAMNR >dy 
*NO 


< r CALL EDDY \ 

to calculate c# Pr^# Lftj / 




1 / CALL C0NTNU \ 


I X t0 inte 9 r8te continuity y 


1 N^eq nation for pv / 

lr 


i 

*1 











GO TO BEGIN 
NEW STATION 
CALCULATION 


CALL PARAMS \ 

to calculate gross boundary 
layer parameters. 



CALL PRINT 
to print final station 


I 


CALL SUMT 
to print summa. Nable 


I 



f RETURN ^ 
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SUBROUTINE GFUNC 

Subroutine GFUNC is called during problem initialization to set up the 
relationship between the normal boundary layer coordinate y and the non-dimensional 
coordinate y, based on information specified in the program input, y and y are 
related as follows: 

y = G(y) where y = 

Here G is a general stretching function which is uniquely specified by input 
quantities; L is the input reference length BLREF; £ is the boundary layer normal- 
ization function. 

The form of G is as follows: 

y = G(y) = b| [log ((e-1) ( 
where b = y g , n is the input quantity 

G (0) = 0 and 

G'(0) = GP0 (input) . 

Subroutine GFUNC begins by calculating y g = 6 j/(L£) , where 6 j is the 
input initial boundary layer thickness DELTA.I, and £ is the initial value of the 
boundary layer normalization function. 

The following arrays are constructed for the maximum possible number of 
points across the boundary layer: 

Y - y at each mesh point; 

YTIL - y at each mesh point; 

(Then y = y L£(s).) 

BGP - G' (y); 

BGFP - G" (y). 


+ a) + 1) l 1/n - 

SN.T, and a and 0are computed such that 
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If SN3 = 0 (nominal value if not input) , then no stretching has been 
specified, and 


where NYI is the input initial number of 
mesh point intervals across the boundary 
layer; 

y - y everywhere; 

G*(y) = 1; 

G"(y) = 0. 



If SN3 ^0, then « and £must be calculated. To find a. a function F is 
defined as follows; 

G 0 (y ) = bQog ((e-1) + l)] 1/n 

- the function G when a = 0 and p - 0 . 

<M» = [1»9 ((e-l) f + l)] 1/n ‘ l 

0 a b (e-1) -£ + 1 

F (y) = Gq (y) - GP0 

Newton-Raphson iteration is applied to solve for y^ such that F (y^) =0. If the 
solution does not converge in 100 iterations, subroutine DEBUG is called to indicate 
the error in GFUNC and terminate execution. Otherwise, when y^is found. 



Correspondingly, 

0 = Clog ((e-1) a + 1 )3 1/n . 

Once aand p are calculated, the precise G-function is known and the 
arrays Y, YTIL, BGP, and BGPP are constructed as follows: 

6j 

y e = G * y e * where y e = TT 
y e 

Ay = -W 
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> 1 = 1,2 


,max. # points 


YTIL(I) = G" 1 ( yi ) 

BGP(I) = G' (y t ) 

BGPP(l) = G M (y t ) J 

Subroutine GFUNC also stores the initial values of the following boundary 
layer quantities and counters: 

YTEDGE = y 

T e 

YEDGE = y 
7 e 

DY = AY (constant throughout the program) 

YTZETA = = 1 , the point where -J- = £ 

YZETA = y^ corresponding to 

NY = NYI + 1 , the initial number of mesh points in the 

boundary layer and the index of the edge point. 

NYI = NY-1 
NY2 = NY-2 


Finally, if a stretching function G has been specified (SN3 ^0), the arrays 
Y, YTIL, BGP, and BGPP are printed out by subroutine GFUNC. 



SUBROUTINE H00DE 


% 


Subroutine H00DE provides the interface between the boundary layer sub- 
routines and the equilibrium subroutines. Whenever the equilibrium properties at 
a given mesh point or group of mesh points are required for a hydrogen-oxygen case, 
subroutine H00DE is called from the boundary layer routines with an argument which 
indicates the type of equilibrium solution to be performed (i.e, , isentropic expansion, 
temperature-pre ’recalculation, etc.). Subroutine H00DE in turn calls the appro- 
priate equilibrium subroutine to obtain the required solution. H00DE also converts 
quantities from units usod in the boundary layer routines (pressure in lb^/ft , 
temperature in °R, etc.) to equilibrium units (pressure in atm. , temperature in 0 K* 
etc.) and back again. 

The calling sequence is as follows: 

CALL H00DE (ICALL) 

where 

ICALL is an input flag equal to 1, 2, 3, or 4 indicating 

the type of equilibrium calculations to be performed. 

The effect of each of the possible values of ICALL is discussed in turn. 

ICALL =1 


This is the initialization call. If IDEAL=0 in the DATA namelist, subroutine 
H00DE is called from the main program with ICALL=1. Subroutine H00DE in turn 
calls subroutine 0DE, which initializes the equilibrium storage areas and causes 
the THERM (3 data and REACTANTS data to be read. H00DE also evaluates the 
necessary unit conversion constants for use in later calls. 

No equilibrium solution results from this call to H00DE. 

ICALL=2 


Subroutine H00DE is called during initialization with ICALL=2 to perform 
an isentropic expansion at the edge of the boundary layer. Using the input values 
of PEDGE, TEDGE, and AFEDGE, H00DE calls subroutine TPCALC to perform a 
temperature -pres sure equilibrium calculation. This establishes an entropy SO and 
an enthalpy h = SHEDGE. Then H , which is constant along the edge of the boundary 

W V 

layer, is given by 

H = HEDGE = SHEDGE + UEDGE * UEDGE/2. 
e 
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Then, holding entropy constant and proceeding through the input pressure 
table PETAB, subroutine H00DE calls subroutine SPCALC to do up to thirteen 
entropy -pres sure calculations at a time, obtaining a value of h = SHE from each 

V 

calculation. From these calculations, an edge velocity table UETAB is obtained; 
for a given point DC in the pressure table, 

UETAB (DO = SORT (2.*(HEDGE - SHE)) 


ICALL=3 


Subroutine H00DE is called every iteration to do a temperature-pressure 

calculation at the wall. Using TWALL, PEDGE, and an 0/F ratio based on the 

latest iterated value of o^, H00DE calls subroutine TPCALC to obtain a value of 

h = SHWALL. Then H = HWALL = SHWALL. 
w w 

ICALL=4 


Subroutine H00DE is called every iteration with ICALL=4 after the main 
dependent variables u, H, h, a, and pv have been obtained across the boundary 
layer. The purpose of this call is to obtain values of p, ji, Pr, T, and h^ and c^ 
for each species at every mesh point. Enthalpy-pressure equilibrium calculations 
are performed only at every lYEQth mesh point (including the wall and edge points) ; 
the equilibrium properties at the intermediate points are obtained by interpolation 
using subroutine PH0ENX. Subroutine H00DE calls subroutine HPCALC to perform 
thirteen enthalpy-pressure calculations at a time until the edge of the boundary layer 
is reached. Subroutine PH0ENX is then called for the arrays RH0, SMU, PRN0, and 
TEM (The latter two arrays sure named PR and T in the other boundary layer subroutines.) 
to interpolate for the missing values. The species mass fractions c^ and species 
enthalpies h^ are not interpolated when Le=l and Le^,=l, because the terms which 
contain them then drop out of the boundary layer equations. 
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Subroutine IG0DE performs the following oalcuiAtKMX*; 


C « *E_ 
P y-1 


where < * 4d?t 1*011 endHt* the input 
weight FM&LWT* 

where y U the input spec MU' heat i*tU' 
GAMMA* 


If h and P are given (ITPHP * 0) , 


T = -4- 

P 
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If T and P are given (ITPHP = 1) , 


$ 


h = C T 
P 

Finally p, fi, and Pr are calculated as follows: 

_ P 
P “rT" 

ft t 1/2 

u = 2.27x10’° TTft-JT 

P (^^98^6) 

Pr - PRI where PRI is the input Prandtl number. 
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SUBROUTINE ITERAT 


This subroutine prepares for an iteration on the boundary layer difference 
equations by calculating averages of the various properties at each mesh point 
across the boundary layer. For each property at a given mesh point, the average 
of the value at the previous station and the latest iterated value is computed and 
stored for use in calculating the difference equation coefficients. For example, 
for dens’ty at mesh point I, 

RH0 (I, JA) =0.5* (RH0 (I,J0) + RH0 (I, JN)) 


Here RH0 (I, JA) is the computed average, RH0 (I, J0) the value at the previous 
station, and RH0 (I, JN) the latest iterated value. At the beginning of a new 
station calculation, RH0 (I, J0) and RH0 (I, JN) will be equal so that the average 
value will initially be equal to the value at the previous station. 


Properties whose averages are calculated and stored by ITERAT are u, H, 

Qfi 1™ Of# h, p, p, Pr , Ee , T, f , Pr#p, and Le^« 

Subroutine ITERAT also calculates, at each mesh point across the boundary 
layer, the auxiliary quantities 


Z = a „ d F = f 


Re 


Finally, ITERAT calculates and stores the temporary quantities 
Oj, a 4# and at the first and second mesh points to initialize the push-down 
storage feature used by the difference equation solution subroutines M0MNTM, 
ENERGY, and ELEMTS. 



SUBROUTINE LCURV 


& 


Subroutine LCURV performs linear interpolation in a table which may contain 
discontinuities and saves the place in the table. LCURV is used to obtain the wall 
t3mperature T from a table of T versus x and the mass aduition rn from a table 

% W W W 

of rn versus x. 

W 

The calling sequence is as follows: 

CALL LCURV (X,XTAB, YTAb, NP f K, Y) 

where 


X is the argument. 

XTAB is the argument table and contains the independent 

variables stored in increasing order. When tabu- 
lating discontinuous functions, the argument is 
repeated at the point of discontinuity. 

YTAB is the function table, where YTAB(I) is the value 

of the dependent variable corresponding to the 
argument value in XTAB (I) . 

NP is the number of entries in XTAB and YTAB. 

IX is a pointer which, upon input, is the position 

at which the argument search is to begin. Upon 
return from LCURV, IX is set such that 

XTAB (DC- 1) < X <; XTAB (IX) . 

By using the output value of IX as the input valua 
tor the next call to LCURV, the place in the table 
is saved. If DC s 0 or IX s NP, the search will 
start at the beginning or end of the table, respec- 
tively. 

Y is the function value obtained by inter plation. 
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Subroutine LCURV defines an arithmetic 

ol: '-k-i- y k-i' v y k> = y k-i + lx ' x k-i 


statement function 0 as follows: 

, y k ~ y k-i 
x k _x k-i 


LCURV then merely searches the table to locate the interva’ containing the argument 
x and uses C to linearly interpolate between the two surrounding table points for 
the function value Y„ If X falls exactly on a discontinuity, Y is set to the correspond- 
ing function value on the high-x side of the discontinuity. 

If X falls outside the table, LCURV linearly extrapolates using the first 
two or last two table entries. 

If NP 1, then LCURV sets Y = YTAB(l). A constant table may thus be 
specified by entering only YTAB(l), 



SUBROUTINE M0MNTM 


i 

$ 


Subroutine Mj&MNTM solves the momentum equation to obtain the velocity 
u at each mesh point in the boundary layer. In normalized, finite-differenced 
form, the momentum equation is as follows: 


A 1 u m+l,n-l + A 2 u m+l,n + A 3 u m+i,n+l " B 


where 


Eu_ 


A , = fii. _ + FG,<T ' 


FG' ffl 


2 as 2 

FG'o 


A y‘ 


A 3 " A i"Ty T 


1 


u 


B = pu 


_F_ 

2 


m.n PvG 


du 


A s 


2” u y, 


e 


m,n 


+ P U “dT + 


G" 


( -p r Oi G'ai ) u + G' a, u 

U 1 »» y m r> 


m,n 


yy 


m, n 


Here 


_ PG'l"9 r. _ G 1 1 _ _ , L 

- 7 ' F ~ 7^ rT '®i-m + P€* 


The momentum equation boundary conditions are 

u = 0 it the wall; 
u = u g at the edge. 

The features of subroutine M0MNTM described below apply equally to 
the energy and element conservation equation solutions in subroutines ENERGY 
and ELEMTS, respectively. 



Subscripting 


£ 


In evaluating the coefficients of the difference equations, the correspondence 
between the program subscripts and the physical points in the boundary layer is 
as follows: 

s - direction 
JN = m+1 

J0 = m 

JA = m+1/2 (average of m and m+1) 

v - direction 
T+l = n+1 

I - n 

1-1 = n-1 


Push-Down Storage 

The auxiliary quantity unlike E and F, is not preca'culated and stored 
at each mesh point. Rather, a push-down storage technique is used. When the 
difference equation at a mesh point I is being evaluated, only the a^'s at I and 
the two surrounding points 1-1 and 1+1 a~e needed because 

a l - *1 
.. n+1 n-1 

a iy n 2 Ay 

The code is set up so that a, and <y. are already available as SIG 1(1) and 

n-1 *n 

SIG1(2), respectively. (Upon entry tc subroutine M0MNTM, for 1=2, SIG1(1) and 
SIG1(2) have been appropriately calculated in subroutine ITERAT.) Then, for eval- 
uating the coefficients at mesh point I, SIG 1(3) is first evaluated as 

SIG 1(3) = SMU(1+1 , JA) + RHO(I+l,JA)*EPS(I+l,JA) 
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All aj' • essary for calculating the momentum equation coefficients at mesh 
point ‘ 4.ow available. Be—.e going on to mesh point 1+1, the o^s are re- 
positioned in their temporary array a* follows: 

SIG1(T) = SIG1(2) 

SIG1(2) = SIG1(3) 

The procedure is then repeated at 1+1. 

Matrix Solqtiog 

The differenced momentum equations, when applied at each mesh point in 
the boundary layer, constitute a syste... of linear equations whose matrix of coeffi- 
cients is tridiagonal. Specifically, the matrix is constructed by applying the 
difference equation at mesh points 2, 3, 4, NY-2, NY-1, where point 1 

is the wall and point NY the edge of the boundary layer. Since the momentum 
equation is not applied at the wall point or edge point, there will be NY-2 rows in 
all in the matrix, and the elements A(I-1, 1), A(I-1,2), A(I-1,3), and B(I-l) will 
contain Aj, A 2 , andB, respectively, for mesh pomt I. 

The wall and edge boundary conditions are figured into the first and last 
equations in .he matrix, so that there will be as many unknowns as equations. For 
the first equacion in the matrix 

A(1 , 1) Uj + A(l, 2) u 2 + A(l, 3) u 3 =B(1) 

The boundary condition at the wall is Uj - 0, so that 

A(l,2) u 2 + A(l, 3) u 3 =B(1). 

Therefore, in the coefficient matrix, we can set 

A( 1 , 1) =0. 

Ton the last equation in the matrix 

A(NY2, 1) u NY _ 2 + A(NY2, 2) Uj^.j + A(NY2, 3) .j^ = B(NY2) 



The edge boundary condition us u^ = u e , so that 

A(NY2, 1) Uny _ 2 + A(NY2, 2) u NY _j = B(NY2) - A(NY2, 3) U 0 

We can then set 

B(NY2) = B(NY2) - A(NY2, 3) *UEDGE 
A(NY2, 3) =0 

Subroutine M0MNTM now calls subroutine TRIM to serve for u 2 , u^,^. . 
U NY-1 anc * sets u j = 0 and u^y = u g in accordance with the boundary conditions. 

Point Addition 


After the momentum equation solution, subroutine MJ&MNTM tests whether 
the velocity profile obtained asympototically approaches u g at the edge of the 
boundary layer. If 


u 


du 

ay 


sEPSLNl (input). 


control returns to the calling routine. Otherwise, subroutine ADDPT is called to 
add one rresh point at the edge «..f the boundary layer and set all ptoperties at the 
new point. Upon return to M0MNTM, the difference equations are applied at the 
last two interior mesh points from the new edge point, the edge boundary condition 
is reapplied, and subroutine TRIM is called to solve the new system of equations. 
This process is repeated until the edge test is satisfied; however, if more than 
five points must be added or the maximum possible number of points in the boundary 
layer is reached, subroutine DEBUG is called to print a message indicating the 
error in M0MNTM and terminate execution. 
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SUBROUTINE NL0UT 


This routine prints the combined data from the DATA and TDKINP namelists 
in a format which is organizes for easy reference by the user. All single-cell 
variables are printed first in their various categories, followed by the tabular 
data. The output from subroutine NL0UT is in genral self-explanatory; however, 
the various input catagories are listed below in the order in which they appear, 
with explanatory comments where appropriate. 

Single-Cell Variables 

All single-cell variables appear on a single page, each in the format 
name = value. 

Where no value was input by the user, the nominal or default value is printed. 
Thrse variables are organized i the following categories: 

(1) Flags and options. 

Possible input values and their meanings are given in the output 
for each flag. 

(2) Problem limits and initial values. 

(3) Reference quantities. 

(4) Input nor. ^ ilzation factors. 

(5) Edge quantities. 

(6) Constants. 

(7) Convergence and edge criteria. 

(8) Counters. 


Tabu l 


p or each table, the table length is printed first, followed by the dependent 
variable array and the independent variable array in adjacent columns. 



The tables are grouped under the following categories: 

(1) Stepsize control tables. 

(2) Wall tables. 

(3) Geometry and edge tables. 

If the wall contour and pressure tables were entered in the TDKINP 
namelist, they will be printed here as RWTAB, XTABRW, PETAB, and 
XTABPE. 

(4) Experimental profiles. 

This pa je will appear only if experimental profiles for u and h 
were input. 
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SUBROUTINE PARAMS 


Subroutine PARAMS is called at each station of the boundary layer solution 
to evaluate the displacement thickness 6*, momentum thickness Q, and other 
profile parameters of inteiest, based on the converged solution at that station. The 
quantities calculated and the method of calculation are summarized here. 

(1) Mass addition correction to displacement thickness: 


( 2 ) 


(3) 


ADDMAS = f m fr 1 ds 
J w * w 

s init 

Trapezoidal integration is used. 
Displacement thickness: 


= ( 1 - -££-) -fe + ftD . PWS . I 

1 o »e u e ^ 0 e u e r w J 

The integral over y is evaluated using Simpson's rule 
Throat radius corrected for displacement thickness: 


r*' = r*L - 6* cos © 

w 

Once the boundary layer solution passes the throat (specified by the 
value input for XSTAR) , the normalized throat radius r* =T*/L and 
6 at the throat are obtained by interpolation from the input wall 

TV 

contour table. 6* at the throat is obtained by interpolating between 
the values of 6* calculated at the last two stations. 

(4) Momentum thickness: 


e = LC I e 


Pu 

Pe u e 


e 



The integral over y is evaluated using Simpson’s rule. 


i 
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(5) Skin friction calculated by integrating the momentum equation: 

y. 


r I = - 


u u 
*r r 


LC 


J. 


r dp - 

?"[f u tr + G ’ " uE) If + "sr j dy 


where 


(7) 


( 8 ) 


(9) 


E = p Gp . JL 


r R* r * 


The integral over y is evaluated using Simpson's rule, (tj is 
evaluated for diagnostic purposes only and does not appear in the 
summary table.) 


(6) Wall shear stress: 


T w ^r u r Lf ^w ay 


au 


w 


A parabolic wall derivative is used. 
Lot al shear stress coefficient: 


C, — 


G' 


Mw-jjS- 

5v lw 


'f £Re 1_ 

T 2 

A parabolic wall derivative is used. 


2 <>e V 


Heat transfer rate: 
a 


- q 


w 


V; U r G' **w 

Lr Pr 
^ w 


“§■“1 + ( Le-1) S h. 


w 


ac. 

k ~ 


A parabolic wall derivative is used. 
Stanton number: 


St = 


r Ke n u Th - H T 
s, r Pee e w 

t(Le - u pk^r| 




Pr. 


w 


[ 


ah. 

ay 


w 




ra . j'lc wall derivative is used. 
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? 

(10) Integial of heat transfer rate: 

O w --(2,)iL) +1 f V J, ds 

s lnlt 

The integral over s is evaluated using the trapezoidal rule. 


Ji 


& 


* 


2-46 


mm 



i 1 


SUBROUTINE PH0ENX 

Subroutine PHJ0ENX takes an array in which every nth value is known and, 
given an argument array, interpolates for the missing values in the array. The 
calling sequence is as follows: 

CALL PH0ENX (V, Y , N , L) 


where 

V is the dependent variable array, and V(l), 

V(N+1) , V(2*N+1) V(L) are known 

(previously calculated) values in this array. 

Y is the independent variable (argument) array 

such that Y(I) is the argument value correspond- 
ing to V(l) , whether or not V(Q is known. 

N is the interval (integer) at which values in 

the array V are given. 

L is the length of the arrays V and Y. 

Subroutine PH0ENX first packs every Nth value (including the last value) 
from the arrays V and Y into the auxiliary arrays F and X, respectively, thereby 
generating a table of V versus Y. Subroutine XNTERP is then called to interpolate 
for each of the missing values in V, using their corresponding Y-values as arguments. 

Subroutine PH0ENX is used as a means of reducing the execution time for 
hydrogen-oxygen equilibrium cases. The equilibrium properties are solved at every 
nth mesh point rather than at every point, and PH0ENX is used to interpolate for 
values of the equilibrium properties at the remaining points. 



SUBROUTINE PRINT 


Subroutine PRINT is called at the end of every station in the boundary 
layer solution to store the summary table quantities for that station, write a 
summary table record on scratch unit IDRUM if it is time to do so, produce a 
printout of the boundary layer protiles if it is time to do so, and advance the 
print station counters before the boundary layer solution proceeds to the next 
station. 

Subroutine PRINT first stores all summary table quantities for this station 
into the proper row of the array SUMARY. If NSTA rows of SUMARY are now full 
or if this is the last station of the boundary layer solution, the contents of 
SUMARY are written onto logical unit IDRUM to be read and printed later by sub- 
routine SUMTAB. (See Description of Program Output for a detailed description 
of the quantities printed in the summary table.) 

The print station counters are then tested. If ILPRNT=NLPRNT, it is time 
to print the GROUP I profiles, and ILPRNT is accordingly set to zero. If ISPRNT= 
NSPRNT, it is time to print the Group II profiles, and ISPRNT is set to zero. (See 
Description of Program Input for explanations of NLPRNT and NSPRNT. See Description 
of Program Output for details of the Group I and Group II profiles.) If this is the 
initial station, the final station, or an XLIM station, ILPRNT and ISPRNT will have 
been set to zero outside of subroutine PRINT. 

If ILPRNT is now zero, the main boundary layer parameters and the Group I 
profiled are printed at this station. For an equilibrium case ISPRNT is tested and, 
if it is zero, the Group II profiles are printed. 

Whether or not the profiles are printed, subroutine PRINT advances the 
print station counters ILPRNT and ISPRNT by one before returning to the calling 
routine, subroutine EXECUT. 



SUBROUTINE PR0FIL 


This subroutine calculates profiles at the initial station of the dependent 
variables u, H, h, a, and qv across the boundary layer. 


If experimental profiles are input for u, h, and a, then the values of these 

dependent variables a f each mesh point are obtained by interpolating in the input 

tables UPR0F vs YBYNU, HPR0F vs YBYNH, and APR0F vs YBYNA, respectively. 

2 

Then H = h+ u /2at each mesh point. 

If experimental profiles are not input, the velocity profile between 0,1 

and 0.9 y is calculated according to the input power law (See PLAW in the Description 
© 

of Program Input): 


u 



_i_i PLAW 


Between the wall and 0. ly , a linear velocity curve is fared into the power law 

V 

profile; u = u between 0.9y and y . 

2 

For an incompressible case, h = h w at each mesh point, and H = h w + u /2. 
For a compressible case. 


H=H » + T- (H e-V- 

e 


Then h * H - uV2. 


The initial a-provile, if not input, is calculated as follows: 


(x + u ^w^' 

e 

where a w and a e are the input quant' ^es AFWALL and AFEDGE. 

Finally, at each mesh point, 

p v = m w + • 

J e 

Subroutine PR0FIL then sets back values of the dependent variables equal 
to their forward values initially. 


SUBROUTINE RDTAPE 


A restart capability exists in the Mass Addition Boundary Layer Program and 
is described in the section on Additional Program Features. 

Subroutine RDTAPE is called at the beginning of a restart case to read 
previously saved data for the proper station from the restart tape (logical unit 
ITAPE). Since data for a series of stations may be saved on the tape, subroutine 
RDTAPE reads the value of the station counter from the beginning of each logical 
record and compares it to the desired station IRSRD. If the correct station is located, 
the tape is backspaced to the beginning of the record and the restart data is read. 

If data for station IRSRD is not present, the following message is printed: 

THERE IS NO RESTART DATA FOR STATION [iRSRD] . 

The case is then terminated. 



I 


li 
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r-.F. ROUTINE STEP 


Subroutine STEP is called at the beginning of a new station calculation to 


determine the stepsize a s to be used and to obtain the wall contour properties 
r 

r..., -ri“ » and 0 at the next station. 


w ' dx 


w 


dr. 


w 


Subroutine STEP first moves the forward values of r , , and 0 to 

W CLX W 


their back positions, since the current x-pcsition will now become the back station. 
Thus 


THW(l) ^HW^) 

RW(1) = RW(2) 
DKWDX(l) = DRWDX(2) 


The stepsize a s is given by 
DX 


AS = 


cos 0 


w 


where DX is determined from the series of criteria explained below. 

(1) i: 


.-5 


X + 10 "’a XLIM(IDX) 


then at the last station the stepsize was set to hit x = XLIM(IDX) , and the next 
stepsize is specified by 


DX = DXLIM(IDX) . 


The pointer IDX is incremented by one to indicate the next XLIM - value to be 
reached. DX, if obtained here, is not altered by any subsequent criteria except 
possibly (5) . 


(2) If criterion (1) above did not govern the stepsize selection, subroutine STEP 
compares X with XTABSK(ISK) . If X < XTABSK(ISK) , then SKTAB(ISK) is the k-factor 
to be used in determining the next stepsize. If X a XTABSK(ISK) , ther .t J j time to 
change to a new k-factor; the pointer ISK is incremented by one, and X is compared 




with the new XTABSK(ISK). This process is repeated, if necessary, until X <.XTABSK(ISK)l 


Then DX is set equal to the previous stepsize time the k-factor: 

DX = DX*SKTAB(ISK) . 

(3) If this is the first station calculation, 

DX = DXI (input) , 

(4) DX, determined from criterion (2), is checked to see whether it exceeds the 
next XLIM-value to be reached. If 

X + DX + 10" 6 < XLIM(IDX) , 

DX is unchanged; otherwise subroutine STEP sets 
DX = XLIM(IDX) -X, 

so that the solution will reach XLIM(IDX) exactly. The print flags are also set 
so that a printout will occur at the forward station. 

(5) Finally, DX is checked to see whether it exceeds XMAX. If 

X + DX + 10" 6 <XMAX, 

DX is unchanged; otherwise, subroutine STEP sets 
DX = XMAX - X. 

After determining DX from one of the above criteria and setting the stepsize 

As = DX/cos 0 , STEP calls subroutine XNTERP, with X + DX as the argument, to 
w 

dr 

obtain RW(2) and DRWDX(2) , the values at the forward station of r w and » 
respectively. Finally, 

THW(2) - e w - tan' 1 . 


SUBROUTIN E SUMTAB 


Subroutine SUMTAB is ..niied at the end of a case or (by subroiXine DEBUG) 
ueu.-e i fetal error is encounterec to read records of summary table data from 
Icr.ica. ...dt 17 and print tho summary table as part of the program output. 

The summary table data is saved up in core and written on unit 17 several 
stations at a time by subroutine PRINT. The summary table contains the values 
of the main boundary layer parameters of interest at each station. The exact format 
of the summary table output and the items contained in it are described in the 
section on program output. 

The integer flags, counters, and constants in subroutine SUMTAB have 
the following interpretations: 


IDRUM 


Fortran logical unit on which summary table 
data is written. (=17) 

NREC 

• 

Number of logical records of data currently 
on unit IDRUM. 

NSTA 

- ’ 

Number of stations per logical record of 
summary table data. (=13) 

ISTA 

- 

Number of stations of data currently saved 
in core and not yet written on unit IDRUM. 

NVAR 

- 

Number of variables per station. (=29) 


Subroutine SUMTAB first checks if there are any stations of data in core 
(in the array SUMARY) which have not yet been written on unit IDRUM. If ISTA 
is non-zero, which may be the case if a fatal error has been encountered, SUMTAB 
writes a logical record of ISTA stations of data before rewinding unit IDRUM. 

Subroutine SUMTAB then reads and prints each record on unit IDRUM in 
turn. One logical record on unit IDRUM corresponds to one page of summary table 
output. 

The quantities (r w - 6 * cos 8) and x, the tenth and twentieth items of data 
for a given station, are divided by RSTPR, the throat radius corrected for displacement 
thickness, before being printed, so that these data may be used for the corrected 
wail contour table input to the Two-Dimensional Kinetics Program (Reference 1) . If 
the nozzle throat was not reached in this case, RSTPR will be unity. 


SUBROUTINE TABLES 

This subroutine is called from the main program at the beginning of a 
case to convert the input tables to their normalized forms and obtain the wall 
and edge conditions at the initial station. Where quantities are needed at the 
back station, subroutine TABLES sets those quantities equal to the values at 
the initial station. 

Subroutine TABLES replaces values of x/x N in the argument arrays 
XTABRW, XTABSK, XLIM, XTABMD, XTABTW, and XTABPE by corresponding x/L 
values. Values of r^/x^ in the dependent variable array RWTAB are also normal- 
ized. Subroutine TABLES then obtains from the input tables the initial values of 


w' dx 


* 0 W , Ax, As, m w , and T w . 


Pressure and edge velocity tables are handled differently for a hydrogen- 
oxygen case and for a perfect gas case. These are described in turn. 


Hvdroaen-( 


ten Equilibrium Case 


For an equilibrium case, a pressure table PETAB vs XTABPE has been input, 
either in the DATA namelist or in the TDKINP namelist as PITAB vs. XTTAB. (Since 1 
PITAB and XHAB are equivalenced to PETAB and XTABPE in the main program, the 
pressure table can always be referenced as PETA" vs. XTABPE.) Subroutine TABLES ■ 
calls subroutine H00DE to do an isentropic expansion at the edge of the boundary 
layer to obtain a normalized edge velocity table UETAB vs. XTABPE. (See subroutine 
H00DE for details.) Subroutine H00DE is also called to perform a temperature- 
pressure equilibrium calculation at the wall, givir g initial values of h w and H w . 

Perfect Gas Case 


For a perfect gas case, subroutine TABLES calls subroutine IG0DE with the 
initial value of T w and the input value PEDGE to obtain h w . Then h w = h w /u a 
and H w * h w . The enthalpy at the edge of the poundary layer is obtained similarly. 
Subroutine IG0DE is called with the input initial values TEDGE and PEDGE to obtain 
the initial values of h and p . Then h = h /u* and 

A “a A A *• 


? 


''’iere u g is the normalized initial edge velocity, obtained from the input value 


ui.DGE: 


u 

u = UEDGE • 
e u 


N 


For a perfect gas case, either an edge pressure table PETAB vs. XTABPE 
or an edge velocity table UETAB vs. XTABPE may be input. If a pressure table 
was input, a velocity table is generated point for point from it as follows (The 
subscript 0 here denotes the initial station.): 


h e = h eO ( f^ ) 
eO 


r-i 

y 

1/2 


u e -<2<H e - h e» 

If a velocity table was input, a pressure table is generated from it as follows: 


P e = P eO + T PeKo ' u ? u r' 


where o e > the edge density, was obtained from subroutine IGODE above. 

Velocity Derivative Table . . 

At this point, for an equilibrium or perfect gas case, both a velocity 
fable and a pressure table are available. The initial value of the edge velocity 
UEDGE is now obtained by interpolation in the table UETAB vs. XTABPE (whether 
input or calculated as described above). The edge velocity table is printed out; 
then a velocity derivative table DUDXT vs. XTDUDX is generated in its place. A 
linear derivative is taken between two velocity table points i and i + 1: 


du 


u 


- u 


_e_ = 


dx 


'i+1 


'i 


x i+l " x i 


du 


e 


This is then stored in the velocity derivative table as the value of -g— at 
1 


x = 


+ X i+ l)* 


i 


J 
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Each set of points in the velocity table UETAB vs. XTABPE is treated in this 
manner. 

du 

Finally, the initial value of is Stained by linear interpolation in 
the new table DUDXT vs. XTDUDX, and 


du„ _ du 
e — e 

ds dx 


cos fl w . 


The initial value of edge pressure PEDGEB is then obtained by interpolation, 
and control returns to the main program. 
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MAIN PROGRAM (TFCBL) 

The Mass Addition Boundary Layer Progre m main driver is responsible for 
reading and processing the program input and obtaining all boundary layer quantities 
at the initial station, X = XINIT, At this point, subroutine EXECUT, the main 
execution control subroutine, is called to integrate along the boundary layer from 
X = XINIT to X = XMAX. 

The main program code also contains a master list of all common blocks 
used throughout the boundary layer subroutines. 

Specifically, the main program does the following (in sequence): 

(1) Sets program constants. 

(2) Sets nominal values of input data. 

(3) Reads the DATA namelist input. 

(4) Reads the TDKINP namelist input if it is present. 

(5) Calls subroutine NL0UT to print the program input. 

(6) Calls subroutine H00DE, for an equilibrium case, to read 
the THERM0 and REACTANTS data and initialize the equili- 
brium portion of the program. 

(7) Sets constants whose values are based on input data. 

(8) Initializes the independent variables x and s and the 
boundary layer normalization function f . 

(9) Calls subroutine TABLES to normalize the input tables, and 
evaluate the boundary layer wall and edge conditions at the 
initial station. 

(10) Initializes the derivative £' of the boundary layer normal- 
ization function. 

(11) Calls subroutine GFUNC to evaluate the boundary stretching 
function G and set up arrays of y , y, G' , and G" at each 
mesh point. 

(12) Calls subroutine PR0FIL to obtain initial values of the main 
dependent variables u, h, H, a, and pv at each mesh point. 

(13) Calls subroutine IG0DE for a perfect gas case or subroutine 
H00DE for an equilibrium case to initialize the thermodynamic 
and laminar transport properties at each mesh point. 

(14) Presets turbulent quantities for a laminar case or calls subroutine 
EDEY to calculate the turbulent transport properties at each mesh 
point for a turbulent case. 
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(15) Moves forward values to back values. (Both are equal at 
the initial station.) 

(16) Calls subroutine PARAMS to calculate the gross boundary 
layer parameters at the initial station. 

(17l Calls subroutine PRINT to print the results of the above 
initial station calculations. 

(18) Calls subroutine EXECUT to solve for ->11 boundary layer 
quantities from X = XINIT to X = XMAX. 

(19) Following completion of the boundary layer solution, prints 
the throat radius corrected for displacement thickness. 

(20) Calculates and prints the thrust loss, if the ambient pressure 
was input. 


SUBROUTINE TRIM 


This subroutine solves the system 
AX = b 

for X where A is tridiagonal, i.e.; 

all a^ = 0 if i > J+l or i < j-i. 

The method used is Gaussian elimination. 

The calling sequence is as follows: 

CALL TRIM (A,X, B, N, NN) 

where 

A Is the input coefficient matrix and must be dimensioned at least 
A(N, 3) . The elements a^ must be input as 


— 

A( 1 , 2) 

A( 1 , 3) 


— 

a ll 

a 12 

A(2, 1) 

A(2,2) 

A(2,3) 


a 21 

a 22 

a 23 

A(3, 1) 

A(°, 2) 
• 

• 

A(3,3) 


a 32 

a 33 

• 

• 

a 34 

A(N, 1) 

• 

• 

A(N, 2) 

— 


a n,n-l 

• 

• 

a 

nn 

— 


The contents of A are preserved by TRIM. 

X Is the output solution vector X and must be dimensioned at least 
X(N). 

B Is the input vector b and must be dimensioned at least B(N) . The 
contents of B are preserved by TRIM. 
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N Is the order of the system. N & 2 is required. 

NN Is the dimension of A (NN,3), X(NN) f and B(NN) in the calling 

program . 

Subroutine TRIM is called (1) by subroutine M0MNTM to solve the 
momentum equation for velocity at each mesh point in the boundary layer, (2) 
by subroutine ENERGY to solve the energy equation for total enthalpy H at each 
point, and (3) by subroutine ELEMTS to solve the element conservation equation 
for a, the element mass fraction of hydrogen, at each - jint. 
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SU BROUTINE XNTERP 

iven a table of n points of independent variables Xj and dependent 
vanabl y^, XNTERP returns a value of y and its first derivative y' for a given 
value ■- 1 x and saves the place in the table. 

Subroutine XNTERP is used widely throughout the program, specifically 
in the following contexts: 


(1) To interpolate in tables of r , versus x, P versus x, and u 

w e 6 

versus x* 

(2) To generate initial profiles across the boundary layer for u and h 
from input experimental data. 

(3) To interpolate in arrays of y and y. 

(4) To obtain thermodynamic and laminar transport properties at mesh 
points intermediate to those for which the equilibrium chemistry 
was solved explicitly. (See subroutine PH0ENX.) 

The calling sequence is as follows: 

CALL XNTERP OCIN, Y0UT, YP0UT, KIN, XAR, YAR, 

IAR, CAR, IP0S) 

where 

XIN is the input value of the argument x. 

Y0UT is the value of the dependent variable 

y returned. 

YP0UT is the value of the first derivative y' 

returned. 

KIN is a pointer which is set by XNTERP upon 

return to the calling program, such that XAR (KIN) 
is the largest entry in the independent 
variable list XAR which is less than or equal 
to the argument XIN. This pointer must be 
initialized at zero prior to the first call to 
XNTERP for a given table. 

XAR is the array of independent variables x^ stored 

in increasing order. 
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YAR is the array of dependent variables y. 

corresponding to the 

IAR is the number of values in XAR and YAR. 

CAR is a vector of six cells used by XNTERP to 

store the coefficients of the parabolas used 
for interpolation. (No user initialization is 
necessary.) 

IP0S is a pointer, set by the user as follows: 

= 0 on the first entry to XNTERP for a 
given table; 

= the subscript, on subsequent entries, 
of the independent variable array XAR 
at which the argument search is to begin. 


Subroutine XNTERP searches the independent variable table until the 
points bracketing x are found (x^ < x < Xj+j)* Using the four surrounding points 
x i-l ' x i' x i+l' anc * x i+2' evaluates the coefficients of parabolas y L , 

passing through the three leftmost points, and y R , passing through the three 
rightmost points. 


y L = C 1 + C 2 + C 3 

y R = C 4 + C 5 (x-Xj) + Cg (x-Xj ) 2 

A cubic weighting function a is then defined as follows: 



IXx) = 


x - X 


i 


x i+l - x i 


a(x) = 3U 2 -2U 3 
a’(x) = (6TT-6U 2 ) U’(x) 


Finally y and y 1 are calculated: 

y - (l-od y L + ay R 

y* = (l-o) Yl - a’ y L + a* y R + ay R 

If the argument x is between the first two or +he last two points in the 
table, only one parabola is used. In a two point table, linear interpolation is 
used. For a single point, y is constant for all x, and only YAR(l) need be specified. 


i 
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Subroutine XNTERP will not extrapolate in a table of four or more points. 

Rather, the following error message is printed: 

XNTERP 0UT P? RANGE ... X - X(l) = X(N) = 

Y(l) = Y(N) = 

XNTERP then calls subroutine SUMTAB to print the summaiy table thus far and 
finally calls EXIT, terminating execution. 

Subroutine XNTERP is programmed to take maximum advantage of the searches 
and calculations performed on the previous entry to it. Through use of the pointer 
IXIN, XNTERP, once it locates the argument in the table, ascertains whether either 
or both sets of parabola coefficients from the last entry can be reused and, if so, 
does not recalcul ? *e those coefficients. The user can further enhance the subroutine's 
efficiency if he has some foreknowledge of where the argument lies in the indepen- 
dent variable table. For example, if it is known that the argument is monotonically 
increasing, then the begin-search pointer IP0S can be set equal to IXIN following 
return to the calling program so that for the next interpolation the argument search 
will begin at XAR (IXIN) , where it left off on the previous entry. An even simpler 
way to do this is to make the actual arguments for IXIN and IP0S the same variable. 

This subroutine i^ the same as that used to interpolate in the contour 
and pressu.a tables in the ICRPG Turbulent Boundary Layer Program (Reference 2) . 
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SUBROUTINE ZFUNC 

Subroutine ZFUNC iteratively calculates the value of the boundary layer 
normalization function £ and its derivative £' . Although the non-iterative portions 
of the C -calculations take place in subroutine EXECUT and not in subroutine ZFUNC, 
the entire method is explained here for clarity. 

The quantities involved in the C. -calculation and their corresponding program 
symbols are given in the following table: 


Ouantitv 

Symbol 

Cti 

ZSTAR(l) 

AS_j 

DSZ(l) 

Cold 

ZETA0 

Cave 

ZETA 

r 

ZETAN 

^new 


c* 

ZETAP 

e* 

ZSTAR(3) 


At the beginning of a new station calculation the quantities C old » V • As _i' 

and a s are known. (Refer to the accompanying diagram.) 

C = slope of line $> ld C ave ^ 


• r* 

T 1 - 1 i 


* f * 


i h~As_ 1 -|- as H 

III! 

old new 

Subroutine EXECUT obtains a value of £ new for iteration 0 by integrating 

Cnew “ Cold + C AS ‘ 


Then 


C wo ^ C old + C new ^ ^ ' 
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i’ • boundary layer equations are solved, using this value of r , to obtain 

new 

u, u h and oat each mesh point. 

subroutine EXECUT then cails subroutine ZFUNC to update f ' and if 

s v new 

full ..*itine ZFUNC proceeds as follows: 

(1) Y r , the value of y where u = 0.99 u , is obtained by interpolation. 

C e 

(2) y^ , the corresponding value of y, is found by interpolating in 
the table of y vs y . 

(3) A corrected value of £is calculated: 

C* = C new y c 

(4) r 1 and C ... are recalculated as follows: 

* new 

{■-({*- CjlAiS.j + as); 

= Cold + f * s - 

Also, 

^ave ~ ^old + C new ^ ^ * 

(5) m is recalculated, based on m normalized by the updated 

w m+l/2 w m+l 

value of If . (pv = m is the wall boundary condition for the 

- new w m +l/2 

continuity equation, which is solved following execution of subroutine 
ZFUNC.) 

If additional iterations are performed at a given station, subroutine ZFUNC is called 

every iteration to recalculate C ' , r . and r . „ based on the latest iterated 

* - new ’ ave 

velocities . 

At the end of a station, subroutine EXECUT updates the remaining £- related 
quantities to prepare lor the next station calculation. Specifically, 


c 

( Cne«-C*l )/(i5 -l + i s) 

e*i = 

Cold 

Cold - 

Cnew 

it 

i 

<0 

<3 

AS 


The boundary layer solution then proceeds to the next station 



3. DICTIONARY OF PROGRAM VARIABLES 


AH program variables of principal importance are defined in this section, 
with the exception of those quantities which are input. The variables are grouped 
according to function. 


SUBSCRIPTS 

Symbol 

NY 

J0 

JN 

JA 

NEL 

NSP 

I 

J 

IEL 


Meaning 

Number of points in boundary 
layer U 250) . 

Back station m. 

Forward station m+1 . 

Average of back and forward 
stations (m+l/2). 

Number of chemical elements (s 2) . 

Number of chemical species, (s6) . 

Mesh point subscript, ranging from 
1 (wall) to NY (edge). 

Station subscript, taking on the 
values J0, JN, and JA. 

Element subscript for equilibrium 
case: 

=1 for hydrogen (fuel) ; 

=2 (NEL) for oxygen (oxidant) . 

Species subscript for equilibrium 
case: 

=1 for H; 

=2 for H 2 ; 

=3 for H 2 0; 

=4 for O; 

=5 for OH; 

=6 (NEL) for Og. 


ISP 


LOGICAL UNITS 


Svmbol 

Meaning, 

IDRUM 

Logical unit on which summary table 
data is temporarily stored (=17) . 

ITAPE 

Logical unit on which restart data 
is saved (=16) . 

COUNTERS 

Svmbol 

Meaning,. 

ISTATN 

Station number. 

ITER 

Iteration number (reset to zero at 
beginning of each station) . 

ILPRNT 

Number of stations since Group I 
profiles were last printed. 

ISPRNT 

Number of stations since Group III 
profiles were last printed. 

LNSPPG 

Maximum number of lines per page 
of output. 

LINE SR 

Number of lines remaining for output 
on current page. 


ANALYTIC QUANTITIES 

The analytic variables are described below in alphabetical order. Dimensioned 
variables appear with subscripts. Variable subscripts take on the meanings given 
above under SUBSCRIPTS. Where a constant appears as a subscript, it represents 
the maximum dimension of the variable. 



J 


Symbol 


Meaning 


A(I, 3) 


AFEDGE 

AFWALL 

ALPHA(I, J, IEL) 

AV(D 

B(D 

BCF 

BGP(D 

BGPP(l) 

BLE(I, J) 
BLET(I,J) 
DLSTAR 
DLSTTH 

DPEDSN 

DRWDX(2) 


Array containing the difference equation 
coefficients. Each column of A represents 
a diagonal of the tridiagonal matrix. 

a , mass fraction of hydrogen at the 

w 

edge of the boundary layer. 
a , mass fraction of hydrogen at 

W 

the wall. 

Element mass fraction. 

Sound speed. 

Right-hand side of difference equations. 
Cf, local shear stress coefficient. 

A/ 

G' (y) , first derivative of the boundary 
layer stretching function. 


G" (y) , second derivative of the 
boundary layer stretching function. 


Le, Lewis number. 

Le T , turbulent Lewis number. 

6*1 displacement thickness. 

6*^, displacement thickness at the 
throat. 


dP 


e_ 

_ 

ds 

m+1 ^ e r 

dr 

dr , 

w 

m and “cbT 

dx 


u 


e 


m+1 


ds 


m+1 


m+1 


DS 

DUEDS 


AS, wetted length stepsize. 



m+1/2 
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DUEDSN 


m+1 


DUEDS0 


ds 



m 


DX 

DY 

E(D 

EPS 

F(D 

H(I,D 

HEDGE 

HWALL 


Ax, axial stepsize. 

AY, spacing between mesh points. 

E, auxiliary quantity in difference 
equations . 

e, eddy viscosity. 

F , auxiliary quantity in difference 
equations. 

H, total enthalpy. 

H e , total enthalpy at the edge of 
the boundary layer. 

H w » total enthalpy at the wall. 


0DYSQ 

02 DY 

04 DY 
PEDGEB 
PR(I,J) 
PRT(I.J) 

REYINF 

RH0(I,J) 

RH0V(U 

RSTPR 


1 

(ayP 

1 

2 ay 
1 

4 ay 

P , unnormalized edge pressure, 
e 

Pr , Prandtl number. 

Pr^, turbulent Prandtl 
number. 

Re r , reference Reynolds number, 
p, density. 

( P v) m+l/2 

r*' , throat radius corrected for 
displacement thickness. 
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RW(2) 


S 

SCI(I, J, ISR 

SH(I.J) 

SHEDGE 


SHI(I, J, ISP) 


SHWALL 

SIG1 (3) ^ 
SIG2 (3) 
SIG3 (3) 
SIG4 (3) 
SIG5 (3) , 

SMDW 


SMDWN 


SMDWJZ) 


SMU(I, J) 
SNTGRL 


SOW 

SQWDS 


r and r , wall radius at the 
w m w m+l 

back and forward stations, 
s, wetted length. 

, species mass fraction. 

h, enthalpy. 

h , enthalpy at the edge of the 

W 

boundary layer, 
h^, species enthalpy. 

h w> enthalpy at the wall. 

Temporary storage fora j# a£« Oy 
o^, and a 5 # auxiliary quantities in 
the difference equations. 

rn .mass addition rate. 

w m+l/2 

m .mass addition rate. 
w m+l 

m w , mass addition rate, 
m 

4 , viscosity. 

m w £r^ ds, part of the mass 

s init 

addition correction to the displacement 
thickness. 

-q , heat transfer rate. 

IV 

W jl ("J 1 %, « 

s init 


STAN 


Stanton number. 


Taj) 

TAUW 

TEDGE 

THETA 

THW(2) 

TWALL 

uaj) 

UEDGE 

X 

Y(D 

YTIL(D 

ZETA 

ZETAN 

ZETA0 


T, temperature. 

T , wall shear stress. 

T , edge temperature. 

6 

0 , momentum thickness. 


0 and 0 , wall contour 

w m m+1 

angle at the back and forward stations. 


T , wall temperature. 

u, velocity in the s-direction. 
u e# velocity at the edge of the 
boundary layer. 

x, axial distance. 


y, non-dimensional normal coordinate. 


y * , normal coordinate. 


^m+1/2' boundary layer normalization 
function. 


^m+1' boundary layer 
normalization function. 


f . boundary layer normalization 
function. 1 


ZETAP 


derivative of the boundary layer 
Layer normalization function. 


4. DESCRIPTION OF PROGRAM INPUT 


Program input for the Mass Addition Boundary Layer Program is in four 
parts: (1) the case title; (2) the DATA namelist; (3) optionally, the TDKINP 
namelist for wall contour and pressure tables which have been punched from the 
Two-Dimensional Kinetics Program (Reference 1); (4) for a hydrogen-oxygen equili- 
brium case, the THERM0 data and REACTANTS data for the one-dimensional equilibrium 
chemistry calculations. These are described in turn. 

CASE TITLE 

The first card contains the case title in columns 1-60, to appear at 
the top of each page of output. 


DATA NAMELIST 


This group of input follows the standard NAMELIST format. The first card 
of this group must contain $DATA starting in column 2. The last card must contain 
SEND starting in column 2. 

All program input may be entered in this group; however, if pressure and 
wall contour tables have been punched from TDK, these must be entered in the 
TDKINP namelist described iater. 

The input variables and arrays in the DATA namelist are described in 
functional groups below. For some of the variables, nominal values are preset 
in the program and will be used if not overridden by an input value. Where 
applicable the nominal values aia indicated. 


Flags and Options 
Name 
IDEAL 

LAM NR 
INC0MP 
J2D 
INTDK 


Description 

Integer flag, set as follows: 

=1 (nominal) for a perfect gas case; 

=0 for a hydrogen-oxygen system 
with equilibrium chemistry. 

Integer flag, set as follows: 

=1 for laminar flow; 

=0 (nominal) for turbulent flow. 

Integer flag, set as follows: 

=1 for incompressible flow; 

=0 (nominal) for compressible flow. 

Integer flag, set as follows: 

=1 (nominal) for axisymmetric geometry; 

=0 for a two-dimensional case. 

Integer flag, set as follows: 

=1 if pressure and wall contour tables 
are input from TDK, in which case 
the TDKINP namelist will follow the 
DATA namelist; 

=0 (nominal) if pressure and wall 

contour tables are input within the 
DATA namelist. (No TDKINP namelist 
will be present.) 


J 


Reference Quantities 


Name 

Description 

BLREF 

L, reference length in it. 
(Nominal value =1.) 

UREF 

u , reference velocity in ft/sec. 
r (Nominal value = 1.) 

RH0REF 

2 y 4 

p , reference density in lb,sec /ft . 
r (Nominal value =1.) 

SMUREF 

2 

a , reference viscosity in lb-sec/ft . 
(Nominal value = 1.) 

Input Normalization Constants 



If data to be input to the Mass Addition Boundary Layer Program exists 
in a nondimensional form (e.g. a nozzle contour table may have been, normalized 
by throat radius), it need not be dimensionalized prior to input. The following 
input normalization constants are used by the program to multiply the input values 
to which they correspond immediately after the data is read. 


Name 


Description 


XN (ft.) 


YN (ft.) 


UEN (ft/sec.) 
PEN (lb f /ft 2 .) 


Xj^j, normalization factor by which x- 

aistances and lengths have been normal- 
ized upon input. (Nominal value =1.) The 
following input variables and arrays are 
assumed by the program to be normalized 
by XN: 

SINIT, XINIT, XMAX 

Stepsize control tables XLIM, XTABSK 

Argument tables XTABTW, XTABMD, 

XTABRW, XTABPE 

y^j, normalization factor by which the 

experimental profile argument tables YBTtNU, 
YBYNH, and YBYNA have been normalized 
upon input. (Nominal value *1.) 

Ujg, normalization factor by which the velocity 

UEDGE and edge velocity table UETAB have 
been normalized upon input. (Nominal value -1.) 

Pjj, normalization factor by which the pressure 

PEDGE and edge pressure table PETAB have 
been normalized upon input. (Nominal value *1.) 
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\ 

\ 

Name 

SMDN (lb f sec/ft?) 

i 

I 

} 

i 

Problem Limits and Initial Values 
Name 
SINIT 

XINIT 

I 

XMAX 

DXI 

DELTA I (ft.) 

ZETAPI 

Edge Quantities 
Name 
UEDGE 

PEDGE 

TEDGE (°R) 

AFEDGE 



Description 

m^j, normalization factor by which the 

mass addition table SMDTAB has been 
normalized upon input. (Nominal 
value =1.) 


Description 

^init^ x N' * n *** a ^ value of wetted length 

s' at which boundary layer solution begins. 
(Nominal value = 0.) 

initial value of axial distance 

3f at which boundary layer solution begins. 
(Nominal value = 0 .) 

x max^ x N* vain- of axial distance If 
at which boundary layer solution terminates. 

A x init 

Ax. = — r , initial stepsize (normalized) . 

init L 

6, initial boundary layer thickness, in 
units of y . 

Initial value of £' , the derivative of the 
boundary layer normalization function £. 

If no value is input, the program will 
calculate one. 


Description 

u e /uj y, edge velocity at initial station. 

P e /Pjr , edge pressure at initial station. 

T , edge temperature at initial station. 

6 

a e> constant mass fraction of hydrogen at 

the edge of the boundary layer for a hydrogen- 

oxygen case. (Do not enter for a perfect 

gas case). Mixture ratio, MR = (1-a )/<* 

e e 
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Constants 


Name 

AFTRNS 

PRI 


GAMMA 

FM0LWT 

PIAW 


PAMB 


GP0* 


Description 

Otrans* mass fr action hydrogen in 

the transpiring gas . (Do not input for 
a perfect gas case.) 

Optional input Prandtl number Pr. If 
input, this constant value will be used 
for either a perfect gas or hydrogen- 
oxygen case. If not input for a perfect 
gas case, Pr = .72 throughout. For the 
hydrogen-oxygen system, a variable Pr 
calculated by the One-Dimensional Equi- 
librium (ODE) subroutines will be used if 
PRI is not input. 

v , specific heat ratio for a perfect gas case. 
(Do not input for an I^-Oj equilibrium case) . 

M, molecular weight for a perfect gas case. 
(Do not input for an Hg-Og equilibrium case). 

Exponent parameter p in velocity power law 
initial profile: 

u » Vp 

u •< oi y > 

If PIAW is not input and experimental initial 
profiles are not used, a power law velocity 
profile with p = 1 will be generated. If ex- 
perimental profiles are input, do not enter a 
value for PIAW. 

P amb^ P N' am hi®nt pressure for final thrust 

loss calculation. If not input, a thrust loss 
calculation will not be performed. 

G' (0) , the value desired for the derivative 

G* (y) = of the stretching function at 

7 = 0. (Do not input if no stretching is 
desired ) . 


*The use of a stretching function for the normal coordinate is required for turbulent 
boundary layers , but is optional for laminar flows . 


Name 

SN3 


XSTAR 

AFWALL 

UEK (ft/sec) 

RH0EK (lb { sec 2 /ft 4 ) 


4 


Description 

Exponent paramc^r r n for the boundary 
layer stretching function 

y = G(y) ~ b { ^loq( 1 . 718( + a) + l)j 1/n - /?! 

If not input, no stretching will be used, i.e. 
y " G(y) = y; G'fy) = 1; G"(y) = 0, 

x/xjj at the throat. 

<v , initial mass fraction of hydrogen at the 

W 

wall. (Do not input for a perfect gas case.) 

For a hydrogen-oxygen equilibrium case, if 
AFWALL is not input, the program will set 
AFWALL=AFTRNS. 

^ek # non-e< l u i^ r i um e£ Mje velocity at the 

final station of the boundary layer solution (from 
TDK) , for use in the thrust loss calculation. If 
UEK is not input, the equilibrium edge velocity 
will be used in the thrust loss calculation. 

p g k# non-equilibrium edge density at the 

final station of the boundary layer solution (from 
TDK) , for use in the thrust loss calculation. If 
RH0EK is not input, the equilibrium edge 
density will be used in the thrust loss calcu- 
lation. 
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Convergence and Edge Criteria 
Name 
C0NVRG 


EPSLN1 


EPSLN2 


EPSLN3 


Description 


Convergence criterion for iterations on the 
boundary layer solution at a given station. 

If the percentage change in ^ u between 

°y w 

successive iterations is less than or equal 
to C0NVRG, the solution will proceed to the 
next station. (Nominal value = .005). 


If, after a given solution of the momentum 
equation. 



du 

dy 


e 


> EPSLN1, 


the boundary layer thickness is increased 
by adding an additional mesh point at the 
edge. 


If, after a given solution of the energy 
equation. 



>EPSLN2 


the boundary layer thickness is increased 
by adding an additional mesh point at the 
edge. (EPSLN2 need not be input for an 
incompressible case.) 


If, after a given solution of the elemen- 
equation. 


1 

da 


a e 

dy 

e 


>EPSLN3, 


the boundary layer thickness is increased by 
adding an additional mesh point at the edge. 
(EPSLN3 need not be input for an incompressible 
or perfect gas case.) 
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Counters 


Name 

MAXIT 


NYI 


NLPRNT 


NSPRNT 


IYPR 


IYEQ 


Description 


Maximun number of iterations on the 
boundary layer solution at a given station. 

If the solution does not converge within 
MAXIT iterations, the program will proceed 
to the next station using the latest iterated 
valuer horn this station. If MAXIT=0, only 
one calculation will be performed at each 
station. (Nominal value =1.) 

Initial number of mesh intervals across the 
boundary layer (integer) . Thus the initial 
number of mesh points will be NYI + 1. 

(Integer) Boundary layei profiles of u/u e , 

h/h , O/P < °v. f , and T versus y and y 
e e 

are output every NLPRNT stations. These 
profiles are also output at the initial and 
final stations and at those x-values given 
in the array XLIM. After each XLIM -value is 
reached, the station print counter is reset 
to zero, and these profiles are not output 
again until NLPRNT stations have been computed 
(unless the next XLIM -value is reached first). 
(Nominal value = 50) . 

(Integer) Boundary layer profiles of O /? , the 
species mass fractions, ^ , and Pr versus y 
are output every NSPRNT stations for a hydrogen- 
oxygen equilibrium case. (NSPRNT should 
not be input for a perfect gas case.) For an 
equilibrium case NSPRNT, if input, should be 
an integer multiple of NLPRNT. I f NSPRNT is 
not input, the profiles under the u-atrol of 
NSPRNT will be output at the initial, final, 
and XLIM stations only. 

(Integer) When profiles are output, the values 
at every IYPRth mesh point will be printed. In 
addition, properties at the wall and edge are 
always printed. (Nominal value =1; i.e. every 
mesh point is printed.) 

(Integer) For a hydrogen-oxygen equilibrium 
case, equilibrium calculations are performed 
at the wall and every IYEQ mesh points there- 
after. (Nominal value = 4.) Interpolation is 
then used to obtain the equilibrium properties 
:v thr. remain inn mesh points. 
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Tables 




The remainder of the input in the DATA namelist consists of tables. In 
general, there will be three variable names associated with each table: 

(1) an integer variable specifying the length of the table: 

(2) an array of independent variable (argument) values in 
increasing order: 

(3) an array of dependent variable values associated with the 
corresponding values in the argument array. 

The input tables are described in functional groups below. Linear 
inter polation/extrapolation is used to obtain values of the dependent variables 
in the tables TWTAB vs XTABTW and SMDTAB vs XTABMD. An averaged parabolic 
interpolation method is used in the tables RWTAB vs XTABRW and PETAB vs XTABPE. 
No extrapolation is performed in these tables, so that the argument arrays must 
cover the entire range of possible arguments for a given case. 

Stepsize Control Tables 


Name 

Description 

LDXLIM 

Number of entries in the arrays DXLIM and 
XLIM U 50) . 

DXLIM j 
XLIM \ 

DXLIM (I) contains the normalized stepsize 
A x 

Ax = — j- — to be used when the boundary 


layer solution reaches x/xjj = XLIM (I) . The 

stepsize at the previous station is reduced if 
necessary to hit XLIM (I) exactly. A printout 
of the variable profiles across the boundary 
layer occurs at every point in the array XLIM . 

The last entry in XLIM (that is, XLIM (LDXLIM) 
must be greater than or equal to XMAX; otherwise 
the case will be terminated at 


XLIM (LDXLIM). The first entry, 

XLIM (1) , must be greater than 

XINIT. (Note that DXI already specifies the 

initial stepsize.) 

LSKTAB 

Number of entries in the arrays SKTAB and 
XTABSK U 50). 
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Description 


4 


Name 

SKTAB | 
XTABSK ( 


Wall Tables 


Name 

LTWTAB 


TWTAB ( 
XTABTW ( 


LMDTAB 


smdtabJ 

XTABMD f 


SKTAB (1) contains *he stepsize multiplication 
factor k. to be used until x/xjj equals or exceeds 

XTABSK (I) . The are used as follows: each 

succeeding stepsize is calculated as times 


its preceding stepsize. Thus the stepsize Ax 
can be steadily increased, decreased, or held 
constant over specified ranges using the-se 
arrays. Since the stepsize can be reset to ar> 
arbitrary value at any point along the boundary 
layer using DXLTM and XLIM, these stepsize 
control tables aive the user considerable 
flexibility in controlling the boundary layer 
solution. The first entry in XTABSK must be 
greater than XINIT. The last entry, XTABSK 
(LSKTAB) , must be greater than or equal to 
XMAX; otherwise the case will be terminated 
at XTABSK (LSKTAB) . 


Description 

Number of entries in the arrays TWTAB and 
XTABTW (<: 100) . 

Table of wall temperature T w vs x/x^. TWTAB (1) 

is the value of T^ corresponding to the value of 

x/xjj in XTABTW (I). Values in XTABTW must be 

increasing order; however, at a discontinuity 
in T , the argument should be repeated once. 

Number of entries in the arrays SMDTAB and 
XTABMD (<; 100) . 

Table of mass addition m w /mjj vs 5 c/xjj. 

SMDTAB (1) is the value of rr^/m^ corresponding 
to the value of x /x^ in XTABMD (1). Values 
in XTABMD must be in increasing order; however, 
at a discontinuity in tn w /mjj# the argument 
should be repeated once. 
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Geometry and Edge Tables 
Name 
LRWTAB 


RWTAB ( 
XTABRW j 


LPETAB 


PETAB j 
XTABPE | 




LUETAB 


Description 

Number of entries in the arrays RWTAB and 
XTABRW (^500), if these are used to enter 
the wall contour table; OR. if the wall 
contour table is input in the TDKINP namelist, 
then LRWTAB must be the number of entries 
in the arrays YITAB and XITAB. 

Table of wall radius r^/x^j vs x/xjj. RWTAB (0 

is the value of F/x., corresponding to the 

value of x/x N in XTABRW (D . Values in XTABRW 

must be in increasing order. If the wall 
contour table is input from TDK, these arrays 
should not be input. 

Number of entries in the arrays PETAB and 
XTABPE (s 500) if these are used to enter the 
pressure table; OR. if the pressure table is input 
in the TDKINP namelist, then LPETAB must be 
the number of entries in the arrays PITAB and 
XITAB. 

Table of P g /P N vs x/x N . PETAB (I) is the 

value of Pg/Pjj corresponding to the value of 

x/xjj in XTABPE (I) . Values in XTABPE must be 

in increasing order. If the pressure table is 
input from TDK, these arrays should not be 
input. If a pressure table is input, either in 
these arrays or in the TDKINP namelist, then 
the edge velocity table UETAB will be generated 
by the program and must not be input. 

Number of entries in the array UETAB (s 500) , 
if an edge velocity table is input. An edge 
velocity table may be input only for a perfect 
gas case. 





t mKasmtsrsmm 
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Name 


Description 


UETAB j 
XTABPE j 


I 


Table of u e /uj^ vs x/x^. UETAB (I) is the 

value of corresponding to the value 

of x/x N in XTABPE (U . Values in XTABPE 

must be in increasing order. If a velocity 
table is input, then the pressure table PETAB 
will be generated by the program, and the 
pressure table must not be input, either in 
PETAB or in the TDKINP namelist. Note that 
XTABPE is the argument array for both PETAB 
and UETAB but that one and only one of the 
latter is input in a given case. 


The edge velocity table, whether input or 
generated from the pressure table PETAB, is 
used during initialization to generate a table 

du e 

of — vs x. Linear interpolation is 

* performed in this latter table, and the edge 

velocity u is obtained at each station by 

du e 

integrating — , rather than by interpolating 
directly in the table UETAB vs XTABPE. 

4 Experimental Profiles 

Experimental initial profiles for u, h, and a may be optionally input and 
will override the initial profiles normally calculated by the program. If experimental 
profiles are to be used at all, the velocity profile UPR0F vs YBYNU must be input, 
f For a compressible case in which an experimental velocity profile is input, the 
j enthalpy profile HPR0F vs YBYNH must also be input. For an incompressible case, 
j the enthalpy profile should not be input. For a perfect gas case, the a~ profile 

APR0F vs YBYNA should not be input; for an equilibrium case, an experimental a- profile 
may be input but is not required. 


Averaged parabolic interpolation is used in all three experimental profile 

tables. 


Name 

i LUPR0F 

1 

i 


Description 

Number of entries in the arrays UPR0F and 
YBYNU (£50) . 


*t- M-i-s T&inuMJV! : 
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warns. 

UPR0F ( 
YBYNU \ 


LHPR0F 


HPR0F j 
YBYNH ( 


LAPR0F 


APR0F 

YBYNA 


Dg^?r.lPtlgn 

Table of Initial u/CT e vs y/y^j. UPROF (I) is 

the value of u/u e corresponding to the value 
of y/y N in YBYNU (I) . Values in YBYNU must 
be in increasing order. 

Number of entries in the arrays HPR0F and 
YBYNH (*50). 

Table of initial h/h^ * r s y/Vjj* HPR0F (I) is 
the value of h/^ corresponding to the value 
of y/y N in YBYNH (0 . Values in YBYNH must 
be in increasing order. 

Number of entries in the arrays APR0F and 
YBYNA (£ 50) . 

Table of initial a vs y/y^. APR0F (0 is the 

value of a corresponding to the value of 
y/y N in YBYNA (Q . Value in YBYNA must be 

in increasing order. 


Replacement of a Table bv a Constant Value 

To replace a table by a constant value, enter the constant value in the 
first cell of the dependent variable array and make no entry for the length of that 
table or the argument array. 

e.g. for a flat plate : 

RWTAB ( 1) = constant , 

(Make no entry for LRWTAB or XTABRW.) 


i 


\ 

i 
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TDKINP NAMELIST 


This group of data is entered only if INTDK = 1 in the preceding DATA 
namelist. Input follows the standard NAMELIST format. The first card of this 
group must contain $TDKINP starting in column 2. The data entries for the 
TDKINP namelist are punched out by the Two-Dimensional Kinetics Program (Reference 
1) . The last card of this group must contain $END starting in column 2. 

If the TDKINP namelist is input for a case, then LRWTAB and LPETAB in 
the DATA namelist must be input as the length of the arrays YITAB, PITAB, and 
XITAB (all are the same length). In addition, the arrays RWTAB, XTABRW, PETAB, 
XTABPE, and UETAB in the DAT’', namelist must not be input. 

The arrays in the TDKINP namelist, although punched out by TDK, are 
described below so that they may be edited by the user if desirad. 

Description 

YITAB (0 and PITAB (0 are the values of 
r^/x N and Pg/Pjf corres P° nc * in 9 to the value 
of x/xjj in XITAB (I) . 

These arrays , although punched out by 
TDK and contained in the TDKINP namelist, 
are not used by the program « They are read 
into a dummy storage area and later overlaid 
by other program arrays. Ineir values should 
be disregarded. 

The TDK punched output does not go all the way back to the chamber . 

The offset option in TDK (IOFF, see Section 6.5. 3. 6 in Reference 1) should be 
used to allow room for points in the pressure and nozzle radius tables upstream 
of the first TDK output point. 


Name 

YITAB 

PITAB 

XITAB 


ZMTABJ 
TITAB I 
CPITAB) 
VITAB I 
R0ITAB 
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i 

THERMO AND REACTANTS DATA FOR 0DE CALCULATIONS j 

The final gro of input data is the THERM0 data and REACTANTS data J 

for the One -Dimensional Equilibrium subroutines. This data is entered only if | 
IDEAL ■ 0 in the DATA namelist, that is, only for a hydrogen-oxygen equilibrium j 

case. The data cards in this group are fixed as to number, order, and content j 

and should not be altered; they are listed below. ! 

} 

The THERM0 data contains thermodynamic data for the six chemical 
species considered in the hydrogen-oxygen system, namely H, ^O, O 
OH, and Og. The REACTANTS data contains information pertaining to the chemical 
reactants hydrogen and oxygen. The REACTANTS data must be terminated by a 

single blank data card. ) 
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5. DESCRIPTION OF PROGRAM OUTPUT 

The various features of the output from the Mass Addition Boundary 
Layer Program are described here. The details of the output format can be 
seen in the sample case output. Portions of the output are described in the 
order in which they appear in a case. 

PRINTOUT OF PROGRAM INPUT 

The combined input from the DATA and TDKINP namelists is printed at 
the beginning of a case. All single-cell input variables appear on the first 
page, each in the format 

name = value. 

Where no value was input, the nominal or default value is printed. These 
variables are printed in functional groups, under the following headings: 

FLAGS AND OPTIONS 

PROBLEM LIMITS AND INITIAL VALUES 

REFERENCE QUANTITIES 

INPUT NORMALIZATION FACTORS 

EDGE QUANTITIES 

CONSTANTS 

CONVERGENCE AND EDGE CRITERIA 
COINTERS 

Tabular input data is printed on subsequent pages of the input printout. 
For each table, the tabie length is printed, followed by the dependent variable 
array and tha independent variable array in adjacent columns. The tables are 
grouped umier the following headings: 

STEPSIZE CONTROL TABLES 
WALL TABLES 

GEOMETRY AND EDGE TABLES 

(If the wall contour and pressure tables were entered 
in the TDKINP namelist, they will appear here as 
RWTAB, XTABRW, PETAB, and XTABPE.) 

EXPERIMENTAL PROFILES 

(These pages will appear only if experimental 
profiles were input.) 
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ONE-DIMEN SIONAL EQUILIBRIUM DATA 

Data for the reactants hydrogen and oxygen are printed if this is an 
equilibrium case, and the chemical species considered for the hydrogen-oxygen 
system (those in the THERM0 data) are listed. 

EDGE VELOCITY TABLE 

The edge velocity table (u vs x) , whether input or generated from the 

3u e 

input pressure table, is printed before generating the table of - ■- ■ ■ vs x. The 

velocities are printed in order from left to right under the heading EDGE VELOCITY; 

the corresponding x's are printed under the heading AXIAL DISTANCE. The velocity 

as printed is normalized by u f ; the axial distance is normalized by L. (Recall 

that the edge velocity u in the boundary layer solution is obtained by integrating 
du 

— =— — rather than by interpolating in the velocity table.) 





G-FUNCTION OUTPUT 

If a stretching function G has been specified by entering values for the 
input variables SN3 and GP0, data pertaining to this function will appear next in 
the output. For each initial boundary layer mesh point and for 30 additional mesh 
points beyond the edge (in case points are later added during the boundary layer 
solution) , the following quantities are printed in adjacent columns: 


y = , under the heading YTEL; 

y = G(50 , under the heading Y = G; 
G*(5) , under the heading GP; 

G”(y) , under the heading GPP* 

BOUNDARY LAYER PROFILES 




As explained in the Description of Program Input, profiles of various 
properties across the boundary layer are printed on a schedule governed by the 
input quantities NLPRNT and NSPRNT and at those axial distances in the input 
array XLIM. Starting with the wall point, every IYFRth point across the boundary 
layer will appear in the profiles. The edge point is always printed as well. 
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For every station at which the profiles are output, the following boundary 
layer quantities are printed above the profiles: 

STATION, X, S, DS, RW, THETAW, ZETA, ZETAP; 

Edge and wall conditions UE, TE, HE, PE, ME, RHOE, MUE, TW, MD0TW; 
Profile parameters DELTA*, THETA, TAUW, QW, CF, ST. 

These quantities also appear in the summary table at the end of a case and are 
explained in d 1 in the next section. 

Profiles printed under the control of NLPRNT appear first at a given station 
and are designated as Group I profiles. These profiles are printed in columns 
adjacent to one another; hence each row of output represents a single mesh point. 
The table below describes the Group I profile quantities printed for a given mesh 
point. Note that the temperature T is solved directly only at every IYEQth mesh 
point for an equilibrium case; temperatures at the intermediate mesh points are 
obtained by interpolation using the known points. 


Name 

Description 

Units 

N0. 

Mesh point # (wall = 1) . 


YEAR 

y , physical normal coordinate. 

ft. 

Y 

y, normalized coordinate. 


uAje 

S/u e 


H/HE 



R0/R0E 

p/Pe 


R0V 

DV = JY + P-'y?. 
P Pr“rC 


EPS 

Eddy viscosity (normalized) 



€=€0//^ 


T 

Temperature 

°R 


i 


i 


Profiles printed under the control of NSPRNT will follow the Group I 
profiles and are designated as Croup II profiles. These profiles will appear 
only for an equilibrium case. Like the Group I profiles, they are printed in 
columns adjacent to one another, so that each row of output represents a single 
mesh point. The table below describes the Group II profile quantities printed 
for a given mesh point. 


Name 

N0. 

YBAR 

0/F 

C(H) 

C(H2) 

C(H20) 

C(0) 

C(0H) 

C(02) 

MU 

PR 


Description Units 

Mesh point #. 

y, physical normal coordinate. ft. 

Oxidant/fuel ratio. 

Mass fraction of H. 

Mass fraction of H2. 

Mass fraction of H20. 

Mass fraction of 0. 

Mass fraction of OH. 

Mass fraction of 02. 

Normalized viscosity 

Prandtl number 


The species mass fractions, viscosity, and Prandtl number, like temperature, 
are solved directly only at every lYEQth mesh point. Viscosity and Prandtl 
number at the intermediate mesh points are obtained by interpolation using the 
known points. Values of zero are printed for the species mass fractions at the 
intermediate points; species mass fractions are not interpolated unless Lewis 
number and turbulent Lewis number are not unity, because the terms in which 
they appear in the boundary layer equations drop out when Le * 1 and Le^ = 1. 

Should points be added at the edge of the boundary layer during the solution, 
the species mass fractions at the added points, although non-zero, are meaningless, 
except at those points where they are solved directly (every lYEQth point, plus 
the edge point.) 
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SUMMARY TABLE 

A summary print of important boundary layer quantities at each station 
appears at the end of a case. Each page of the summary table contains data for 
thirteen stations. At the top of each page are three lines indicating the names 
and format of the variables printed for each station. The following table explains 
the quantities printed in the summary table: 


Name 

STATION 

X 

s 

RW 

THETAW 


Description 
Station number. 

Normalized axial distance x = x/L. 
Normalized wetted length s = s/L. 

Normalized wall radius, 

r -T / L 
w w 


Units 


Wall angle. 


K = tan 


-1 dr w 


radians 


DS 

A s, normalized integration stepsize. 


UE 

u e , edge velocity. 

ft/sec 

TE 

T , edge temperature 

°R 

HE 

h , edge enthalpy. 

ft® /sec 8 

ME 

Mach number at edge. 


PE 

P , edge pressure. 

lb f /ft® 

RH0E 

p , density at edge. 

lbjSec®/ft 4 

MUE 

viscosity at edge. 

lb^sec/ft® 

MDOTW 

Mass addition at wall. 

lbjSec/ft 3 


* w * w 

_ 1 , m . w m+ 1 v „ 

' 2 ( C m C m+1 > Wz 
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Description 


Units 


Name 

TW 

CF 

TAUW 

ST 

OW 

SUMQW 

DELTA* 

THETA 

ADDMAS 

ZETA 

ZETAP 

DUEDS 

RWPR 

XPR 



wall temperature. 

local shear stress coefficient. 


c I =, . / l»eV 

t ,, wall shear stress. 

w 

Stanton number. 


s« = -O w /o e iT e (H p - H w ) 


-q , heat transfer rate. 

W 

l (F/ Vis 

s init 

6*, displacement thickness. 
0, momentum thickness. 


°R 


lb f /ft 8 


ft-lbj/ft 2 sec 

ft-lbj/sec 

ft 

ft 


J m w ds ' of the 

s init 

mass addition correction to the 

displacement thickness. 

Value of the boundary layer normalization 
function £. 

Value of r‘ / the derivative of the boundary 
layer normalization function.. 

d(u /u ) 

/.i - , normalized edge velocity 
<WL ' derivative. 

Displacement thickness corrected nozzle 
radius, normalized by corrected throat 
radius. 

Axial distance normalized by 6* corrected 
throat radius. 



J 
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AFTER THE SUMMARY TABLE 

Following the summary table printout the program prints out the 
displacement thickness corrected throat radius and the thrust loss (in lbs) 
due to boundary layer effects. 

ERROR MESSAGES AND DIAGNOSTICS 

Following are the error messages and diagnostics which may be output 
from the Mass Addition Boundary Layer Program. 


( 1 ) 


( 2 ) 


(3) 


POINT WAS ADDED TO BOUNDARY LAYER IN SUBROUTINE 
AT STATION AND ITERATION 


MOMNTM] 
ENERGY 
ELEMTS 


This is a diagnostic message rather than an error message. It indicates 
that the edge criterion for velocity (subroutine MOMNTM) , enthalpy 
(subroutine ENERGY), or mass fraction of hydrogen (subroutine ELEMTS) 
has not been satisfied, resulting in the addition of a mesh point at the 
edge of the boundary layer. The criteria for adding points are given in 
the Description of Program Input under the input variables EPSLN1, 
EPSLN2, and EPSLN3. 


MOMNTM 

EXIT CALLED FROM SUBROUTINE ENERGY 

ELEMTS 

AND ITERATION 


AT STATION 


This fatal error message results (a) if six mesh points must be added 
to the boundary layer in a single iteration in the subroutine indicated, 
or (b) if another mesh point must be added when the maximum possible 
number of mesh points in the boundary layer (250) has been reached. 


XNTERP OUT OF RANGE 


X = 


X(l) = 
Y( 1) = 


X(N) 
., Y(N) 


This message is printed by the parabolic interpolation routine when an 
argument, given by X, exceeds the bounds of the independent variable 
array, given by X(l) and X(N) . The dependent variable values corresponding 
to X(l) and X(N) are given by Y(l) and Y(N) . Recall that averaged parabolic 
interpolation is used in the tables RWTAB vs XTABRW, PETAB vs XTABPE, 
and in the experimental profile tables UPR0F vs YBYNU, HPR0F vs YBYNH, 
and APR0F vs YBYNA. 
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6. OPERATING PROCEDURES 


This section contains information pertaining to the hardware and software 
requirements and the execution of the Mass Addition Boundary Layer Program. 

CORE REQUIREMENTS 

The UNIVAC 1108 version of the program requires overlay; with the 
cv - 1; v structure currently implemented using the Memory Allocation Processor 
and /vith the current program dimensions (250 points allowed in the boundary 
layer), the program requires 47800 words of core. 

FORTRAN LOGICAL UNITS AND I/O DEVICES 

The Mass Addition Boundary Layer Program accesses the following 1/0 
devices through their corresponding FORTRAN logical units: 


Logical Unit 

I/O Device 

Function 

5 

Card reader 

Input 

6 

Printer 

Output 

16 

User PCF 

Scratch 

17 

User scratch 

Temporary summary table storage 

29 

Last fifth of user 
drum 

Temporary THERM0 data storage 


If any of these logical unit assignments are non-standard at a particular UNIVAC 
1108 installation, appropriate changes to the system NTAB$ routine will be 
necessary. 

SPECIAL SYSTEM ROUTINES 

A system subroutine ERR0RX is called at the beginning of subroutine 
EXECUT to allow the summary table to be prlnt^gl before case termination in the 
event of a transfer to the system error exit routine MERR$ . Subroutine ERR0RX 
may not be available at other installations; if not, the call to ERR0RX may be 
removed, with the result that a negative argument tor the square root function, for 
example, will abort the run without a summary table print. This could be remedied 
by assigning logical unit 17 to a save tape and processing the tape independently in 
a manner similar to that in subroutine SUMTAB. 
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Execution 

The Mass Addition Boundary Layer Program may be executed from the 
absolute element on the program tape via the control card 

q XQT TBLABS 

The program may be executed from the relocatable elements on the 
program tape by referencing the processed MAP element via the control card 

I XQT MAPTBL 

O 

The input data deck which follows the XQT card will have the following 
general form: 

Case title card \ 


b$DATA 


DATA Namelist 

b$END j 

b$TDKINP | 

TDKINP Namelist) 

b$END I 

THERM0 DRUM 

Species thermodynamic 
data 


Always present 


Present only if INTDK = 1 in DATA namelist 


Present only if IDEAL = 0 in DATA namelist 


REACTANTS 

Reactants data 
Blank data card 

A perfect gas case may be run without overlay provided a dummy version 
of subroutine H00DE is used to override the version of H00DE on the program tape, 
as follows: 

SUBROUTINE H00DE (J) 

RETURN 


The program will then be loaded and executed, minus all of the One-Dimensional 
Equilibrium subroutines, via the control card 

l XQTTFCBL 
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Two sample cases are presented in order to facilitate program check-out. 

One case is a perfect gas solution while the other is the beginning of an Hg-Og 
equilibrium solution. Both cases are rocket engine flows and make use of the same 
TDK generated nozzle shape and pressure distribution. The input supplied with both 
cases has been set up to allow a complete solution to be obtained out to a normalized 
distance, x, of 21.39, however, for check-out purposes the perfect gas solution is 
terminated at XMAX = .052 and the equilibrium case is terminated after the first 
twenty steps. The length of the sample cases was restricted so that the program 
could be checked out without large expenditures of computer time and money. 

INPUT 

The input required to run the sample cases is contained in the fourth file of 
the tape containing the program itself. The output part of this section contains a 
copy of the input for each case as printed out by the program, so that the accuracy of 
the data punched from the tape can be verified. 

OUTPUT 

A sufficient sample of the output from the perfect gas case and essentially 
the complete output of the short equilibrium test case have been reproduced and are 
presented below. Both case*, use the same TDKINP Namelist Data, thus to avoid 
repetition the PETAB and RWTAB tables are presented only once, with the perfect gas 
sample output. 


PERFECT GAS SAMPLE CASE 


:-fcr;frrCT r,t«: ' A"Pi.E 


C'.br 


FlaRS :;rT!0' b 


mt*u * i 

l A -\P * 0 

J2'- • 1 

INTO* * 1 


(■1 PCR f’fc’ ,r e'C.T GAF » PT» MyrHOCFN-OXYCf. IG'l|Lt“RIU"> 
(»1 PC* L 4“ I*’A * Ptc*» «0 F"R T'MPliLiNT) 

<•1 POR 1'vfOMPRfc' S’-'LE FLO-f. P«» C0*PRF.S8!bll> 

(*1 PC*? AxlSY v ‘PT * GFCPITHY, «o POP T*0- , 'lME , ' ) SIC’-AL) 
(•1 IP JSP' T Taj; is CO*“t FRO* To* OUTPUT, *0 OThE*«IS£> 




PPO B LE« Ll* AM? INITIAL VALUES 


sjnit • .’.t^oro^p 
PtLTAI* 6.0099999-C& 


XIN1T * -3,5516000400 
2ETAFI* P.POOOnOP 


XMAX * 52POOOCO-0? nxi ■ 6,5499499.05 


RFFF«t£v.*. :<HAVTITIEP 

„RfcF ' J,u5025n0-0l 

iK-pt.-T mortal i zat i cm factors 
ys « 3,0502500-01 

irtF fiLANTITIFS 

i EOCP « O.CCOOOPO 

CCAS’iAkTS 

aFTRV8« 0,wC0P0?C 
PA*b • O.COOPOPO 
i.’tx » c.oPorocG 


URfcP » l.OOPOOOO*^ RHOREP* 3,0000000-03 8*UR£F» 5,0000002*0? 


VN 


• l,nOOOOOC*PO UE* * 1,0000000400 PIN ■ «,4996»00405 SNON • 4,4?60000*03 


PEOP.E « l,OQOOO0O*"O TEPOF » 6,58949*9*03 AFE06I* 1,3?1?00'-01 


TO 

m 

"U 

70 

o 

o 

c 

n 

r~~ 
— ■ l 

-< 

o 


— 1 

x 


o 

TO 

o 

z 

> 


PRI » 5,9999999-''! 

6P0 ■ 1,00000004*3 

RNDE* * P,COOOOO r 


P.AMbA • 1,1900000*00 

S\3 * 4, 0000001*00 


FMOLWT* l,4500nOP*01 
X8TAR • 3,000000'' 


plan • 1,0000000*00 > 

AFWALU* *9, 9990000403 Q 


CCNVER.PNCE AMO ED6F CRITERIA 

C0NVR6* 4,9999999.03 fcPSLNl* 3, nP0000P**'2 


FPSL*2» 3,000 n 000-02 ERSU3» 3,000000'’-P2 


CCu'TEhf 


I 


go 

T) 

O 

O 

73 


V AYJT • 1 

IY*, » 4 


NY | » I2v 


2" 


N8P*NT«9999 


ivpr » 




r MW * — * 






- 



STEPS I Z? CONTROL T4BUES 

IPXLl"* tO 

CXLIM XUN 


8 KTAB 


LSKTaB* « 

XTA88K 


2.73POOOO-C? 

2,5"0n0n0-02 

l.ijOCOOO-02 

2.0000009*03 

9.9999999- 03 

2.0000000- 03 

2.0000000- 02 

9.9999999- 03 
2. 0000900-01 
2.0000000-01 


-3,3994000*00 

-1,9120000*00 

-3,0000000-01 

-3,0000000-02 

0,0000000 

2.0000000-03 

3,2000000-02 

1.9120000*00 

3,0100000«00 

2 . 2000000*01 


1,1309000*00 

1,0300000*00 

1 . 0000000*00 

1 , 0100000*00 

1.1300000*00 

1.0000000*00 


• 3,39T400'-'*00 
-2,962400°*00 
3,2000000*02 
1,9120000*00 
3,0190000*00 
2 , 2000000*01 








wall tables 


ltwtab* 12 tW T *8' if 


TWTAB 

XTABTW 

8*DTAB 

XTABHD 

i.9oonooo*o3 

•3,5516000*00 

6,6000000-03 

-3,5516003*00 

1,9000000*33 

-1,4120000*00 

6,6000000-03 

-3,415000^*00 

t, 7000000*03 

•1.3660000*00 

6,6000000-03 

-3,1420000*00 

1,7000000*03 

0.0000000 

6,6000000-03 

•2.6690000*00 

1 , 9n00000*03 

0,0000000 

7,3000000-03 

-2,5950000*00 

1. >000000*03 

1.9120000*00 

4 , 6000000-03 

•2,3220000*00 

1.4400000*03 

1.9120000*00 

4,6900000-03 

-2,0490003*00 

1 . "!600000*03 

3.0050000*00 

5,9000000-03 

-1,7760000*00 

1,0*00000403 

4,1000000*00 

4,8000000-03 

"i,503000C*00 

1.0000000*03 

6,6300000*00 

6,6999999-03 

-1,2290000*00 

1.1000000*03 

8,4700000*00 

P, 9999999*03 

•9,9699999*01 

1.6200000*03 

2,1390000*01 

1,3600000-02 

-8,3300000*01 



1,4600000-02 

-6 , 9700000-01 



2,2600000-02 

•5,5999999*01 



2,7600000-02 

*4 , 2300000-01 



3,9700000-02 

-2,6700000-01 



3,6*00000-02 

0,9000000 



7,2200000-03 

4,1000000-02 



5,1299999-03 

1,3700000-01 



9,5999999-03 

2,6000000-01 



4 , 66000C9 ’03 

3,5500000-01 



2,2500003-05 

5,4600000-01 



2,3300. j’ *'3 

6,0999999-01 



1.97000C -f 

1,1340000*00 



7.999994V-04 

1,9650000*00 



0,0000000 

1,6300000*00 



0,0000000 

2,2000000*01 
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RWT4B 

m 

BRW 

9ETA9 

XTA99E 


1.7227*00*00 

-3,931 

4000*00 

9,4*99999*01 

• 3,99U009«00 


1, 7227*00*00 

•3,009000"*00 

9,9*40000*01 

-3,r**Q000*00 


1.7227*00*00 

-2,499 

000*00 

9,»999999*01 

-2,49*0000*00 


1.7227900*00 

•l,*7l 

200*00 

9,7479499.01 

•1,9719200*00 


1,721*200*00 

•1,*0* 

000*00 

• ,7*10900*01 

•1,*0*4000*00 


1.7l*««00*00 

•1,730 

700*00 

•,79*2299*01 

•1,790*700*00 


1,70*7*00*00 

•l.*92 

400*00 

9,7991199.01 

•1,4927400*00 

1 

1, ****700*00 

•1,374 

100*00 

9,7499300*01 

•1,974*100*00 


1, **14*00*00 

•1.49* 

•00*00 

9,7390900*01 

•l,49***00*00 

1 

1. **21*00*00 

-1.41* 

900*00 

9,7213300*01 

•1,41*9900*00 


l.«3S*300*00 

•1,341 

200*00 

9,7094900*01 

•1,3410200*00 


1. *11*000*00 

•1.2*5 

*00*00 

9,ft*09700*01 

•1,2*30*00*00 

.s' 

1,B7***00*00 

•1,1*9 

900*00 

9,4910*01*01 
9, *127*99*01 

•1,1*91900400 


1.3437*00*00 

•l,10» 

200*00 

•1,1072200*00 


1,902*200*00 

•1.02* 

900*00 

9,9*23200*01 

•1,0292900*00 

ifr 

1,49*1700*06 

•9.913 

300*01 

9,4949700*01 

•9,913*300*01 

s' 

1.403*700*00 

•*,734 

200*01 

9,3999799*01 

•*,7343200*01 


1,3490300*00 

•7,999 

200*01 


•7,9990200*01 


1,2701*00*00 

•7 i 179 

100*01 

•7,1797100*01 


1,2142400*00 

•*,39* 

099-01 

*••029999.01 

•4,3944099*0l 

!• 

tt 

1.1*07100*00 

•9, *17 

999.01 

••9111300*01 

•9,4170999*01 

1,11*9*00*00 

•4, *37 

000*01 

*,1919000*01 

•4, *97*000*01 

1.0*0*100*00 

1,091**00*00 

:5:«S 

900*01 

900*01 

7, *147099*01 
7,4131100*01 



1, 029*400*00 

•2,499 

*00-01 

*,9*00900*01 

•2,499M00*0l 


1,0132400*00 

•i,*7l 

700*01 

*, 4911400*01 
9,9799099*01 

•1, *71*700*01 

i 

1.0033*00*00 

•*,43* 

799.02 

•*, 43*9799*02 


•3,37* 

000*02 

9,3902*99*01 

•3,9749000*02 

t 

1.0000100*00 

•3, *32 

100*03 

4,4999000*01 

•9,*32*100*03 

U 

ti 

1.0000000*00 

0,000 

600 

4,429*900*01 

0,0000000 

i 

1,0000100*00 

*,479 

399.04 

4,2214000*01 

9,4799399*04 

1.0000400*00 

1,92* 

900*03 

4,0343900*01 

1,92*7900*03 

* 

1,0000*00*00 

2,93* 

900*03 

3, *42*700*01 

3,93*4900*03 


1,0001*00*00 

3,970 

900*03 

3,**00700*01 

9,9702900*03 


1,0002900*00 

»,033 

399.03 

3,4777900*01 

9,0332999*03 

i 

1,0003*00*00 

♦ ,ljl 

099.03 


4.1212099*03 

P 

1,0009300*00 

7,232 

099.03 

3,199*300*01 


, 

1,0007100*00 

*,3*7 

799.03 

2,9711400*01 

«, 3*79799*03 

‘i 

1.000*200*00 

9,924 

•00*03 

2, *109700*01 

9,9249*00*03 

<0 

1,0011*00*00 

1,070 

700*02 

2 ,994*900*01 

1,0704700*02 


1.0014400*09 

1.190 

700*02 

3.9024200*01 

1,190*700*02 

n 

1.0017900*00 

1,313 

•00*02 


1,9130*00*02 


1 . 0021100*06 

1,*39 

900*02 

2,2099000*01 

1,437*900*02 


1.0029100*00 

1.9*4 

900*02 


1,9*44900*02 

S- 

1,002**00*00 

1,499 

900*02 

1,9927000*01 

1, *934900*02 

i 

1,0034900*00 

1,*24 

*00*02 

1,9001900*01 

1,*24**00*02 

iSr 

w 

1,003**00*00 

1,99* 

000.02 

1, *714*00-01 

1,99*1000.02 

r. 

1.004*000*00 

2.093 

•00*02 

1.94*7400*01 

2,0937*00*02 


1.0092*00*00 

2,231 

100*02 

1,42*0*00*01 

2,2917100*02 

■f 

1, 0099*00*00 

2,371 

100*02 

1,3094900*01 

2,9719109*02 


1.00*7*00*00 

2,914 

000*03 

1.19*9900*01 

2,91*4000*02 


if 






RWT4B 

XT4BP* 

f> ft*b 

XT*3 p t 


1.0076*00*00 

2,6*00700-02 

1,1233900-01 

2,6600790-02 

V /O » 

1.0065300*00 

2,605*200-02 

1,1246900-01 

2,605*200-02 

1.006000*00 

2,64*1*00-02 

1,126*900-01 

2,64*1400-02 

rn 

1.0103100*00 

3, 0*30000-02 

1,1262900-01 

3,0*30000-02 

"O 

■ o . 

1 , 0111600*00 

3,23*6700-02 

1,1266600-01 

3,23**700-02 

1.0120700*00 

3,3*0**00-02 

1.131*500-01 

3,3604600-02 

1.0126300*00 

3,5235500-02 

1.1337*00-01 

3.5235500-02 

•' O 1 

1.0136300*00 

3,6**2196-02 

1,135*500-01 

3, *6*2199-02 

C 

1.01*7000*00 

1.0161000*00 

3,80*7*00-02 

4,03*6500-02 

1.1375000-01 

1,1224700-01 

3, *0*7*00-02 
4,03*6500-02 
*. 3714*66-02 

0 

1.01*1300*00 

4,371*996-02 

1,111*500-01 

1.0212200*00 

4,6732300-02 

1,1010700-01 

4, *732300-02 

T 

1.0233600*00 

3,5664200-02 

1,1101600-01 

5, §6*4200-02 

— S 

1.0324400*00 

6,7034666-02 

1,0*60300-01 

*,7034*99-02 

X 

1,0*24600*00 

8,33*7*66-02 

1, 0*0*600-01 
1,0756300-01 

*,33*7*96-02 

o 

1,0374600*00 

1.07*1100-01 

1,0761100-01 

1,0765600*00 

1,4366700-01 

1,0647300-01 

1. *3**700-01 

-n 

1.1123300*00 

1.6742100-01 

1,040*100-01 

1,6742100-01 


1.15*5100*00 

2,6632000-01 

1.0110100-01 

2, *652000-01 

1,1*67600*00 

3,4015400-01 

6.61160*6-02 

3,4015400*01 

1,2424400*00 

1.2642500*00 

4,0674600-01 

4,77*6800-01 

9,4811000-02 

6,2007469-02 

4,0674*00-01 

4,776**00-01 

rr ' 

o 

1,3255600*00 

5,433*000-01 

*,*221766-02 

5, *536000-01 

1,3*72100*00 

6,1326*00-01 

6,5560000-02 

*,132*600-01 

73 

1,4067100*00 

6.6066000-01 

6,23*1369-02 

*,*0*9000-01 

o 

1.4466600*00 

7,4617166-01 

7,634*666-02 

7,4617166-01 
*, 1511600-01 


1.4606200*00 

6,1511600-01 

7,636*669-02 

z 

1,531*500*00 

6, 61*0100-01 
*,5066766-01 

7,36057*6-02 

8,61*0100-01 

> 

1,5741500*00 

7,0367069-02 

*, 50667*6-01 

- r - 

1,41*2*00*00 

1,01*6100*00 

6,7442600-Ul 

1,016*109400 

"'O 

1, *563100*00 

1,0*01700*00 

6,4403769-02 

1,0*01700400 

1,702*100*00 

1,1606400*00 

*, 1452*66-. 2 

1,1606400400 

> 

1,74*7300*00 

1,2330100*00 

5.6513366-02 

1.2330100400 

; O 

1,7613600*00 

1,63*6600*00 

1,3056*00*00 

1,3602000*00 


1,305**00400 

1,3*02000400 

m 

1,6630300*00 

1,4551100*00 

4,6*67*00-02 

1,4551100*00 

Th 

1.6302400*00 

1,5317500*00 

4,7366600-02 

If 5317500*00 

1.6771300*00 

1,6076400*00 

4,5126600-02 

1,6076600400 

"O 

2.02*6190*00 

1,6*54200*00 

4,2*67600-02 

1,4664200*00 

O 

2,0*13000*00 

1,777*200*00 

4,0260400-02 

3,6*77200-02 

1,7774200*00 


2,1264630*00 

1,65*2900*00 

1,6542600400 

2.153*100*00 

1.6411100*00 

3. 7236600-02 

1,6411100*00 

2.231*600*00 

2,0259700*00 

3.6063599-02 

2,0259700*00 


2,2966500*00 

2.1357400*00 

3.3635*00-02 

2,1557400*00 


2,33**000*00 

2,2042*00*00 

3,3157100-02 

2,2042*00*00 


2.3673300*00 

2,30*6100*00 

3,150*300-02 

2,30**100*00 


2.4**1200*00 

2,3674600*00 

3,0393000-02 

2,3674*00*00 


2.50*0700*00 

2,499*600-00 

2.6163400-02 

2,7723500-07 

2,49***00*00 


2.372*000*00 

2. *165900*00 

2, *1*5*00*00 


2.727*600*00 

2,6030000*00 

2.4636200-02 

2,6030000*00 


2.6317100*00 

3,1379500*00 

2,2967100-02 

3,1376500*00 


3.0073400*00 

3,4414*00*00 

2,0502400-02 

3,441**00*00 


3.1366300*00 

3,749**00*00 

1.6437200-J2 

3,74*4*00*00 


3,3369*00*00 

4,12*2000*00 

1,6650*00-02 

4,1242000*00 



St r ■■* It* r^fliwwwwi; 



RWT4B 

XTABRtt 

PETAB 

xtYspe 

3, 

,9199000*00 

4,9195700*00 

1,4939400*02 

4,9199700+00 

3,7300*00+00 

4,9961600+00 

1.3074100*02 

4,9941400*00 

3 

9921000*00 

9,9241100*00 

1,1993700*02 

5,9241100*00 

4 

,2'»17700*00 

6,1469200*00 

1,0110900*02 

4,1444200+00- 

4 

4 

4*0*000+00 

7932200+00 

6, 6346800+00 
7,4993299*00 

4,9749300*03 

7,4111999*03 

4,9344900+00 

7,4993299*00 


,C*05*. 00 

*,9663900*00 

6,4753799*03 

*,9**3900+00 

* 

2201100+uS 

9,0992499+00 

3. 9974600*03 

9,0992499*00 


<■33400*00 

9,4414900+00 

5,6111900*03 

9,4414900*00 


•474199*00 

1,0292400+01 

9,1299999*03 

1,02*2400+01 


>431300*00 

1,0997400+01 

4,6366300*03 

4,2604099*03 

1,09*7400*01 


*427*00+00 

1,1499000+01 

1,16«*000+01 

*. 

,1193200*00 

1,2399900+01 

3,9794200*03 

1,2399*00*01 

6, 

3021900+00 

1 .3179100*01 

3.6994400*03 

1,3179100+01 

\ 

,4* n +000+00 

1,4009900+01 

..3569700*03 

1,4009900+01 

s 

64. 1900+00 

1,4912900+01 

3, 0*94400*03 

1,4912900+01 

6, 

,9643300 >00 

I«»9l9l00+01 

2,9342900*03 

1,9*1*100+01 

6, 

,9607100+00 

1,4319400+01 

2.7179100*03 

1,6319400+01 

7, 

,0926600+00 

1,6902900+01 

2,4197000*03 

2,9339900*03 

1,6902900+01 

7, 

,1429700+00 

U' 299700+01 

1,72*9700+01 

7, 

,2301000+00 

1,7797900+01 

2,4999900*03 

1,77*7*00+01 

7 , 

,3191600+00 

1,9307900+01 

2.39934*9*03 

1,4307*00+01 

?: 

,4091300+00 

,4933399*00 

1,9929400+01 

1.9397300+01 

2.3214700*03 

2,2940900*03 

1,9*29400+01 

1,9347300+01 

7, 

,9792199+00 

1,9949900+01 

2,1979200*03 

1,9949900+01 

7, 

,6969699*00 

2,0497000*01 

2,1941300*03 

2,04*7000+01 

7, 

,7339100400 

2,1099000+01 

2,1190900*03 

2,109*000+01 

7, 

,7769600+00 

2,1390900+01 

2,0733900*03 

2,13*0*00+01 
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H 
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VELOCITY TABLE CE‘'E l 'ATEr 


EDwE VELOCITY 

• 368458-00 

• U9037-01 
,166*16+01 
. 397392+01 
.633122+01 
,74*4SB+oi 
,66*786+01 
,*ll949+oi 
,9 1 754S*C1 
. 944934*01 
,*95u3+0l 
.105100*0? 
,19**62+02 

»U‘56C*0? 

«1257t3*C2 

.139*04+02 

,13?6**+0? 

AXIAL DISTANCE 

■. 33516^*01 
■.1496*8*01 
■••73432-00 
*.2499*8»00 
,293649.0? 
,U9o67*Cl 


t 4oiA74*»oO 

.1120*3+91 

,l8i43p*ji 

• 434421*ij1 
,647826* ji 
,764007*'’! 
.66676P' . 
,9n654*i i 
,916343*01 
,*5 1 *22+oi 
.1001*7*02 
.105730*02 
.11000**02 
.U764i*n2 
,1J4?25*o2 
,l3ll67*..<2 
,132«64*o2 


.223i7i*oi 

.336348*01 

,670349*01 

,477966*00 

.1C1»*1*01 

.16ci7«4*C1 

,230681*01 

.412420*31 

.*0*625+01 

.149125*32 

.1*3673*02 


-,3oc500*ni 
*, 1*1*93*01 
•,7953-2-po 

*.1671*7-90 
.397023 -., 2 

,131306*01 

.237i9i-oi 

.362353-01 

,«338’9*ni 

,5453 - ’0* , l0 

.lt)9ci7*ol 

,16*342*91 

,239746*01 

,43i5S7* 0 i 

,964145*01 

.1381*1*02 

,1*9495*02 


,673692-00 

,116494*01 

,210656*01 

,472333*01 
.462772*01 
,778649*01 
, *02o39*01 
,*11362*01 
,*20733*01 
,938602*01 
,100665*02 
,104563*02 
,110479*02 
,116*27*02 
,126922*02 
,131474*02 
,132**1*02 


». 243*00*01 
•.134102*01 
'.7i757i«oo 
*, 843858*01 
.303324*02 
,143765*01 
,251440-01 
,366622*01 
.107*11*00 
»6i32***00 

,24*966*01 
,49*616*01 
.102*24*02 
• 1.63l54*o2 
.<04970*02 


,1P3003*01 
. 1?0*7?*01 

,236626*01 

,5l566?*oi 

.677277*01 

,7934.7+ni 

,9l326r*01 

,*11077*91 

,*22533*01 

,964629*01 

,101563*02 

,106*62*02 

,111076*02 

.120053*02 

,1276674 0 2 

,13i7i8*o2 

,133106*02 


,106448+01 

,2661*4*01 
,582400*01 
,6*1724*01 
,80"305*01 
,*l3027*01 
,*13403*01 
,*2*3*3*01 
,970489*01 
,1022*0*02 
,107552*02 
,112410*02 
. 1212*8*02 
,1262*2*02 
,131*39*02 
,133234*0? 


,10504**91 
,133386*91 
,2*8764*«h 
,3B7742*6i 
,7o6l36*ri 
,8?3363*ri 
,*12746*01 
,*13106*31 
,931323*91 
, *76342*01 
,102**8*12 
,107*3**92 
,1133*1*02 

,122447*02 

,128623*02 

,132i34*-2 


’.1671*2*01 

■,l263o6*ni 

>,63*641*00 

■,337636*01 

.612121*02 

,156443-91 

,2640fl7*01 

,38o874*oi 

.143**7-00 

,68o**P«90 

.123301*«’1 

.163*2**01 

,2*l65**Ol 

,552411*01 

.10» 9 74*02 

,1*802**02 

.210560*02 


1*0660*01 

ll85l»*01 

561710*00 

3632*1*02 

7232*1*02 

16*345*01 

2*05*2*01 

4g36*5*01 

1*7421*00 

74*172-00 

130586*01 

1*4111*01 

2*0300*01 

*146*2*01 

itittKg 

213*06*02 


»ii73o87*9t 
• , 110722*C1 

*,463760*00 

: jsbb .„ 

,182466*01 

,2*4*l4*oi 

,437i5o*01 

,26*320*00 

•8l3ll**00 

.13*020*01 

,2025*7*91 

,3137*5«01 

,663468*01 

,123*59*02 

,177678*02 


,10*185*01 
,141*67*01 
,330561*01 
,603420*01 
,720545*01 
,836414*01 
,*12500*01 
,9i675n*oi 
, *34337*01 
,* 82026*01 
,10373**0? 
.1067***02 
,114407*0* 
,123756*02 
,l2*4i8*o* 
,132335*02 


■,145274*01 

•.102*2**01 

•,4o5*4**oo 

, *47**4«03 
,*52466*02 
,l*98l0»01 
.3g*3og*!n 
,417323*01 
,340154*0- 
,881601*09 
,145511*01 
,213574*01 
, 344146*01 
,769*33*01 

•i2i 7 *ito? 

,1*3o7**02 


,107822*01 

.152*17*01 

,3*3*58*01 

,*17*63*91 

,734*76*01 

,69407**01 

,912238*01 

,*15320*01 

,94201**01 

,*6*741*01 

,10*455*02 

.108*72*02 

,11591**02 

,125084*02 

,13002**02 

,132511*02 


•.197481*01 
■•*5l3*3>00 
’.327*i*-00 
,1*2875-02 
,107047-01 
,20*378*01 
.323**7-01 
•598842*01 
, 4o*74*-9o 
« *90**8*00 
.153179*01 
,220628*01 
, 374*46*01 
,698639*01 

:ias %8 






y 


/ 
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^fc«FECT GAS SA*PU C^SE 


STATION 


«3tSbt60C0 4 f j0 


iinonoooc 


ns 

6,54999 J9-n5 


1#72279 qa^ 0 q 


T^FTAa 

0100:000''' 


7ET A 

1«’12«US-C4 


ZETAP 

4,6821*3*-02 


POC-t *\D -ALL CONDITIONS 


P*»OP!lE PA^ETER*. 



CF 

■ 

4*Hn*; 3 >!F 


* 9 7in7^5^»rfj > 


DELTA* « 

.1909092-04 

CP 

B 

.3630303 


tr 

■ t 65A64394-r.4 mHuE 

« # i?6a56**ru 


THETA • 

,1745367-04 

,5118137*91 

ST 

B 

,31*4418 


*F 

* .1414A28*; 1 * nlJfc 


« ,l7ftft391«n5 


?*U- » 





PF 

» .4446301-06 TW 


• ,l90000 r >f)4 
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l3o97«9* 0 i 

■oft 

1 ***. 

49 

If 
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4 

6,62B9758«o9 

2,7736664-02 

1,5372457-04 
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3157477*01 
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4*24614*01 

5373366*01 

5976896*01 

6666576*01 

82227*4*01 

0006*9g*01 

1010635*01 

261*295*01 

2367110*01 
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3.595nn57-i2 

5.675*654-11 

6.o33*7o7-io 
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2.69463o6-0* 

1.32*2*22-07 

5.6o*6o44«o7 

2.q665**6-o6 

6,9759*69-06 

2,il9*8S7-o5 

5.9246656-05 

6, 4826131-0* 

7.H02969-03 

1.2*56345-02 

2,242*9*2-02 

3,69*o4o6-o2 

5,6363261-02 

6,6465362-02 

2,5945789-01 

3,542*029-01 

4,74*21*3-01 

6.2602347-01 

6,1232335-01 
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2,9132573-d 

4 

2,6o94d3l-o7 

2,7736664-02 

7 

7,8668856-07 

5,5473326-02 

6.5232934-02 

3,072*053-01 

10 

1,7147849-06 

«.320999i-n2 

1,3588o89-oi 

3,3383472-01 

13 

3,2119251-06 

1. 1094666-01 

2,l3o53o7-oi 

3,6793720-01 

1* 

3.4741129-06 

1,3866332-01 

2,78o9369-oi 

*.0212528-01 

1* 

8,727o368-q6 

1,6641996-01 

3,2934483-01 

*.3lJl2o’-01 

*.559i4i9-oi 

22 

1,3226622-oS 

1 «-■'« 1 

3,7o26554-oi 

25 

1,9239631-03 

.1 

*,0*19693-01 

*,756q675-01 

2* 

2,?l437S5-o5 

2.4962997*01 

*,33559n7-oi 

*,9199563-01 

31 

3,7237ip7- 0 3 

2,7736664-31 

*,5994417-01 

5,060*986-01 

?4 

4,992?i6i-o5 

3,0310330*01 

4,8432449-01 

5,1*67213-01 

37 

6,5627260-05 

3,3283997-01 

5,rj7245ro*01 

5.31232*8-01 

40 

8,4819998-05 

3.6n57683-0l 

5,2698929-01 

5, *112433-01 

43 

1,0801320-0* ' 

3.8831329-01 

5,497o65o-oi 

5,5155949-01 

4A 

1,3576966-04 

4.1604996-01 

5,6948924-01 
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49 

1,687o923-04 

4,4378663-01 
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9,Jl522l8-01 

5,6ll7»4/-oi 
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4,7i52329-01 

6,0652754-01 

55 

2,3294027-04 

4,9925996-qI 
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5,9o679p9-oi 

5» 

3,03*2668-04 
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6,40392x4-01 
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61 

3.67io932-04 
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6,566x359-01 
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6,869o56x-oi 
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5,l9i8969-04 
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6,27788i9-oi 
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6.3794327-01 
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73 
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6.6347993-ox 

7,x5o437i-01 

6.46252*4-01 

6,5557543-01 
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6,4l3o599-o* 

6,934166q-01 

7,2842l98-oi 

79 

9,8052400-04 

7,2115326-01 
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62 

1.1392531-03 

7,4888992*01 

7,5381631-01 

7,6638363-01 
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7,7662639-01 

6,8472534-01 

66 

1.3263734-03 

8.0436325-01 

7,7993433-01 

6,9630758-01 

9 1 
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8.3209991-01 

6.5963636-01 

7,9474302-01 

J.0* 7 »l43-01 

94 

2,0299409-03 

•»10829ll-01 

7,2532769-01 

97 

2.3366943-03 

8,8737325-01 

8,281*463*01 

7,4333008-01 

ICO 

2,6878519-03 

9.1930991-01 

8.4680766-01 

7.6411056-01 

103 

3, 0*06626-03 

9.4304658-01 

8,66785*6-01 

7, 8SiB6o7-oi 

1C* 

3,5538410-03 

9,7j78324-ci 

8,8846842-01 

•, 1649*06-01 

10* 

4,0679197-03 

9,985i99q-qx 

9,1270878-01 

9,5o88o97*q1 

11? 

4,70569*2-03 

1.0262566*00 

».4077ii5-01 

8,9417867-01 

113 

5,4228163-03 

1.0539932*00 

9, 7i86666-ol 

9, «84*l79»oi 

11* 

6,2584943-03 

1. 0*19299*00 

9,9532393-01 

9,6996689-01 

12! 

7.2364914-03 

1,1094666*00 

1.0000194*00 

1.00002*0*00 

12* 

8,3863570-03 

l,l372o32*C.1 

1,0000000*09 

1.0900000*00 
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tmftak 
■A i 7 lAo 8 jj 3*32 


TETa 

1, < >77c2i''-02 


zetap 

•2»i283o59-g3 


profile parameters 
OElTa* ■ 1 77630 »C -P 3 C e 

THETA ■ i lDl 73 l*i -02 ST 

TAU W » , 328341^*03 

UW ■ , 3287923*07 


,210’673-02 

, 8126967-03 


RO/ROE 

5 ? 799 ?l*oO 

4323839*00 

2343532*00 

9953148*00 

7i77 0 73* D o 

4867872*00 

J&ti&S 

1023774*00 

0323385*00 

9758946*00 

*2*0002*00 

8859637*00 

84aon4l*oo 

8i3Q410*or 

78o4l64*oo 

7*97x34*00 

7206389*00 

6929666*00 

6665q79*oo 

6411013*00 

5iSS8l?:gS 

5698542*00 

5473820*00 

5233775*00 

5o37446*q0 

48242q3*00 

46o43’6*oo 

436i469*oo 

4o**44o*00 

3786872*00 

3452974*00 

3087111*00 

2687359*00 

2247411*00 

1752525*00 

ll83447*oo 

0365443*00 

0101348*00 

9997297-01 

0000000*00 


RCV 

2,60669* 7*oo 
2,8o66oo5*oo 
2,8o59io6*oo 
2,8033574*00 
2,7986639*00 
2,79o86P5*00 
2,77954i4*oo 
2,76367io*oo 
2,7420070*00 
2,7i312««*0C 
2,67548i8*oo 
2,6274239*00 

2,567274 0 * 00 

2,4933507*00 
2,4039810*00 
2,2974845*00 
2,1721553*00 
2,0262497*00 
1,6579832*00 
1,6655221*00 
1, 44696?2*qo 
1,2003018*00 
9.2340121-01 
6,1393965-01 
2.6936267-01 
1,13229^4-01 
*,3717996-01 
1,0083434*00 
1.32513*0*00 

2 , o97689i*Qg 

2.72774p3*00 

3,4i96768*oo 

4,17879-6*00 

5,0113945*00 

5,92476i7*oo 

6,92642*0*00 

8,0221414*00 

*,2l36ii5*00 

I,38l52 ri 4*0i 

1.6010734*01 


EPS 

coooooo 

1255540-04 

C2l5io6-o3 

7odi2l0-02 

B290284-01 

0885532-01 

64373 * 8*00 

5*35652*00 

1931652*00 

7309301*00 

4303696*01 

9937797*01 

6820334*01 

7341284*01 

i63737q*oi 

8413388*01 

0793353*02 

3047337*02 

5632155*02 

8577651*02 

l9t5oo4*o2 

5676498*02 

9897785*02 

46l8i#4*o2 

9434353*02 

5354768*02 

7919462*02 

7784129*02 

77oi66i*o2 

6106575*02 

4n36466*02 

9937235*02 

3o»«86i*o« 

3800*32*02 

4325802*02 

6032429*01 

5589660*01 

4342751*01 


1842,9 

2002.3 
220»,9 
24U.6 

2589.3 

2734.4 

2852.5 
295C.8 

3033.4 

3110.8 

3180.1 

3245.4 

33 0 8. 0 

3368.6 
34*7,6 

3415.7 

3542.6 i 
3396 9 

3654.6 

SU:l 

3620.5 

3676, g 

3931.9 

$}:• 

4256.8 

4350.2 

4458.2 

4697.0 

3103.2 
338*,9 
3676, 
3937, 4j 
3997,71 
3997,81 
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,lt - 




THETA* 

5 , 4979727-01 


ZE TA 

3 ,** 42359-02 


PR0PU.E PARAMETERS 
OEUTA* ■ , 1401**4-02 CP 

THETA » ,1791396-02 ST 

TAyW * ,420*545403 

OW ■ ,2*99923407 


PO/ROE 

4l47i444 00 

JiSHMtSS 

1357550*00 

0131370*00 

*243230*00 

0401493*00 

033 lo 54 *oo 

# 100 * 22*00 

903540**00 

7*05*22*00 

7*»337**oo 

75074*5*00 

722*379*00 

7o9*g59*oo 

*9357s2*00 

*75*393*00 

•55*n3t*00 

#3*2701*00 

0113**2*00 

S07«269* O o 

5*35254*00 

5395*19*00 

4&2*Ul*00 

47g3**2*00 

44*2o*2*00 

425*930*00 

40001***00 

3723943*00 

5 ; s ? s ;:: g ? 

2*5*024*00 

222**74*00 

j ;« 8 ? s:ss 

0570727*00 

masss : 

0000000*00 


ROV 

i*3j0430-ol 

,*343715-01 

,«30259*> o i 

i*46354*-qi 

i*5054i2"oi 

,6779967-qi 

, 7 o 747 o 4 -oi 

i74*5q34*01 

«007«H»01 

•« 12255 -oi 

i 9235 o 32-01 

97*9 1 9*- 01 

005l59**oi 

,9040*37-01 

,* 7924 * 9 -oi 

*43'9*1—01 

, 2 i 77 q * 9 - 0 i 

529*§52-oi 

02*99n*-o2 

3403524-02 

»397i?9-oi 

,242»l97-oi 

192*606-01 

996*155*00 

5o34*o7*oo 

q75*73**00 

7202512*00 

443140**00 

1*5*319*00 

20233*3*00 

50*0654*00 

0*20960*01 

2250470*01 

92337n4*01 


looooooo 

,2l2l479- 0 3 
1 105*501-02 
6297*76- 0 i 
244*53o*or 
,*00**30*00 
140*302*00 
540o*3q*o1 
39i93»9*oi 
*44*430*01 
*0*1*51*01 

i00-**3**02 

2*20402*02 

*335107*02 

0*2*030*02 

ssasitjs 

53*3*0^*02 

45*5352*02 

09133***02 

*443020*02 

tnsa% 

3«045i3«o3 

ssma 

47*9«5i*03 

443*924*03 

43*Q4l**o3 

370*240*03 

2454449-03 

032*1*0*03 

4*0*402*02 

54*77*3*02 

VSMS SI 

333 0 1?I*01 
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PERFECT 3AS ^a*PLF CASE 


STATION 

Q < 

Ufc 

X 

THETA* 

Tc 

S 

ns 


.0 


3,6846-?? 

-3.55i«-oc 

o.on?c 

6.56*8*03 

o.aooo 

6,55Q0*n& 

1.4144^06 

1 1 J 

1 • 7228*00 

3.6846*02 

•3.5315*90 

0*0000 

6.3668*03 

6. 5500-05 

6.5500^5 

1.4146*06 

2.0 

1.7?2»^00 

3.6647*02 

*3.33i3*00 

0*0000 

6,5868*03 

i.4082-0* 

7,3325-0* 

1.4146*08 

3.C 

1.722«*00 

3.6847*02 

•3,5314*00 

o*ooon 

6,5666*03 

2.2743-04 

6,6624-0^ 

1.4l46'-08 

4.0 

1. 7228400 

3.6848*02 

-3.3313*00 

8:82S9- o5 

6.3866*03 

3.2707-04 

1.4146*08 

9.0 

1.7228*oo 

3.6849*02 

«3,35i2*C0 

0,0000 ' 

6,3868-03 

4,4113-04 

1.1438-04 

1.4146*08 

6,9 

1.7228*00 

3,8649*02 

-3,53io*00 

0.0000 

6.3868*03 

3,7337-04 

J..3: T 4-04 

1.4146*08 

7,0 

1,^228*00 

3.6850*02 

•3.5i05*00 

0.0000 

l.»l3i-o4 

6.5866*03 

7,2488-04 

1.4146*0« 

6,0 

1,7228*00 

3.6832*02 

•3.3So7*00 

3.0090 

6.5868*03 

8.9*11-04 

1,7423-04 

1.4146*06 

9.0 

1,7228*<ij 

3.-«53*02 

-3.55o’*00 

0.0000 

6.5868*03 

1.0*93-03 

2.0037-04 

1,4146*08 

10.0 

1.7228*00 

3.6854*02 

•3.55o3*00 

o.ooro 

6.5668*03 

1,3299-03 

2.3c •'04 

1.4146*08 

11.0 

1.7i^t*oo 

3.6855*02 


0.0000 

2.64*8-04 

6.5866*03 

1.4146*06 

12*: 

1.7r2**9L 

1 156*P2 

•3,*>497^iC 

0 . 0 u 0 11 

6.5868*03 

1.69?6-g3 

3.Q473'-M 

!, 4146*0* 


“E 


HUE 


CP 

PE 


. "t>Tv 



RHCE 


T» 


mmmm 

7.1071-02 

1 

7984-r 6 

3t 

83o'-r3 

4.4*63*05 

0 

O-v.r 

5. 

1183-03 

1.9686-02 

1 

9oc«*o3 

:• 

ecor 

7.1071-02 

1 

78*4-06 

3. 

4' - 33 

4.4463*05 

2 

9581-02 

4. 

6*i- roo 

1,9686-02 

1 

90CP*ft3 

c * 

oroo 

7.1072-02 

1 

7894*06 

3. 

6796-H3 

4.4463*05 

2 

9601-02 

4. 

9l72*00 

1,9686-02 

1 

9000*03 

0. 

onoo 

7.1073-02 

1 

7484-06 

3. 

7485-03 

4.4463*05 

2 

9569-02 

9. 

0095*00 

1.9686-02 

1 

9000*03 

c. 

0003 

7.1075-02 

1 

7884-06 

4. 

3059-33 

4,4463*05 

2 

9544-02 

9. 

7546*03 

1,9686-02 

1 

9000*03 

Ct 

0000 

7.1076-02 

1 

7884-06 

4. 

6646-03 

4,4463-05 

2 

9542-02 


2342*00 

1,9686-02 

1 

9000*03 

0. 

0000 

7.1078-02 

1 

7884-06 

9. 

1787- n 

4. 4 163*05 

2 

9542-02 

6, 

*21**0J 

1.9686-02 

1 

9oc9*t3 

0. 

0000 

7,1080-92 

X 

7864*06 

9. 

2888-03 

4.4463*05 

2 

9546-02 

7. 

0424*00 

1,9686-02 

1 

9co"*n3 

c. 

0000 

;. 8082*02 

1 

7864-06 

b. 

5^75-03 

4.4463*05 

2 

9563-02 

7. 

T5ii*00 

1.9686-02 

i 

9ncc*03 

0. 

oooc 

7,10®4-(’2 

1 

7864-q* 

5* 

5513-02 

4.4463*05 

2 

9573-02 

7. 

42o®*00 

1,9686- 02 

1 

9nc0*03 

0. 

0000 

7,1087-02 

1 

7884*06 

9 ( 

88q9-p3 

4,4463*05 

2 

9565-32 

7. 

•*22*00 

1 • 9686-02 

1 

9fiao*o3 

D« 

0000 

7. U *0*02 

1 

7864-n* 

6. 

0*6o-c3 

4.4463*05 

2 

9555-02 

4. 

0*37*00 

1,9686-02 

1 

9oon*o3 

0* 

OPOO 

7 . 1094-02 

1 

7*84-06 

6 « 

3241-53 

4.4463*05 

2 

955n*02 

e. 

4564*0C 

1,9666-02 

1 

9ooO**3 

0. 

0O0P 


ST 


OECTA* 


zeta 

G- 


TWET* 


ZETA? 

SUC.* 


ADO'US 


DUEDS 

3.1944-P3 

1 

,9o9i*CiS 

I® 

9l2«-0* 

2,3338*06 

1 

,745«-"5 

4. 

*82^-r^ 

O.OCO" 

0 

.0000 

-Oi 

000' 

3«o9i3»03 

1 

,9377-05 

2% 

0569-04 

2.2583*06 

1 

, 7958-05 

2* 

2061-01 

1.5146*02 

1 

.1112-0* 

6 • 

4362-02 

3,l7ij3-03 

1 

,0532-05 

2. 

2497-04 

2.3163*06 

1 

B4io-05 

2« 

5593-01 

3,2498*02 

2 

,3890-06 

6. 

4392-02 

3,3082-03 

1 

.9014-05 

2. 

39on*f)4 

2,4l7l*o6 

1 

,83*3-05 

1 8 

6192-01 

5,3145*02 

3 

,8566-0* 

6. 

4422-02 

3,7095*03 

1 

, 8l68-o5 

2. 

4315-04 

r.iiiini 

1 

5 

•o9l"05 

,5486-06 

4. 

6. 

1697-02 

4396-02 

4,0439-03 

1 

,7395-05 

2. 

4116-04 

2, 9547*06 

1 

,7912-03 

•1. 

7356,02 

1,1155*03 

7 

,4920-06 

6. 

4419-02 

4,3580-03 

X 

; 7023-05 

2. 

4012-04 

3.1*42*06 

1 

,6o9l*o3 


9334-03 

1.5227*03 

9 

.7270-0* 

6 « 

4414-02 

4,4417-03 

1 

.7021-05 

2. 

4565-04 

3,2455*06 

X 

,8565-05 

3. 

7887-02 

2.0133*03 

X 

,2297*05 

6, 

4*22-02 

4.5218-03 

1 

,7192-05 

2. 

5936-05 

3.3041*0' 

2,5879*1,- 

1 

,*l73-05 

7. 

7647-02 

1 

,5253-05 

6. 

4422-02 

4.5523-03 

X 

.7329-05 

2, 

7692-04 

3.3265*06 

X 

,97t9- 0 5 

8. 

752n-02 

3.2569*03 

X 

,6652-05 

6. 

4423-02 

4,7067-03 

X 

.7372-05 

2. 

93 0 »-04 

3.4395*0* 

2 

, 0143-05 

7, 

016--02 

4.0420*03 

2 

.2361-05 

6. 

44 2 >-o2 

4.»2o4*o3 

X 

,7426-oS 

3. 

g36i-o4 

7999*02 

3.5226*06 

2 

»o544-p5 

4, 

4,«7io*C3 

2 

,*057-o5 

6 • 

4416-02 

4.»2l3-0J 

1 

.762!-r5 

3. 

l*2«-04 

3.5967*36 

2 

,1057'-c5 

3. 

4229-02 

6,0*35*03 

3 

.2220-05 

6 1 

4425-02 


R-Pfi 

x»a 


1.727i*Cfl 


1.7271*00 

’'3.3*05*90 


1.7271*00 

"3,3*o**00 


1,7271*00 
•3, 36o3*00 


1,7271*00 

•3,S6o2*00 


1.7271*00 

-3,5*01*00 


1.7271*00 

-3.3*00*00 


1,7271*00 

-3.55*8*00 


1.7271*00 

-3.3397*0(3 Ifc. 


1.7271*00 

-3.5394*00 


1.7271*00 

•3.3592*00 


1.727i*oo 
• 3. 559o*00 


1.727i*oo 

-3.3386*00 
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4, 

Or 
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CE 

ST 
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zeta 


k »P« 

X 


THfTiW 

TS 
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Gw 
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s 


04 

Hfc 

9H0£ 


TW 
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SUMO* 


AOO*A& 


otios 



52.0 

1 

,7?2«*n:> 

4.2033*02 

6, 1093*02 

1 

,7663-06 

1 

,9963-03 

1,0783-03 

7 

.4267-04 

1 

021C-02 

1* 

7247*oo 

•3.0340403 

0 

.oooo 

6.5859*n3 

4,4449*05 

2 

,9558-02 

3 

,4743*00 

8,9659*05 

6 

.6182-0* 


10*9-02 

•3. 

06i7-00 

4,9758-01 

4 

.2351-02 

1.4144*03 

1.9682-02 

1 

.9000*03 

0 

.0000 

5,0769*05 

S 

.4*12-03 


2*33-01 



33, o 

1 

,722»*oo 

4.3183*02 

8.3302-02 

1 

,7882-06 

2 

1 0*32-03 

i.0’72-03 

7 

.6242-0* 


0*23-02 

it 

7246*00 

•3.0096*00 

p 

• 0000 

8.3857*03 

4,4438*03 

2 

.9553-02 

3 

8p92*00 

9,39n*o5 

6 

,61*1-04 


3039-03 

•3. 

0l7t*00 

5.4235-01 

4 

,4489-02 

1.414**08 

1,9678-02 

1 

.9000*03 

0 

.0000 

9,4660*05 

9 

.1*36-03 


5*51-01 



34,0 

1 

.7226*00 

4,4558*02 

8,5955-02 

1 

7B8?"fl6 

2 

,1644-03 

1.119C-03 

7 

,7353-0* 


0*8--o2 

!• 

7266*00 

•2.9*29*00 

0 

• 0000 

8.5454*03 

4.4423*05 

2 

,9550-02 

4 

,2267-oq 

9,6799*05 

6 

,9694-0* 


6434-03 

"2. 

*7o3*00 

5.4874-01 

4 

.8892-02 

1,4143*08 

1.9872-02 

1 

.9oon*o3 

0 

• 0000 

5,9*12*05 

9 

,9877-03 

2 

•441*01 



93,- 

1 

,7228*pc 

4 ,8l66*02 

8,9o99*02 

1 

,786i-o6 

2 

26o<-o3 

1.1*41-03 

7 

,8224-04 


1273-02 

1* 

7246400 

•2.9138400 
8, 3778-01 

3 

.0000 

6,S*Sn*03 

4,44n7*05 

1,9666-02 

2 

.9547*02 

4 

.7411*00 

1,0*67*06 

7 

,0664-04 


0096*03 

•2* 

*212*00 

4 

.9027-02 

1.4142*04 

1 

,9000*03 

C 

• 0000 

6,4435*05 

l 

0*19-02 


320*-01 



38.0 

X 

.7224*00 

4,8o99*o2 

9,2792-02 

i 

,786l-o6 

2 

,36q3-q3 

1.1696-03 

7 

,8303-0* 


l5o2»02 

It 

7*66*00 

•2,8824*00 

D 

• 0000 
.1473-02 

8,5846*03 

4,4387*05 

2 

,9382-02 

f > 

,3673*00 

1.113**06 

7 

.11*1*04 

,16*4-02 


41*2-03 

*2. 

6696*00 

8.8924-01 

5 

1.4142*08 

1,9638-02 

1 

,9oon*o3 

0 

,0000 

7,0039*05 

1 


7164-01 



37.3 

1 

.7229*00 

5.0221*02 

9,68*9-02 

1 

,7880*06 

2 

i44o7-03 

1,1874-03 

7 

,78*o-o4 


1656-02 

1* 

7266*00 

•2.9lo9*00 

0 

• 0000 

6,3841*03 

4.4364*05 

2 

,99ii-o2 

6 

,068n*oo 

1.1*02*08 

7 

, 1313-0* 


9*23-03 

•2# 

6110*00 

7.4072-01 

5 

,1476*02 

1. 4140*0® 

1,9650-02 

1 

,9000*03 

0 

• OOOC 

7,5*82*05 

1 

,2578-02 


1217-01 



98,5 

1 

,7228*0u 

5.2351*02 

1.0139-01 
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.8299-02 

-2,13o5-o2 

7.3420-01 

5.0128*02 

9,2359-04 

1 

,4oi7-o3 

3,9642-02 

1,021*800 

2,6999+06 

9,0*28+06 

i 

,7*84-03 

••323-02 

-5.9498-02 

-7,4694-01 

5,2111*02 




EQUILIBRIUM TEST CASE 


FLAGS AND OPTIONS 


IDEAL » 0 

LANNR * 0 

INCOMP* 0 
J2D » 1 

INTD* » 1 


, . ms PERFECT GAS# *0 FOR HYDROS! N-OXYGEN EQUILIBRIUM) 

:sr.Si t r 0 i#M.., a L 5 . 

T.ssnSSi' ‘JSS ® '«! A#*. ««“"»• 


PROBLEM LIMITS AND INITIAL VALUES 

KS.: W8SU 585,: -S:SS»~ 


XMAX 


. •3,3448903*00 D*i * 6,3499999-09 


REPERENCE QUANTITIES ftnnnnnnAf ,, r H0 REFp 3.0000000-0 »MUR|F» 3.0000000-09 

BLREF * 3.0502500-01 UREP • 1.0000000*03 RHOREF* 3.0°° 


input normalization factors 

XN s 3.0502300-01 YN 


. 1.0000000*00 U6M • 1.0000000*00 PEN 


* 4,4595800*03 SMON * 1.1190000*00 


EDGE QUANTITIES 

UEDCE * 0.0000000 

CONSTANTS 

AFTRNS* 9.9990000-01 
PAMB * 0.0000000 

UE« * 0.0000000 


PE05E . 1,0000000*00 now ■ #. ««««-.# *«»«' 


PR] » 0*0000000 

$P 0 s l. 0000000+03 
RHOtK * 0. 0000000 


CONVERGENCE ANO EDGE CRITERIA 

CONVRG* 4.9999999-03 EPSLN1* 9,9999999-02 


SAMMA • 0.0000000 

SN3 • 4,0000001*00 


FMOLVT* 0,0000000 
X8TAR * 0,0000000 


EPSLN2* 9.9999999-02 EPSLN3* 9.9999999-02 


PLAW * 7,0000000*00 

AFVAI.u* 9.3999999-01 


COUNTERS 


MAX IT • 1 

IYEQ « 8 


NY I * 1*0 


NLPRNT* 20 


NSPRNT* 20 


IYPR * 1 


. REPRODUCIBILITY OP I HE ORIGINAL PAGE IS POOR 



I 

I TEPSIZE CONTROL T ARLES 

LOXLIMs a HKTAP* 

XL IH SKTAB 

-3.1407*00*00 1*1503000*00 

-5.0000000-^1 1.0300000400 

-5 ,0000000-02 1*0000000400 

0*0000000 1*1000000*00 

2.00000U0-C3 1*0100000*00 

1.9120000*00 1.1500000*00 

3,0100000*00 1*0000000*00 

2,2000000*01 

I 

I 

I 


| ?.:>r>ooom}-o? 

i.!>'Tonooo-o2 

il.i.-1000P0-03 

9.9999999-09 

I ?.ooooono-o 3 
9.9999999-03 
2.0000000-01 
I 2.0000000-01 





XTA3SK 

- 1.53) 5$0'"i,G0 

-3.14' 75CK ♦OS 

3,2 VH00--02 

V,4.44i941?-n; 

l,9l2D!)00*nc 
3,nu ooor+oo 
?,2o?ooor*oi 




wall tables 


i 

i 

i 


pnmm-' ■ w 


LTWT4B* 12 LW0T48» 2? 


TWTAB 

XT48TW 

SM0T4B 

XT43MQ 

1.4000000+03 

-3.5516000+00 

6.6000000-03 

-‘S,55U<1Qr’+0Q 

1 .4000000+03 

-1.9120000+00 

6.6000000-03 

-3.4140000+00 

1.3000000+03 

-1.3660000+00 

6,6000000*03 

-3,1420000+00 

1.3000000+03 

0.0000000 

6.6000000*03 

-2,8690000+00 

1.4000000+03 

0.0000000 

7.3000000-03 

-2,5950000+00 

1.4000000+03 

1.9120000+00 

4.6000000-03 

-2,3220000+00 

1.3000000+03 

1.9120000+00 

4.6000000-03 

-2,0 490000+00 

1.2000000+03 

3,0050000+00 

5.9000000-03 

-1,7760000+00 

1.3000000+03 

4.1000000+00 

4.6000000-03 

-1.503000O+00 

9.5000000+02 

6. 0300000+00 

6,6999999-03 

-1. ,2290000+00 

1 .0600000+03 

8.4700000+00 

6,9999999-03 

-'1,9699999-01 

1.5000000+03 

2.1390000+01 

1,3600000-02 

-i, 3300000-01 



1.4800000-02 

-4 ,9700000-01 



2.2600000-02 

-5,5999999-01 



2,7600000-02 

-4,2300000-01 



3.9700000-02 

-2,8700000-01 



3.6600000-02 

0.0000000 



7.2200000-03 

4,1000000-02 



5.1299999-03 

1,3700000-01 



9,5999999*03 

2,6000000-01 



4,6600000-03 

3,5500000-01 



2.2500000-03 

5,4600000-01 



2.3300000-03 

8,0599999-01 



1,5700000-03 

1,1340000+00 



7.9999999-04 

1,5850000+00 



0,0000000 

1,8300000+00 



0.0000000 

2,2000000+01 







1 


•• I ' M 


•- *!f 


IMERHO DRUM 


EACTANTS 
2. 0003 

-.0000 

-.onoo 

-.nooo 

-.0000 

100.0000 

2.0000 

-.0000 

-.0000 

-.0500 

-.0000 

100.0000 

0 INPT2 VALUE 

GIVEN FOR 

OF, EQRAT. 

FA, 09 FPCT 




FUEL OXIDANT MIXTURE 


'.CAL/G 


.nooooooo 


.00000000 .00000000 


.99206350-00 .00000000 


.49603175-00 


.00000000 -.12500000-00 -.62500000-01 


TQMS/G 

.99206350-00 

.00000000 

.49603175 00 


.00000000 

.62500000-01 

.31250000-01 


.00 G 29f< , 1 5 .3 p 

.00 G 2 9 p . J 5. 0 


-.000'; J 

-.COO'- 


REC1ES BEING CONSIDERED IN THIS SYSTEM 

J 9/65 H J 3/61 H2 J 3/61 H20 J 6/62 n w 3/66 OH 

J 9/65 02 


I 

I 

I 


I 

I 

I 

I 

I 

I 





377228-00 

.413012-00 

,688972-00 

.105301*01 

112484*01 

.115393*01 

.119067*01 

.123650*01 

170540*01 

.189601*01 

.215414*01 

.244290*01 

4n7247*01 

.445469*01 

,484676*01 

.529796*01 

652017*01 

.667421*01 

.683096+01 

.698330*01 

7 74440*01 

.789845*01 

.805364*01 

.821040*01 

902533*01 

,9t967l*ol 

.937168+01 

.949211*01 

947803*01 

.947484*01 

.947169*01 

.946864*01 

953824*01 

.954682*01 

.957279*01 

.959194*01 

985513*01 

.990898*01 

,998116*01 

.100*63*02 

103765*02 

.104498*02 

.105234+02 

.105992+02 

109837*02 

.110523*02 

.111430+02 

.111886*02 

114719*02 

.115183*02 

.115696*02 

.116347*02 

122319*02 

.123495*02 

. 124891*02 

.126117+02 

132217*02 

.132766*02 

.133511*02 

.134329*02 

137428*02 

.138043*02 

.138347+02 

.138604*02 

139634*02 

.139807*02 

.139941*02 

.140063*02 

355160*01 

-.300500*01 

-.245900+01 

-.187192*01 

149686*01 

-.141895*01 

-.134102*01 

-.126308*01 

873432-00 

-.795502-00 

-.717571-00 

-.639641-00 

249988-00 

-.167187-00 

-.643858-01 

-.357650-01 

293849-02 

.397025-02 

.503326-02 

.612121-02 

119067-01 

.131306-01 

,143765-01 

.156445-01 

223171-01 

.237i9l-ol 

.251440-01 

.266007-01 

338048-01 

.352355-01 

.366622-01 

.380876-01 

670349-01 

.833879-01 

.107911*00 

.143987-00 

477968-00 

.545380-00 

.613299-00 

.680990-00 

101981*01 

.109017*01 

.116094*01 

.123301+01 

160784*01 

.168542*01 

.177762*01 

.185629*01 

230*81*01 

.239746*01 

.249988*01 

.261859*01 

412420*01 

.451557+01 

.499618*01 

.552411*01 

909825*01 

.964145+01 

.102824*02 

.109974*02 

149125*0? 

. 15El9l*o2 

.163154*02 

.16802**02 

193873*02 

.199495+02 

.206970+02 

.210580*02 


.106759+01 

. 107420*-! 

.10*560*01 

.110255*01 

,129335+01 

.136362*11 

.145134*01 

.156239+01 

,274405+01 

. 305801+01 

. 336452*01 

.372300*01 

.599017*01 

.604593+01 

.620935+01 

.636179*01 

.713519+01 

.728701+11 

.743892*01 

.789130 ‘01 

.8366*4+01 

.852924*)! 

, -69i99+oi 

.885720+01 

,948964*01 

,948658*11 

.948396*01 

.946113*01 

,9+9366+01 

.951201*01 

.952968*01 

.951633*01 

,963354*01 

,968670*01 

.974074*01 

.980204*01 

,101097*02 

,101731*02 

,10?346*02 

.J ;3074*02 

,136770*02 

.107551*02 

.108357*02 

.1-9136*02 

,112507*02 

.112929*n2 

.113866+3? 

.114055+02 

,117801+02 

.118837*02 

,1l9976*c? 

.121187*02 

,127459+02 

.128701*02 

.13011**02 

.131544*02 

.13*9*3+02 

,135539*02 

.136172*02 

.136619*02 

.138836*02 

,140193*02 

-.180860+01 

.139042*02 

.139252+02 

.139437*02 

-.173067*01 

-.165274*01 

-.157481*01 

-.118515*01 

-.110722*01 

-.102929*01 

-.951363-00 

-.5*1710-00 

-.483780-00 

-.405849-00 

-.3279i9-oo 

-.3*32*1+02 

,000000 

.947994-03 

.192875-02 

,723291-02 

,836758-02 

.952488-02 

.107047-01 

,1*93*5-01 

.182466-01 

.195810-01 

.209378-01 

,2*0562-01 

,294914-01 

.309300-01 

.323697-01 

,403695-01 

,437i50-01 

.487323-01 

.558642-01 

.197421-00 

,269520-00 

.340154-00 

.4-9749-00 

,748172-00 

.815119-00 

,•*81601-00 

.950*98-00 

,130588+01 

,138020*01 

,145511*01 

.153175*01 

.194111+01 

.202597*-! 

.213574*01 

.220428*01 

,290300+01 

.313795*01 

.34*146*01 

.374948*01 

,614682+01 

,683468*01 

.765933*01 

.856835*01 

.116880+02 

.123959*02 

.131751*02 

.140095+02 

.172957*02 

.213908+02 

,177878*02 

.183079*02 

.188294+92 




EQUILIBRIUM TEST CASE 


STATION X S DS RW 



C 

-3. 55t6000+00 

0.0000000 

6.5499999-05 

1.7227900+00 


UE 

* 

EDGE AND WALL 
.3772285+1)3 

CONDITIONS 

ME 

s 

.7202565-01 

DELTA* 

3 

TE 

9 

.6647*80+04 

SHOE 

s 

.1804808-01 

theta 

3 

HE 

X 

-.7968727+07 

MUE 

s 

.2117189-05 

TAUW 

X 

PE 

* 

.44*6301+06 

TW 

s 

.1400000+04 

OW 

r 




MDOTW 

3 

.0000000 




:0. 

YBAR 

Y 

U/UE 

H/HE 

eO/HOE 

x 

0.0000000 

-o.ooooono 

0.0000000 

-8.2290472+00 

7.1139076-01 

2 

4.4802745-09 

9.2455547-03 

6.2901413-05 

-8.2284673+00 

7.1139459-01 

3 

1,C>78445-O0 

1.8491.109-02 

1,4599036-0* 

-8.2277012+00 

7 . Il39i93-Ul 

4 

1.00300*0-08 

2.7736604-02 

2.5313510-0* 

-8.2267133+00 

7 .11382*1-01 

5 

2.7077740-00 

3.6982219-02 

3. 8058534-c* 

-8.2254644+00 

7.1136720-01 

6 

3.9668914-08 

4.6227773-02 

5.5693703-0* 

-8.2239123+00 

7,1134512-01 

7 

5.4355053-08 

5.5473328-02 

7.6313347-04 

-8.2220110+00 

7,1131056-01 

a 

7.2114003-08 

6.4718082-02 

1.0124654-03 

-8.2197123+00 

7.112*151-01 

9 

9,33*80*9-08 

7.3964437-02 

1.3105/19-03 

-8.2169630+00 

7.1123999-01 

10 

1.10481/7-07 

8.3209991-02 

1.663**06-03 

-®. 2137101+00 

7.Hl93no-01 

11 

1.4796750-07 

9.2455546-02 

2.0774 0 95-o3 

-8.2098933+00 

7.1114383-01 

12 

1.8220159-07 

1.0170110-01 

2.5591608-03 

-8. 2054514+00 

7.1109257-01 

13 

2.2192553-07 

1. 11)94666-01 

3.1157533-03 

-8.2003190+00 

7.1103610-01 

8 1« 

2,6742599-07 

1.2019221-01 

3.7545630-03 

-8.1944297+00 

7.1090979-01 

15 

3.1933492-07 

1.2943777-01 

4.4833453-03 

-8. 1*77100*00 

7.1(108097-01 

16 

3.7822970-07 

1.3060332-01 

5.31020*5-03 

-8.1000802+00 

7.1379054-01 

17 

4.4471333-07 

1.4792900-01 

6.2430123-03 

-8.1714002+00 

7. 1067400-01 

1® 

5.1941452-07 

1.5717443-01 

7.2923898-03 

-8.101*102+00 

7.1054570-01 

19 

t. 0290801-07 

1.6641998-01 

8.4657309-03 

-8.1509916+00 

7.10*0912-01 

i 20 

6.901J468-O7 

1.7566554-01 

9.7731952-03 

-8.1389364+00 

7,10264*8-01 

21 

7.9950183-07 

1.8491109-01 

1.1224713-02 

-8. 1255531+00 

7,1010909-01 

22 

9.1388348-07 

1.9415605-01 

1.2830590-02 

-8.1107403+00 

7.0993574-pi 

23 

1. j|4OO2n7-O0 

2.0340221-01 

1.4*01511-02 

-0.0944170+00 

7,09730*8-01 

24 

1.1787020-06 

2.1264776-01 

1.6548545-02 

-8.0764653+00 

7.0951454-01 

25 

1.3307434-06 

2.2189331-01 

1.8683151-02 

-8.0567832+00 

7.092*286-01 

26 

1.4969893-06 

2.3113087-01 

2.1017102-02 

-8.0352625+00 

7.0*98938-01 

27 

1.6783126-06 

2.403*442-01 

2.3562895-02 

-8.0117897+00 

7.0*69984-01 

2» 

1.8756154-06 

2.4962997-01 

2. 6332954-02 

-7.9062401+00 

7.0*39437-01 

29 

2. 0890294-06 

2.3087553-01 

2.9340439-02 

-7.9585173+00 

7.0*0*877-01 

30 

2.3219164-06 

2.6812109-01 

3.2598858-02 

-7.9284725+00 

7.0771670-01 

31 

2.5720697-06 

2.7736664-01 

3.6122149-02 

-7,8959852+00 

7.0733199-01 

32 

2.8437130-06 

2.8661220-01 

3.9924690-02 

-7.8*39226+00 

7.0691030-01 

33 

3.1355050-06 

2.9585775-01 

*.*0213*3-02 

-7.8231479+00 

7.06*51*0-01 

34 

3.4493340-0* 

3.0510330-01 

*.8427387-02 

-7.7825202+00 

7.0596x78-01 

35 

3.7863249-06 

3.1434886-01 

5.3158610-02 

-7.738*933+00 

7.05**435-01 

36 

4.1476371-06 

3.2359441-01 

5.8?" • 'i6-0? 

-7.092J j75*oo 

7.0489959-01 

37 

4.5344659-06 

3.3203997-01 

* 1-02 

• 7 . 6*20377+00 

7,0*32492-01 

3® 

4,9480443-00 

3.4?0®5S3”0t 

6 j-02 

-7.5*84943+00 

7.0371043-01 

39 

5.3896432-0* 

3.5133108-01 

7 , . .>(..1618-02 

-7.5313225+00 

7.0307020-01 

40 

5.8005736-06 

3.6057663-01 

8.2200308-*? 

-7.4703526+00 

7.0230367-01 

41 

6.3621877-06 

3.6982219-01 

8.9322788-02 

7.4054094+00 

Oi ->5092-01 

42 

6.8958804-06 

3.7906774-01 

9.6815637-02 

-7.3363120+00 

7.0009109-01 

43 

? .463o9o3-06 

3.8831329-01 

A. £477906-01 

-7.2*20739+00 

7 . 0008201-01 

44 

8.0653027-06 

3.9755805-01 

1.1323390-01 

-7.1849028+00 

6.9922957-01 

45 

8.7040505-0* 

4. 0680441-01 

1 . 2220168-01 

-7.1021999+00 

6.9*33048-01 



d- „ , , 


THE TAW 7ETA 

000000" 1.36A0062-03 3. 

P“OP I LE PARAMETERS 
.1353381-03 CP * 

,31*9647-04 ST ■ 

, 18*1192+01 
-.52*2559+06 


RCV EPS 

1.8?2189.<+0O 0.0000000 

1.8^21943+03 7.6749931-16 

1.6')22061 + 00 2.6735604-13 

1. «C222l4+oo 1.2oS3#n5-i2 

1.802240,7*00 4.5o7i6r>l-l2 

1.8022646+00 1.4716925-11 

1.6022940+00 4.3304865-n 

1.8023295+00 1.17245*4-10 

1.8023720+00 2.9637622-10 

1.8024223+00 7.0*96713-10 

1 . 8*1 24812+00 1.6n4i8l8-09 

1.6025499*00 3.4842619-09 

1.8026292+00 7.2796903-09 

1.8027202+00 1.4489673-0* 

1.80282«0+00 2.87253P7-06 

1.8029418+00 5.459?767-08 

1.0030747+00 1.0109171-07 

1.80322*1+00 1.82/6402-37 

1.6033913+00 3.2328513-07 

1.8035775+00 5.6oi.7o00-Q7 

1.8037843+00 9.5215766-07 

1.8040131+00 1.5894650-06 

1.0042654+on 2.6006263-06 

1 , 0045427*00 4.2135566-06 

1.8048468+00 6.7f)5l78O-O0 

1.0051793+00 1.0521834-05 

1.8055419+00 1.6293465-05 

1.80593*5+00 2.4913298-05 

1 . 0 n 6365o + OO 3 . 7 S3l2'’8-0S 

1 ,0068291+00 5, 6i 7 4945-05 

1. 0073310+00 8.2909239-05 

1.6370727+00 l.?l05014-0* 

1.8084563+00 1.7494795-04 

1.0090840+00 2.50*2600-04 

1.80975*0+00 3.551«l66-04 

1.8104806+03 4.9925143-04 

1. 8112542^00 0 . 955*467-04 

1-8120814+00 9.60*7245-0* 

1.0129040+00 1.3147786-03 

1 .0139c -5+00 1.7847720-03 

1.0149097+00 2.4*39007-03 

1.0159771+00 3.2145748-03 

1.8171115+00 4.2605692-03 

1.8103159+00 5.6287342-03 

1.0195934+00 7.3691250-03 


ZETAP 

02 


.1372945-02 

.1011224-02 


T 

1 * 00.0 
1399,9 
140C.0 
1 * 00,1 
1*00.3 
1 * 00,6 
1*00,9 
1*01.4 
1*01,9 
l*n2. 5 
1*03. 1 
14*3.7 
1*0*. * 
1*05.2 
1*06.2 
14*7.5 
1 * 0 ’. 0 
1 * 10.6 
1*12.3 
1*1*. 1 
1 * 10.1 
1*1® , 3 

1420.8 
1*23,7 
1*27,1 
1*30.7 
1*3*. 5 
1*38.5 
1*42.7 
1**7, 3 

1452.6 
1*56.5 

1 * 6 - a 1 

1*72.2 

1*79,5 
1*87,2 
1*95.3 
l5o* • 1 

1513.9 
152*.9 

1537.1 

1550.1 

1563.6 
1577.4 
1592.0 




M0 + 

YBAP 

Y 

U/UE 

H/ME 

PO/ftOE 

PC/ 

EPP 

T 

46 

9.3809158-0* 

4. 16Q4996-Q1 

1.3170463-01 

-7.0145602400 
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3.7786353*00 

1.2218203+OO 

1.7232295+00 

8.2032107-07 

1 4?7 , 9 


1« 

20 

21 

22 

23 

24 

25 
2» 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 
<0 
41 
4? 

43 

44 

45 


6.3322426-07 
7.3102067-07 
8.3959206-07 
9.597o9?»-o? 
1. 0921/16-06 
1,2378069-06 
1. 39747*3-06 
1,5720544-06 
1.7624699-06 
1 .9696*63-04 
2. 1946*19-06 
2.4383467-04 
2 . 7 018* 38-06 
2.9863089-06 
3.2927319-06 
3.6222974-06 
3.9761864-06 
4.3556163-06 
4.7618423-06 
5.1961591-06 
5.6599316-06 
6.1544442-06 
6.6812134-C6 
7.2416676-C4 
7.6373197-C6 
8 . 4697294- C 6 
9.1405065-C6 


1.6641998-01 
1.7566554-01 
1.8491109-01 
1.9415665-01 
2.0340221-01 
2. 1264776-01 
2.2189331-01 
2.3113887-01 
2.4038442-01 
2.4962997-01 

2.5.187553-01 

2.6812109-01 

2.7736664-01 

2.8661220-01 

2.9585775-01 

3,0510330*01 

3.1434866-01 

3.2359441-01 

3.3283997-01 

3. 4208553*01 
3.5133108-01 
3.6057663-01 
3.6982219-01 
3.7906774-01 
3.8831329-01 
3.9755885*01 

4. Q680441-01 


1.5886240-02 
1.8330316-02 
2.1040820-02 
2,4036126-02 
2.7334949-02 
3.0956320-02 
3.4919616-02 
3.9244746-02 
4.3952326-02 
4.9063550-02 
5.4599838-02 
6.0582451-02 
6. 7fj32ll4-02 
7.3968757-02 
8.1411353-02 
8.9377986-02 
9.7885720-02 
1 . 0694989-01 
1.1658317-01 
1.2679454-01 
1.3758030-01 
1.4096341-01 
1.6091334-01 
1.7342634-01 
1.8648513-01 
2.0006603-01 
2 . 1414316-01 


3.7625051*00 
3,7445246*00 
3.7245537*00 
3, 7024497*00 
3.6780681*00 
3.6512645*00 
3.6218952*00 
3,5898157*00 
3.5548780*00 
3.5169329*00 
3,4758412*00 
3.4314829*00 
3.3837652*00 
3.3326244*00 
3.2700211*00 
3,2199239*00 
3.1582914*00 
3.0930678*00 
3.0241971*00 
2.9516286*00 
2.8753177+00 
2.7952202+00 
2.7112051*00 
2.6234469+00 
2.5316387*00 
2.4358319*00 
2.3360773*00 


1.2186611+00 

1.2156986*00 

1.2122647*00 

1.2065339*00 

1.2043611*00 

1.1997099*0? 

1 . 1945818*00 

1.1090372*00 

1.1830962+00 

1.1767589*00 

1.1700238*0? 
l . 1628686*00 
l.l*535lC*on 
1.1474095*00 
1.1790643*0? 
1 . 1 "'02926 + 00 
1.1210160*00 
1.11123*9+00 
1. l? 10206 +o? 
1.Q90484V+00 
1.0797391*00 
1.0*88566*00 
IOH/8354 + oo 
1. 0465521*00 
1.0348355*00 
1.0226820*00 
1> 0102297*00 


1.7233135*00 
1.7234154*00 
1.7235496*00 
1.7237195*00 
1. 7239325*00 
1.7241970*00 
1.7245217*00 
1.7249157*00 
1.7253885*00 
1.7259506*00 
1. /266140+C0 
1.7273886*00 
1.7282624*00 
1.7292973+00 
1.7304201*00 
1.7316616*00 
1.7329770*00 
1.7343477*00 
1. 7357440*00 
1.7371339*00 
1.7384628+00 
1.7397520*00 
1.7400946*00 
1.7418474*00 
1.7425278+00 
1.7428483*00 
1.7427248*00 


1.4648926-06 
2.5020044-06 
4.21/6227-06 
6.9759472-06 
1.1332439-05 
1.8100821-05 
2. 8457943-05 
4.4n8l067-05 
6.7323246-05 
1 .01 4 3 i 86-04 
1.5?6i477-04 
2.2136313-04 
3.2105266-04 
4.5945219-04 
6.5?30593-04 
9.1^16658-04 
1 . 2600585*03 
1.7?592?3-03 
2.3391845-03 
3.13/2798-03 
4.16417*7-03 
5.4708190-03 
7.1150722-03 
9.1617189-03 
1.1682339-02 
1.4754886-02 
1.8462940-02 


1432.4 

143/, 0 

1442.1 
144/. 7 

1454.1 

1461.5 

1469.8 

1478.9 

1488.4 

1498.5 
lS r 9.3 

1520.9 

1533.8 

1547.9 

1563.4 

1580.0 

1597.5 

1615.9 

1635.4 

1656.1 

1678.2 

1702.0 

17?7.3 

1754.0 
1 / 82,1 

1511.6 

1842.4 




\ 






am r; 



f -'ll 


I 


*!0 . YR4Q Y U/UP H/ME 

46 9,«M3l?6-06 4.1604996-01 2.2868()36-0l 2.2323237*00 

4? 1.0803864-03 4.2529552-01 2.4363895-01 2.12541*0*00 

48 1.1399935-05 4.3454107-01 2.5897437-01 2.0150769*00 

49 1.2241357-05 4.43786*3-01 2.7463853-01 l,»m8790*0O 

50 1.3l3ou?4-05 4.5303218-01 2.9o582l8-01 1.7862093*00 

51 1.4067897-05 4.6227773*01 3,0675592-01 1.6*84372*00 

52 1.5057002-05 4.7i52329*01 3.2311158-01 1.5489189*00 

53 1.6099437-05 4.81)76884-01 3.3960320-01 1.4880212*00 

54 1,7197376-05 4.9Q01440-01 3.5618665-01 1.3061310*00 

55 1.8353070-05 4.9925996-01 3.7281896-01 1.1836607*00 

56 1,9568852-05 5,0850550-01 3.8945835-01 1.06106*5*00 

57 2.0847142-05 5.1775106*01 4,0606579-01 9.3876727-01 

58 2.2190449-05 5.26996*1*01 6,22*0641-01 5.1731294-Q1 

59 2.3*01381-05 5.3*24216-01 4,3904895*01 8.9717290-01 

60 2.5082645-05 5.4543772-01 4.5536577-01 5.7882588-01 

61 2.6*37051-05 5.5473328*01 4.7i53265*01 6.6275734-01 

62 2.8267524-05 5.6397883*01 4.875273Q-01 3.6944451-01 

63 2.7977106-05 5.7322439*01 5.0332816-01 2.3935197-01 

64 3.17689*1-05 5,8246994-01 5.1891527*01 1.3894292*01 

65 3.36463*8-05 5.9i7i550-01 5,34271*4-01 3.0*74889-02 

66 3.5612819-05 6.0096105-01 5.4938032-01 -6, 7035187-02 

67 3.7*71836-05 6.1020640-01 5.6422614-01 -1.5978400-01 

*8 3.9827178-05 6.1545216*01 5,7B800«9-01 -2.4717134-01 

69 4,2082738-05 6.2869772-01 5.9309557-01 -3.2887003-01 

70 4.4442589-05 6.3794327-01 6.0708608-01 -4,0466870-01 

71 4.6910989-05 6.4718883-01 6.2077820-01 -4.74098*2-01 

72 4.9492378-05 6.5643436-01 6.3419490-01 -5.3*65640*01 

73 5.2191412-05 6.6567993*01 6.4720737-01 -5.9383404-01 

74 5,5012962-05 6.7492549-01 6.6002532-01 -6.4800722-01 

75 5.7962124-05 6.84l7io4-Ql 6,7239302-01 -6,8549397-01 

76 6,1044243-05 6.93+1660”01 6,8451998-01 -7.2036240-01 

77 6,4264923-05 7.02**216-01 6.9*26652-01 -7.4903986-01 

78 6,7630034-05 7.u9o77l-Ql 7.0749463-01 -7,7053071-01 

79 7.1145748-05 7.2H5326-01 7.18781*1-01 -7.6521458-01 

60 7.4818547-05 7.3 0 398R2-Gl 7,29530*3-01 -7.93240*0-01 

61 7.8655230-05 7.3964437-Q1 7.3994223-01 -7.9*79699-01 

82 8,2662956-05 7.4888992-01 7.5001698-01 -7.9014765-01 

83 8.6849257-05 7.5813546*01 7.5975*07-01 -7,7959773-01 

84 9.1222051-05 7.6738103*01 7,691*0*7-01 -7, 6349515*01 

85 9,5789*86-05 7.7662659-01 7.7823345*01 -7 . 4222882*01 

86 1.005*096-04 7.8587215-01 7.8*97**7-01 -7.1*221*6-01 

87 1.0554513-04 7. 9511770*01 7.9539427-01 -8,8592425-01 

88 1.1075198-04 8.0436385*01 8,D349o74-oi -6,51*0848-01 

89 1.1619185-04 8.1360881-01 8.li27i52-01 -6.143*075-01 

90 1,2187562-04 8-2285436-01 8.1*74134-01 -5.7408150-01 

91 1.2781481-04 8.3209991-01 8.2590084-01 -5.3149536-01 

9? 1.3402161-04 8.4134547-01 8.3275800*01 -4.870*757-01 

93 1,4050890-04 8.5o59io3-01 8.39350*8-01 -4,4115758-01 

94 1.4729030-04 8.5983658*01 8.4572928-01 -3.93*79*9-01 

95 1,543*026-04 8.690*214-01 8.5l93752*01 -3.4539931-01 

96 1.6179407-04 8.7832769-Q1 8.5802422*01 -2.9577225-01 

97 1.6954798-04 6.8757325-01 8.6403536-01 -2.4501944-01 

98 1 , 7765917-04 8.96818*0*01 8.6999469-01 -1.9327247-01 

99 1,8614590-04 9.0606435*01 8.7591393*01 -1 . 4071698-01 

100 1.9502760-04 9.i53o99i-Ql 8.6180743*01 -8.74924*2-02 

101 2.0432485-04 9.2455547-01 8.8768864-01 -3.3715957-02 

102 2.1405958-04 9.3380102-01 8.9357140-01 2.0*32455-02 






/1'TtVF • 


•» o/Rnt 

9,9 768 7*6-0 1 
9.8526430-01 
9.7309681-01 
9.6118568-01 
9.4934lp9-oi 
9.3739340-01 
9.2534U2-01 
9.1330088-01 
9,0144923-01 
8.8996577-01 
8.7896746-01 
8.6*45188-01 
8.5*27500-01 
8.4834432-01 
8.30*5793-01 
8.2926805-01 
8.20253*7-01 
8.11*9380-01 
8.03640*1-01 
7.96092*1-01 
7.8899257-01 
7.8231535-01 
7.7606100-01 
7.7023640-01 
7.6465195-01 
7.5991805-01 
7.5544158-01 
7.5142238-01 
7.4786101-01 
7,447*074-01 
7,4218221-01 
7.4004138-01 
7.3»32l99-oi 
7.3698768-01 
7.3601449-01 
7.3540300-01 
7.3519564-01 
7.35454*3-01 
7.361*215-01 
7.3732598-01 
7.3880887-01 
7.4055316-01 
7.4250757-01 
7,4467341-01 
7.4710644-01 
7.4982747-01 
7.52*3659-01 
7.5*13011-01 
7.5970256-01 
7.6354838-01 
7.6766394-01 
7.7204937-01 
7,7666*28-01 
7.6137945-01 
7.8618561-01 
7.9121298-01 
7,9465555-01 


°C.' 

1.742n*l5*0" 
1 . 7 + 0+5 32+flP 
1.7389836*00 
1.7364112+03 
1 . 7330396*0 ; 
1.7287559*00 

1.7234711*00 

1.7171285400 
1.7097010400 
1 . 7011861400 
1.6915953*00 
1 . 6809240*00 
1.6691234*00 

1.6561377*00 

1.6419623*00 

1.6246354*00 

1.6102391*00 

1.5989081*00 

1.5747903*00 

1.5560651*00 

1.5368701*00 

1.5174401*00 

1.4980683*00 

1.4790103*00 

1.4607042*00 

1,4435743*00 

1.4275449*00 

1.4135981*00 

1.4017204*00 

1.3925034*00 

1.3871343*00 

1.3857392+00 

1.3893057*00 

1.3980583*00 

1.4123253*00 

1.4324110*00 

1.45885i7*oo 

1.4925910*00 

1.5346244*00 

1.5857258+00 

1.6462744+00 

1.7161993*00 

1.7950850*00 

1.8824698*00 

1.9783373*00 

2.0826083*00 

2.1951558*00 

2.3154423*00 

2.4425260+00 

2.5754735+00 

2.7134515*00 

2.8556096*00 

3.0008995*00 

3.1*73441*00 

3.2934590*00 

3.4389669*00 

3.5651406+00 


EP5 

2.2895133-02 
2.81**1*7-02 
3 . '• 30+593*02 
+ • 1+ ^ 1216-02 
«.97J9i+o- 02 
5.9?0+575-02 
6.996533+-02 
8.2121998-02 
9.5779o9i-02 
I.ll0*+e6-01 
1.2802313-01 
1. + 682122-01 
1.675+366-01 
1.9n29557-01 
2.1517986-01 
2. +229680-01 
2.7i74992-ci 
3.0363951-01 
3. 3805372-01 
3.7509523-01 
+.1+87887-01 
+.57++59+-01 
5.0286787-01 
5.51+8+75-01 
6.0271052-01 
6.5723792-01 
7.1+89517-01 
7.75+3U3-01 
8.3921394-01 
9.0612027-01 
9.7604152-01 
1.0+90093*00 
1.1253615*00 
1 • 2r+5479*00 
1.287m i*oo 
1.3725587*00 
1.4608958+00 
1.5519167*00 
1.6454977*00 
1.7415890*00 
1.8402351+00 
1.9414613*00 
2.0+50400*00 
2.1507263*00 
2.2595195*00 
2.3742*45*00 
2.4941564*00 
2.5914577*00 
2.6655997*00 
2.7146278*00 
2.73+5622*00 
2.7?50537*00 
2.7126023*00 
2.6963460*00 
2.6752021*00 

2.6478326*00 

2.6126199*00 


Ml'.. ; 




wo. 

Yf?4R 

Y 

u/ue 

H/*E 

-70/WOfc 

ROV 

EP5 

T 

*03 

2. 242551 0-'IM 

9.430*658-01 

8.9947027-01 

7.5212088-02 

8 . 0 •’.70734-0 1 

3. 73464?9*c>-» 

2.5*76657+00 

5387,6 

104 

2.3493623*04 

9.5229213-01 

9.0540025-01 

1.3031860-01 

8. 0949657-qi 

3. *90*6*2*03 

2.5lJ«322*00 

5464,2 

105 

2,4612941-04 

9.6153768-01 

9.1137710-01 

1.8588514-01 

8 . 1 7Q1998-C t 

4 . 0365731*01 

2.4398495*00 

5542,7 

10ft 

2.57862*1-04 

9.7q7832*-01 

9.1741793-01 

2.4t96iiB-01 

8.2517773-01 

4 . 2331 1*3*00 

2.3525132*00 

56?3 , 5 

107 

2.7016647-04 

9. 8002879-01 

9.2354020-01 

2.9*69869-01 

8 . 3404619-ut 

4. *?297?9*() r 

2.2469867+00 

57-9,1 

10» 

2.8307245-04 

9.8927435-01 

9.2976067-01 

3.5613858-01 

8.4362824-01 

4.6275327 + 0'' 

2.1221875+00 

5799,4 

109 

2.9661499-04 

9,9851990-01 

9.3609538-01 

*.1**0456-01 

8.5381245-01 

4 , *4675*2*00 

1.9782270*00 

5892,1 

110 

3.1083062-04 

1.0077655400 

9.4255835-01 

4.7358663-01 

8.6443029-01 

5.c79*l?6 + 0 ’l 

1 .8167640*00 

5983,5 

111 

3.2575845-0* 

1 , 01 70110*00 

9 , 4915829-01 

5,3372607-01 

8.7531318-01 

5, 325*5i 4+oo 

1.6411839*00 

60*9,9 

112 

3.4144031-04 

1.0262566*00 

9.5589253-01 

5.9475549-01 

8.5634977-01 

5 . 5*4 < 3+00 

1.4564718*00 

*148,8 

113 

3.5792094-04 

1.0355021*00 

9.6273693-01 

6.56*0400-01 

8.9754291-ni 

5.8582152*00 

1.2*87504*00 

6220,4 

114 

3.7524824-04 

1.0447477*00 

9.6962906-01 

7 . 1*05550-01 

9. 09*0591-11 

6.156*616*00 

1.0**5525*00 

*2*8,1 

115 

3.9347363-04 

1.0539932*00 

9.7644334-01 

7.7857286-01 

9. 2325225-01 

6.508*732*00 

9.0999068-01 

*355,7 

lift 

4.1265217-04 

1.0632388*00 

9.8296540-01 

8.3615263-01 

9.3911271-01 

6.9246454*00 

7.5002475-01 

6*23.3 

117 

4.3284295-0* 

1.0724843*00 

9.8867143-01 

8.8*27468-01 

9.5578737-01 

7.3*885*4*00 

6.0799435-01 

64*8,1 

11« 

4.5410950-04 

1.0817299*00 

9.9373836-01 

9.3191971-01 

9.71*6113-01 

7.8*18436*00 

4 . 8549009-01 

6545.8 

119 

4,7652004-0* 

1.0909754*00 

9.9716400*01 

9.6*38735-01 

9.8431882-01 

8.2*53011*00 

3.8252952-01 

6592.3 

120 

5.001*783-0* 

1.1002210*00 

9.9904405-01 

9.8474374-01 

9.9316058-01 

8,6014114*03 

2.9793011-01 

*6?4,7 

121 

5.2507190-0* 

1.1(194666*00 

9.9977495-01 

9,9489826-01 

9.9301718-01 

8.7882631+00 

2.2974926-01 

6643,1 

122 

5.5l37??o*0* 

1.1187121*00 

9.9996568-01 

9.9878411-01 

1,0001654*00 

8.8711301*00 

1.6409747-01 

6650.5 

123 

5.7915526-0* 

1.1279577*00 

1.0000000*00 

1.0000000*00 

1.0000000*00 

8 ,8567821+uD 

1.2406790-01 

6647.9 
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EQUILIBRIUM TEST CASE 


WO. 

YRAR 

0/F 

C( H ) 

C(W2) 

C<W20) 

-CO) 

CCOH) 

C ( 02 ) 

MIJ 

PR 

1 

Q.C000300 

8.1699194-01 

.000000 

.493706 

.506294 

•onoooo 

. oonooo 

.000000 

4.36348-07 

.55793 

2 

4.7049335-09 

8.1705880-01 

.000000 

.000000 

.000000 

.ocoooo 

.oonooo 

.uonr. :> 

4.36314-07 

.55793 

3 

1 ."919464-0. 

8.1714701-01 

.000000 

.000000 

.oooono 

.oooooo 

. 000000 

• uooooo 

4.36335-07 

.55793 

4 

1 .A934139-0" 

8 . 1726059-01 

.000000 

.000000 

i 0 0 C 0 0 

.oaooon 

.oonooo 

.OOOOC 0 

4.36411-07 

.55793 

5 

2 . 9Q65613-04 

8.1740384-01 

.000000 

.000000 

.oonooo 

.OOOOOO 

.oonooo 

.oooooo 

4.36540-07 

.55793 

6 

4.165»073-0« 

8.1756154-01 

.000000 

.000000 

.oooono 

•OOOOOO 

•oooooo 

.oooooo 

4.36725-07 

.55793 

7 

5.7081264-06 

8.1779895-01 

.000000 

.000000 

.000000 

.oooooo 

.OOOOOO 

.oooooo 

4,36963-07 

.55794 

A 

7.5730930-08 

8.1806159-01 

.000000 

.000000 

.000000 

•OOOOOO 

. oooooo 

.000000 

4.37256-07 

.55794 

9 

9,8028854-08 

8.1837557-01 

.000000 

.493234 

, 5067*6 

.oooooo 

. oonooo 

.001)000 

4.37603-07 

.55794 

10 

1.2442292-07 
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4.4463+05 

0.000° 

1,8412+00 


0. J000 

6.5500-05 

-7.9687+06 

1.88 *■’-02 

1.4000*03 

6.6750+00 

— 

1 .0 

1.722"+00 

3.7723+02 

7.2626-02 

2.1172-06 

1.2800-03 


-3.5515+Ov 

o.ooor 

6.6479+03 

4.4463+05 

7.3854-03 

1.7246+00 


6.5500-05 

6 . 5500*05 

-7.9687+06 

1.8848-02 

1.4000*03 

1.6116+01 

Sr 

2 • ' 

1.7228+00 

3.7724+02 

7.2627-02 

2.1172-06 

1.5716-03 


-3.5515+00 

o . noon 

6.6479+03 

4.4463+05 

7.3854-03 

2.04JT+00 

ir 

1.4082*04 

7.5325-05 

-7.9687+06 

1.8848-02 

1.4000*03 

8.3434+00 

3.1 

1.7228+00 

3.7724+02 

7.2629-02 

2.1172-06 

1.6018-03 

i 

-3.5514+OC 

C.OPOO 

6.6479+03 

4.4463+05 

7.3854-03 

2.1483+00 


2.2745-04 

8.6624-05 

-7.9687+06 

1.8848-02 

1.4000*03 

-3.8329+00 


4.0 

1.7228+00 

3.7725+02 

7.2630-02 

2.1172-06 

1.8751-03 


-3.5513+00 

0.0000 

6.6479+03 

4.4463+05 

7.3854-03 

2,5149+00 


3.2707-04 

9.9617-05 

-7.9687+06 

1.8848-02 

1.4000*03 

-6.8475+00 


5.C 

1.7228+00 

3.7726+02 

7.2631-02 

2.1172-06 

2.0285-03 


-3.5512+00 

0.0000 

6.6479+03 

4.4463*05 

7.3855-03 

2.7207+00 


4.4163-04 

1.1456-04 

-7.9687+06 

1.8848-02 

1.4000*03 

-3.1961+00 


6.0 

1.7228+00 

3.7727+02 

7.2633-02 

2.1172-06 

2.3113-03 

I s 

-3.5510+OC 

0.0000 

6.6479+03 

4.4463*05 

7.3854-03 

3.1002*00 


5.7337-04 

1.3174-04 

-7.9687+06 

1.8848-02 

1.4000*03 

2.5542+00 


7.0 

1.7228+00 

3.7728+02 

7.2635-02 

2.1172-06 

2.4-Q8-03 


-3.5509+00 

0.0000 

6.6479+03 

4.4463*05 

7.3854-03 

3.3006+00 


7.2488*04 

1.5151-04 

-7.9687+06 

1.8848-02 

1.4000*03 

6.2525+00 


8*0 

1.7228+00 

3.7729+0? 

7.2637-02 

2. 1172*06 

2.6922-03 


-3.5507+00 

0.0000 

6.6479+03 

4.4463*05 

7.3854-03 

3.6116+00 


8.9911-04 

1.7423-04 

-7.9687+06 

1.8848-02 

1.4000*03 

6.2979+00 


9.C 

1.7228+00 

3.7730*02 

7.2640-02 

7.1172-06 

2.7992-03 


-3.5505+00 

0.0000 

6.6479+03 

4.8463*05 

7.3854-03 

3.7553+00 


1.0995-03 

2.0037-04 

-7.9668+06 

1.8848-02 

1.4000*03 

4.5766-i 00 


to.c 

1.7228+00 

3.7732+02 

7.2642-02 

2.1172-06 

2.9689-03 


-3.5*03+00 

0.0000 

6.6479+03 

4.4463*05 

7.3854-03 

3.9433+00 

f 

1.3299-03 

2.3042-04 

-7.9688+06 

1.8848-02 

1.4000*03 

2.7471+00 


U'O 

1.7228+00 

3.7733+02 

7.2646-02 

2.1172-06 

3.0278-03 


-3.5500+00 

0.0000 

6.6479+03 

4.4463*05 

7.3854-03 

4.0627+00 

4 

1.5949-03 

2.6498-04 

-7.9688+06 

1.8848-02 

1 . 400°*03 

2.3555+00 

If 

ir.o 

1.7228+00 

3.7735*02 

7.2650-02 

2.1172-06 

3.1283-03 


-3.5497+00 

o.onoo 

6.647V+03 

4.4463*05 

7.3854-03 

4.1980+00 

i 

1.8998-03 

3.0473-04 

-7.9688+06 

1.8648-02 

1.4000*03 

3, 0371+00 
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ST 
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RrfPR 
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S* IMQvt 

TmETA 

Anr'iAS 

ZpTaP - 

O.iEDS 

XPR 

1.0112-03 

1.3534-04 

1 . 366"- f j3 

5.2536*01 

5.2626+05 

0.0000 

3 , 1 8 'i 6 - 0 5 
L' ,00C n 

3.1023-02 

-1.0«33+°0 

1.7543-03 

1.35««-04 

1.3725-03 

5.2536-01 

7.3415+05 
6. /912+00 

3,2039-05 
2, 7740-07 

9.5934-02 

6.5463-02 

-1.0833+00 

2.1211-03 

t.35i7-04 

1.3784-03 

5.2536-01 

7.2913+0:: 

6.2295+01 

3, 214^-05 
5,9726-0/ 

8.C7i?-02 

6,5431-72 

-1.0833+30 

2.7537-03 

1.3277-04 

1 .3795-03 

5.2536-01 

8.5447+05 

1.3137+02 

3. 21'6-OS 
9,6465-07 

3.1187-03 

6.5454-02 

-1.0833+00 

.5614-03 

1.3154-04 

1 . 37 18-03 

5.2536**1 

9.3642+05 

3.2096-05 

-6.0702-02 

-1.0832+00 

2.2121*02 

1.3871-06 

6.5443-02 


4,7324-03 

1 . 3044-04 

1 ,36?4-03 

5.2536-01 

1.0734*06 

3.3715*02 

3.2095-05 
1 , 87 J • -06 

-8.1953-02 

6.5459-02 

-1.0832+00 

6.4794-03 

1.2945-04 

1.3557-03 

5.2536-01 

1.1756*06 

4.8635+02 

3.2177-05 

2.43l«-C6 

-5.0943-02 

6.5461-02 

-1.0832+00 

8.9632-03 

1,2762-04 

1.3549-03 

5.2536-01 

1.2929+06 
6. 7464+02 

3.2337-05 

3,0743-06 

-5.4184-03 

6.5455-02 

-1.0831*00 

1.325H-02 

1,2964-04 

1.3594-03 

5.2536-01 

1.3590*06 

9,0735*02 

3,25j' -05 
3.8133-06 

2.6035-02 

6.5444-02 

-1.0831*00 

2.0797-02 

1.2967-04 

1.3659-03 

5.2536-01 

1.4344+06 

1.1892+03 

3,2674-05 

4,6630-06 

3.2376-02 

6,5464-02 

-1.0830*00 

4.4045-02 

1,2965-04 

1.37 0 5-03 

5. J536-01 

1, 4560+06 
1.5244+03 

3,28i3-05 

5,6403-06 

2.0178-02 

6.5462-02 

-1.0829+00 

3.4442-01 

1.2962-04 

1.3714-03 

5.2536-01 

1.4953*06 

1.9184+03 

3.2959-05 

6.7641-06 

4.79 0 9-03 

6.5457-02 

-1.0828+00 

5. 7147-02 

1.2975-04 

1.37U-03 

5.2530-01 

1.4733+06 

2.3739+03 

3.3151-05 

8,0566-06 

-2.1966-03 

6.5464-02 

-1.0827+00 
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STATION 


UE 

ME 

MUE 

CP 

X 

tketaw 

TE 

PE 

MDOTW 

nuw 

S 

OS 

HE 
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TW 


13.0 

1.7226+00 

3.773«*02 

7.2654-02 

2.1172-06 

3.1289-03 

-3.5494*00 

0.0000 

6.6479*03 

4.4463*05 

7.3854-03 

4.1993+00 

2.2501-03 

3.5044-04 

-7.9688*06 

1.8848-02 

1.4000*03 

4.3558*00 

14.0 

1.7228*00 

3.7740*02 

7.2659-02 

2.1172-06 

3.1631-03 

-3 . 5489*00 

0.0000 

6.6479*03 

4,4464*05 

7.3854-03 

4.2458*00 

2.6531-03 

4.0301-04 

-7.9688*06 

1.8848-02 

1.4000*03 

5,1280*00 

15.0 

1.7228+00 

3.7743*02 

7.2665-02 

2.1172-06 

3.1285-03 

-3.5485*00 

0.0000 

6.6479*03 

4.4464*05 

7.3855-03 

4.2000+00 

3. 1165-03 

4.6346-04 

-7.9638+06 

1.8848-02 

1.4000*03 

5.3339*00 

16.0 

1.7226*00 

3.7747*02 

7.2672-02 

2.1172-06 

3.1*0*-03 

3.5480*00 

0.0000 

6.6479*03 

4.4464+05 

7.3855-03 

4.2041*00 

3.6495-03 

5.3298-04 

-7.9688+06 

1.8848-02 

1.4000*03 

4.8701*00 

'.7.0 

1.7228*00 

3.7751*02 

7.2679-02 

2.1172-06 

3.0889-03 

-3.5473 >-00 

0.0000 

6.6479*03 

4.4464+05 

7.3855-03 

4.1486*00 

4.2624-03 

6.1292-04 

-7.9688*06 

1.8848-02 

1.4000*03 

4.4332*00 

18.0 

1.7228*00 

3.7755*02 

7.2688-02 

2.1172-06 

3.0755-03 

-3.5466*00 

0.0000 

6.6479*03 

4.4464*05 

7.3855-03 

4.1316*00 

4.9673-03 

7.0486-04 

-7.9689*06 

1,8848-02 

1.4000**3 

4.1262*00 

19.0 

1.7228+00 

3.7761*02 

7.2698-02 

2.1172-56 

3.0310-03 

-3.5458*00 

0.0000 

6.6479+03 

4.4464+05 

7.3655-03 

4,0729*00 

5.7779-03 

8.1059-04 

-7.9689*06 

1.8849-02 

1.4008*03 

4.1953*00 

20.0 

1-7228+00 

3.7767*02 

7.2710-02 

2.1172-06 

3.0926-03 

-3. 5-49+00 

0.0000 

6.6479*03 

4.4464*05 

7.3855-03 

4,0361-00 

6.7100-03 

9.3213-04 

-7.9689+06 

1.8849-02 

1.4000*03 

4,3061*00 


THRUST LOSS * -7.0205680*02 



i* ^ r 



ST 

OELTA* 

2ETA 

RrfPB 

QW 

TM=TA 

zetap 

XPR 

SUMQW 

ADCSAS 

DUEDS 


3.0075-02 

1.3022-04 

X. 3716-03 

5.2536-01 

1.4731+06 

3,3424-05 

1.3457-03 

-1.0*26*00 

2.8939*03 

9.5426-06 

6.5461-02 


1.9299-02 

1.3092-04 

J .3756-03 

5.2536-01 

1.4146*06 

3,3772-05 

9.9046-03 

-1.0825*00 

3.4799*03 

1,1252-05 

6.5464- 2 


1.5250-02 

1.3161-04 

1 . 3634-03 

5.2536-01 

1.3903*0* 

3,4163-05 

1.6836-02 

-1.0*24*00 

4.1345*03 

1,3216-05 

6.5462-02 


1.1917-02 

1.3200-04 

1.3935-03 

5.2536-01 

1.3158*06 

3,4627-05 

1.8930-02 

-1,0*22+00 

4.8608*03 

1.5476-05 

6.5464-02 


1.0363-02 

1.3383-04 

1.4u >9-03 

5,2536-01 

1.2881+06 

3.5107-05 

1.6885-02 

-1.0820+00 

5.6645+03 

1.6076-05 

6.5*65-02 


6.6451-03 

1,3495-04 

1.4135-03 

5.2536-01 

1.2112*06 

3,5634-05 

1.3640-02 

-1.08l8*00 

6.5516+03 

2,1067-05 

6.5462-02 


7.8335-03 

1.3625-04 

1.4231-03 

5.2536-01 

1.1834+06 

3,6246-05 

1.1874-02 

-1.0816*00 

7.5290*03 

2,4505-05 

6.5462-02 


6.7214-03 

1.3784-04 

1.4345-03 

5.2536-01 

1.1010*06 

3,6961-05 

1. 2241-01 

-1.0»13*00 

8.6012+03 

2.8456-05 

6.5463-02 
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ADDITIONAL PROGRAM FEATURES 


This section describes several features of the Mass Addition Boundary 
Layer Proqram not discussed elsewhere in this report. 

RESTART CAPABILITY 

The Mass Addition Boundary Layer Program contains a restart cap j~i y 
which can save a significant amount of computer time when rerunning cases with 
altered wall and edge tables, different stepsizes or convergence criteria, etc. 

If the case in question yields satisfactory results up to a certain XLIM-station, 
the case may be restarted from that station with the appropriate changes to the 
input data. Implementation of the restart capability is described below. Note, 
however, that this feature has not been adequately checked out for all types of 
cases; therefore, results of a restarted case should be analyzed carefully. 

Writing a Restart Tape 

Logical unit 16 in the Mass Addition Boundary Layer Program is the unit which 
will contain the restart data and must be assigned to a save tape for a restart case. 
Saving of the necessary data at each XLIM-station is initiated by including the 
following variable in the DATA namelist: 

IRSWR Flag (integer) 

=0 (nominal) if no restart data is to be saved; 

=1 if restart data is to be written on logical 
unit 16 at each x/x^ value in the input 

array XLIM. 

Restarting a Case From a Previously Saved Tape 

To restart a case using previously saved data, all data from the original 
case must be reinput, with the following necessary changes in the DATA namelist: 

IRSWR Should oe set to zero if the restart data used 

to restart this ca e is to be preserved; otherwise, 
if IRSWR=1 , restart data at every XLIM -value 
encountered in this case will be written on the 
tape over the previous data. 
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XINIT 


SIN IT 


= the value cf the station counter at the initial 
station for this restart case. This value is used 
to locate the appropriate XLIM -station on the 
restart tape. 

Initial value of x/x N for this case, the same as 
the value in the array XLIM corresponding to 
station IRSRD above. 

Initial value of sT/x^ for this case. 


DXI Initial stepsize ix/L for this case, the same as 

the value in the array DXLIM corresponding to 
station IRSRD above. 


SKTAB \ Must be changes such that XTABSK(l) is greater 

XTABSK > than XINIT above. 

LSKTABj 


DXLIM' 
XLIM » 
LDXIM 


Must be changed such that XLIM(l) is greater 
than XINIT above. 


If station IRSRD is not present on the restart tape, the following message is printed: 

THERE IS NO RESTART DATA FOR STATION [iRSRD] . 

If the restart case is set up to duplicate the results of the original case 
for the first ten stations or so, the summary table printouts for these stations may 
be compared for possible discrepancies. Small errors in the restart case should, 
in general, damp out. 


DIFFERENCE EQUATION DEBUG CAPABILITY 


The ability to dump the difference equation coefficients and the current 
values of the dependent variables u, H, h, and pv has been provided (See SUBROUTINES 
DUMPIT for further details.) and is implemented by entering values far the arr$y 1 

IDEBUG in the DATA namelist as follows: 


IDEBUG(l) 


Flag (integer) 

=1 if the momentum equation coefficients are 
to be printed; 

-2 if the energy equation coefficients are to be 
printed. 
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IDEBUG(2) 
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Integer specifying the first station for which 
the coefficients and variables are to be printed. 

IDEBUG(3) Integer specifying the last station for which 

the coefficients and variables are to be printed. 

Coefficients and values for the f ir^t five and the last five mesh points are printed 
for every iteration over the range of stations specified by IDEBUG(2) and IDEBUG(3) . 

VARIABLE LEWIS NUMBERS 


The Lewis number Le and the turbulent Lewis number Le,p at each mesh 
point are currently set equal to ALEWIS and TLEWIS, where the latter are program 
constants set to unity in the main program TFCBL. However, back, average and 
forward values of Le and Le^, for each mesh point are stored in the program arrays 
BLE and BLET, so that a method of calculating variable Lewis numbers could be 
provided without difficulty. If this were done, the following areas of the current 
code would have to be changed accordingly: 

(1) D0 40 loop in main program TFCBL, which sets Le=TALEWIS across 
the boundary layer. 

(2) DQf 600 loop in subroutine EDDY, which sets Le^,=TLEWIS across 
the boundary layer. 

(3) D0 4300 loop in subroutine H00DE, which sets Le=l across the 
boundary layer for an equilibrium case. 

(4) The D0 4910 loop at the end of subroutine H00DE, which is 
currently bypassed, must be executed fc .e^l; the c^'s and h^'s 
must be obtained by interpolation for mesh points intermediate to the 
equilibrium calculations, since the terms containing them in the 
difference equations are now nonzero. 
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ADDITIONAL NOTES ON INPUT 


REFERENCE QUANTITIES 

The reference quantities are used to nondimen sionalize the dependent and 
independent variables. By using reference values as close as possible to the 
values of each variable encountered in the solution, *he normalized variables be- 
come more uniform in magnitude (all of order unity in an ideal case) and the possibility 
of round-off error is reduced. As long as significant round-off error is not encountered 
the calculated results are independent of the choice of reference values. 

PROBLEM LIMITS AND INITIAL VALUES 

DXI - It is desirable to start the solution with a series of small steps in order 

to damp out the initial truncation error inherent in starting from initial profiles that 

do not satisfy the differential equations. It appears that 30, or so, steps of the 
-5 -4 

order of 10 ' to 10 feet is usually sufficient. The acutal number of small initial 
steps or their size can be varied over a considerable range without significantly 
affecting the results. DXI itself is the value of the initial step size. These initial 
steps can all be the same size or the step size can be slowly increased in a geometric 
progression. The controls for changing the step size are discussed later on. 

DELTAI - See discussion of PLAW in constants section. 

2ETAPI - The initial value of is normally not known and therefore the 
program has been designed to calculate one. The results are not sensitive to the 
value of 2ETAPI as long as the solution is begun with a series of small steps. 

CONSTANTS 

PLAW (and DELTAI) - If one is starting from known initial conditions 

(experimental cr theoretical) DELTAI is known and the known initial profiles can 

be input using the UPROF and HPROF arrays. In this case, PLAW is not needed. 

Normally, if initial conditions are not known, the solution should be started as close 

as is practicable to the origin of the boundary layer flow in order to reduce the overall 

effect of initial profile errors. If the boundary layer is started near its origin (DELTAI 
-4 -5 

«10 or 10 depending on Reynolds number) then the effect of initial profile errors 
is quickly damped out as the boundary layer grows. Under such conditions the choice 
of PLAW is not important. 
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If, for some reason, a calculation must be started at a location where the 
boundary layer has appreciable thickness and the thickness and boundary layer 
profiles are not known, then the choice of DELTAI and PLAW can have a significant 
effect on the calculated results. Under such conditions one can do nothing more 
than use best estimates for DELTAI and PLAW. 

GPO,SN3 - In most of the calculations performed to date, values of GPO 
of about 1000 a. d Sn3 = 4 have been used successfully. If very high Reynolds 
numbers are encountered larger values of GPO may be needed. GPO and SN3 have 
not been systematically varied to date so that a more detailed set of recommendations 
cannot be made at this time. 

CONVERGENCE AND EDGE CRITERIA 

EPSLN1, EPSLN2, EPSLN3 - It is suggested that values of the c's in the range 
-1 -2 

10 to 10 be used. The solution is not sensitive to the values of the c's. Normally 
when the solution terminates for adding too many points either the problem has been 
poorly posed (e.g. an input error) or the step sizes require modification. On rare 
occasions the program adds points extraneously due to a local numerical problem. If 
this should happen, increasing the value of c can usually circumvent the problem. 

COUNTERS 

MAXIT - In the majority of cases MAXIT should be set equal to 1. Occasionally 
the oscillations in r w and q w do not damp out (or amplify) with MAXIT =1. If that is 
the case increasing MAXIT to 2 usually solves the problem. 

NYI - The initial number of mesh points in the Y direction for turbulent 
boundary layer flows has been varied from 80 -* 160 in the solutions obtained to 
date. Most of the computations have been carried out with NYI = 120. The exact 
number of mesh points to use is not critical since the computed results are usually 
fairly insensitive to the mesh spacing over a considerable range. The few laminar 
solutions which have been obtained used 60 or less mesh points (usually with no 
stretching - GPO and SN3 not input) . 
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IYEQ - For ^ 2~®2 equilibrium solutions a large part of the computer time is 
used in solving for the equilibrium concentrations of the species. To reduce the 
amount of computer time required to solve equilibrium boundary layer flows 

the species concentrations are solved for only at every IYEQth mesh point and the 
values at other mesh points are found by interpolation. For laminar flow solutions 
appreciable error can be encountered if IYEQ is much larger than 4. Based on 
current experience it appears that turbulent flows can be computed with IYEQ = 8 
without introducing more than 1% errors. 

STEP SIZE CONTROL TABLES 

The XLIM, DXLIM arrays are used to effect discontinuous changes in step 
size. Sudden large changes in step size do not usually cause problems as long as 
they are not encountered too frequently. (Equilibrium solutions can be more sensitive 
to step size changes) . A large step size change can cause a local perturbation in 
the solution which normally quickly damps out. If the step size is repeatedly changed 
by a significant factor, before the perturbations from the previous changes have 
disappeared, then erroneous results can be computed. 

The SKTAB, XTABSK arrays allow continuous step size changes to be made. 

The step sizes form a geome.ric progression of common ratio SKTAB(K). If SKTAB 
is specified as one over a given X range, the step size will remain constant. SKTAB' s 
less than one cause step size reductions and SKTAB' s greater than one increase the 
step size. If the last value of X computed is denoted X Q , the latest step size is £X q 
and SKTAB = K the following formulas can be used to compute the step size and distance 
covered after n additional steps 

AX = aX K n 
n ° o 

K** - 1 

X = X + AX^ f 

n o o K - i 

The two methods of varying the step size can be used separately or together 
to suit the requirements of a particular problem. A large number of problems have 
not yet been solved with the present program so only general guidelines , rather than 
exact step-size recipes, can be currently suggested. The size of the X steps to be 
used should be governed by the size of the gradients in X direction, e.g. , dr w /dx, 
dT /dx, dP /dx, dm/dx, etc. The higher the gradients the smaller the step size 

W w 

required to maintain a given level of truncation error. Many times discontinuities in 
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one or several of the boundary conditions are encountered. Disco.. tinuities should 
be handled essentially as if the solution were being started with approximate initial 
profiles. An XLIM value should be set equal to the value of X at which the dis- 
continuity begins, the corresponding DXLIM entry should be used to reduce the step 
size to a small value (say « 10 - ^ of the total length to be computed) . This small 
step should then be maintained (or slowly increased by a K factor) for 30 or so steps 
until the perturbing effect of the discontinuity has been damped. 

In regions where rapid gradients or discontinuities (including the initial 
"smoothing out" region) are not present, a step size on the order of 1/100 of the 
length to be computed has typically been used. Much smaller steps have been used 
without changing the results by more than 1 or 2%, however, the effect of larger 
steps has not been adequately determined. In regions with sizeable gradients smaller 
(than 1/100 the length) steps should be used. Experience will be the only guide as 
to how small a step one need take in order to achieve the desired accuracy. 
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